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Executive Summary 

This report documents an experimental and analytical study of the active stabilization of 
surge in a centrifugal engine. The aims of the research were both to extend the operating range of a 
compressor as far as possible and to establish the theoretical framework for the active stabilization 
of surge from both an aerodynamic stability and a control theoretic perspective. In particular, 
much attention was paid to understanding the physical limitations to active stabilization and how 
they are influenced by control system design parameters. 

Previously developed linear models of actively stabilized compressors were extended to 
include such nonlinear phenomena as bounded actuation, bandwidth limits, and robustness criteria. 
This model was then used to systematically quantify the influence of sensor-actuator selection on 
system performance. Five different actuation schemes were considered, along with four different 
sensors. Sensor-actuator choice was shown to have a profound effect on the performance of the 
stabilized compressor. The optimum choice was not unique but rather shown to be a strong 
function of some of the non-dimensional parameters which characterize the compression system 
dynamics. Specifically, the utility of the concepts were shown to depend on the system 
compliance to inertia ratio ("B” parameter) and the local slope of the compressor speedline. In 
general, the most effective arrangements are ones in which the actuator is most closely coupled to 
the compressor, such as a close-coupled bleed valve inlet jet, rather than elsewhere in the flow 
train, such as a fuel flow modulator. 

The analytical model was used to explore the influence of control system bandwidth on 
control effectiveness. The relevant reference frequency was shown to be the compression 
system’s Helmholtz frequency rather than the surge frequency. The analysis shows that control 
bandwidths of three to ten times the Helmholtz frequency are required for larger increases in the 
compressor flow range. This has important implications for implementation in gas turbine 
engines since the Helmholtz frequencies can be over 100 Hz, making actuator design extremely 
challenging. 

One result of the modelling was to explain the power flow in an actively stabilized 



compressor. In particular, it was shown that the power the control system must input (or extract) 
scales not with the power of the compressor but with the power of the perturbations in the system, 
which are generally 10 4 to 10 6 smaller. The analysis indicated that the source of the perturbations 
are important. For example, three times more actuation authority is required if the destabilizing 
fluctuations come from perturbations in heat release in the combustor than if they come from 
perturbations in combustor outflow (turbine flow), and 6 times more if the source is unsteady 
pressure rise in the compressor. The actual structure of the perturbations in gas turbine engines is 
completely unknown. Since it is only important in the context of active control, it had not been 

previously studied. 

The experimental investigation centered on a small (40 KW) turbocharger configured in a 
pumping system resembling that of a gas turbine engine, including compressor, plenum 
(combustor), and throttle (turbine). The actuation consisted of a high speed throttle valve which 
could be positioned either at the plenum exit or between the compressor and plenum. Sensors 
included mass flow, total and static pressure at the compressor inlet, and pressure in the plenum. 

Open-loop actuation of the valve permitted measurement of the transfer function of the 
compression system. The flow development process was modelled as a one-dimensional lag. 

With this addition, a lumped parameter model was then shown to model the measured compressor 

response. 

Closed-loop control of surge was shown to decrease the mass flow at surge inception by 
50% compared to the uncontrolled case. The modelling was shown to do quite a good job in 

predicting the stability boundaries of the stabilized compressor. 

Overall, the results of this research are very encouraging. Active control of surge has been 
shown to give larger gains when applied to centrifugal compressors. The performance of the 
controlled compressor has been shown to be well predicted by theory, indicating that we now have 
the tools to extend this concept to new geometries and applications. Work in this area could be 
very profitably pursued in applying these concepts to small gas turbine engines with centrifugal 
compressors. Much would be learned from such research. 
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Chapter 1* 


Introduction 

1.1 Introduction 

The useful range of operation of many turbomachine type compressors is limited 
by the onset of a system instability, known as surge, which occurs when the flow 
is reduced below some critical value. The unstable operation is characterized by 
periodic variations in the net through flow and output pressure of the compressor. 
The severity of the oscillations depends upon the overall dynamics of the system 
in which the compressor is installed and can be quite energetic, with flow reversal 
through the system. The forces involved may be great enough to cause structural 
damage [29] and, in addition, there may be consequences such as loss of power and 

overheating in an aircraft gas turbine engine. 

Traditionally, such undesirable behavior has been avoided by using control systems 
which prevent the operating point from entering the unstable regime delineated by 
the surge line ([32, 47]). The approach described here is fundamentally different than 

* The remainder of this document is condensed from the Ph.D. thesis of J.S. Simon. 


19 



these previous concepts in that it makes use of feedback stabilization to allow stable 
operation in the naturally unstable regions of the performance map [10]. This idea 
is illustrated in Figure 1-1 which shows an extended region of stable operation made 
possible with feedback stabilization. 

Such an approach may have significant benefits for particular applications. The 
peak efficiency region is often quite close to the surge line [48] and feedback stabiliza- 
tion could thus allow stable operation in a more efficient region of the performance 
map. Other benefits include greater operating range and greater flexibility in design- 

ing compressors, if stability is not a primary concern. 

Several researchers to date [12, 20, 40] have experimentally demonstrated that 
significant extensions in mass flow range can be achieved using various feedback sta- 
bilization mechanisms. Those studies focused on modeling and demonstrating partic- 
ular implementations. In contrast, the overall purpose of this thesis is to report the 
results of a fundamental study of some of the more generic aspects of the compression 
system feedback stabilization problem. 

It should be noted that any such attempt to artificially stabilize a compression 
system is not without associated risks and costs which must ultimately be considered 
and weighed against the benefits. In particular, the consequences of a control system 
failure could be be catastrophic in an aircraft application. The expense of additional 
hardware, added complexity and a probable resulting reduction in overall reliability 
are all admitted drawbacks. Before any final judgement can be passed, however, it 
is necessary to obtain a thorough understanding of the basic physics involved, the 
fundamental limitations to stabilization, and the available implementation options. 
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This thesis represents a first step in this direction. 


1.2 Previous Work 

i 

1.2.1 Modeling of Compression Systems 

Low order, lumped element, models for the surge process have been proposed by 
many researchers. Classic papers in which such models are derived and analyzed 
include those of Emmons et al [9] and Taylor [48]. Other references in which such 
models may be found include [2, 17, 22, 30, 49]. A comprehensive treatment is given 
by Greitzer [19] including an extensive literature review. 

The approaches taken by these authors have many similarities and include assump- 
tions of one dimensional, incompressible flow in the ducts, spatially uniform pressure 
in plenum chambers and quasi-steady pressure/flow relations for the throttling valves. 
The compressor is treated as an actuator disk which either follows its steady state 
performance characteristic (the quasi-steady assumption) or, in some cases (see for 
example Greitzer [17]), a first order lag is used to model the flow development process. 

The system model is obtained by performing momentum balances for the ducts, 
and mass balances for the plenum chambers with the ideal gas law and some form of 
polytropic process used to relate the plenum density to plenum pressure. The models 
proposed by these authors vary primarily in the number of elements which are used, 
with some including additional plenum chambers and/or ducts at the compressor inlet 

or exit. 

All the authors provide linearized stability analyses but some examine some of the 
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nonlinear aspects of the problem as well. Most authors recognize the importance of the 
slope of the compressor pressure-flow characteristic at the equilibrium operating point 
as a prime determinant of the compression system stability. Greitzer [17] also found 
global conditions on the compressor slope which must be satisfied to maintain surge 
oscillations. In addition, Greitzer showed the importance of the relative amount of 
fluid compliance to inertia (which he captured in his nondimensional B parameter 1 .). 
More recent nonlinear analyses include those of Oliva and Nett [36] who looked at 
Hopf bifurcation phenomena and those of Simon and Valavani [45] who performed a 
Lyapunov stability analysis of a simple compression system. 

The validity of these models has also been examined by a number of authors 
who all show that this simple modeling approach can capture both the qualitative 
and quantitative aspects of the surge phenomena for a broad class of turbomachines. 
Using a centrifugal compressor, Emmons [9] found reasonable agreement between 
linear theory for both the experimentally predicted onset of instability (the neutral 
stability point occurred near the peak of the compressor characteristic) and for the 
dependence of the system instability frequency on the system parameters. Greitzer 
[18] performed tests using a low speed axial flow fan and found good agreement for 
the predicted nonlinear behavior when the system exhibited large amplitude surge 
oscillations. In particular, he verified the theoretical prediction of the importance of 
the B parameter in determining whether a system would experience sustained surge 

J The B parameter is defined by B = where: is the tip speed, a p is the speed 

of sound in the plenum, V p is the volume of the plenum, C e is the equivalent inertial length of the 
compressor duct, and A c is the compressor inlet area. 
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oscillations in annulus averaged mass flow or operate steadily (without surging) in a 
condition known as rotating stall. He showed good agreement between the predicted 
and actual time histories of the system when it was experiencing sustained oscillations. 
Hansen [21], using a small centrifugal compressor, also showed that a model similar to 
Greitzer’s was able to accurately predict the nonlinear surge behavior. Fink [13, 14] 
performed an in depth study of surge, also using a centrifugal compressor, and showed 
impressive agreement between predicted and measured behavior. Finally, Bons [5] 
found a good comparison between experiment and theory using the same modeling 
approach for a centrifugal pump, thus demonstrating that the usefulness of the model 

is not limited to gaseous compression systems. 

The models proposed by the above authors are all lumped parameter models, 
that is the system is represented by a finite number of ordinary differential equations. 
A distributed parameter model of a fan-duct-plenum system, forced by modulating 
dampers, has been developed and validated by Goldschmied and Wormley [16]. Their 
model uses a one-dimensional lossless acoustic transmission line representation for 
the duct, rather than a lumped fluid inertia. This is a more appropriate model 
for their system in which the acoustic wavelength at the higher forcing frequencies 
was less than the length of their duct. The fan and plenum were modeled using an 
approach similar to the other authors. The overall agreement between prediction and 

experimental measurements was found to be acceptable. 

The most important general finding which is common to all of these studies is that 
it is not necessary to model the detailed, and very complex, flow phenomena occurring 
in the compressor itself in order to accurately predict the overall system behavior. 
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The pressure-rise versus flow characteristic of the compressor apparently provides a 
sufficient description for predicting the overall system dynamics. This was true even 
though in many of the experiments the flow in the compressor was three dimensional 
and locally unsteady, with rotating stall occurring in the axial flow machines and 
various diffuser stalls occurring in the centrifugal machines. 

1.2.2 Modeling of Compressors 

The previously cited work provides indirect evidence that a quasi-steady or simple 
first order lag model for the compressor is adequate. There is also some direct ex- 
perimental support for this modeling approach provided by researchers who have 
attempted to measure the transfer function of the compressor. Ohashi [35] performed 
both an analytical and experimental study of the frequency response of turbo-pumps. 
He made a careful distinction between the unsteady behavior which is due only to 
fluid inertia effects and those due inherently to the unsteady pump performance. 
He found experimentally that the pump performance deviated significantly from the 
quasi-steady predictions when the reduced frequency (based upon convection time 
through the blade row) of the pump was greater than 0.1. Above this frequency the 
unsteady response of the compressor could be only roughly approximated by a first 
order lag. Abdel-Hamid [1] measured the frequency response of a low speed, single 
stage centrifugal blower with a tip diameter of 660mm operating in air at a tip Mach 
number of 0.36. His basic result was a polar plot showing the real and imaginary parts 
(in-phase and quadrature components) of the sinusoidal transfer function relating flow 
perturbations to pressure perturbations (This is the impedance of the compressor al- 
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though he did not call it that.) Over the range of frequencies and flows which he 
studied, the real part of the impedance varied only slightly while the imaginary part 
increased strongly with forcing frequency. He attributed this to a fluid inertia ef- 
fect. However, the equivalent compressor length required to explain this data was a 
function of the mean flow rate, which does not seem physically plausible given that 
equivalent compressor length is a purely geometric quantity. The variation m the real 
part of the transfer function also would not be expected from quasi-steady theory 
which predicts that the real part remain constant at the value of the quasi-steady 
slope corresponding to the mean flow. Another attempt to characterize the unsteady 
response of a compressor was reported by Paulon [38] who studied a low speed single 
stage axial flow compressor, but no clear conclusions were reached. 

Review of the available literature thus indicated the need for the direct assess- 
ment of the validity of simple lumped element models of the compressor and such 
an assessment was performed as part of the current research effort described in this 

thesis. 

1.2.3 Feedback Stabilization of Compression Systems 

There have been a number of demonstrations that the stable operating range of a 
compressor can be extended by modifying the overall dynamics of the system. Such 
modification can always be interpreted as providing a stabilizing feedback mechanism. 
In terms of implementation, these methods can be divided into two general categories: 
unilaterally coupled and bilaterally coupled (see for example [53]). 

In a unilaterally coupled feedback system, power flows essentially one way; from 
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the actuator into the controlled plant. A typical example of a unilaterally coupled 
feedback consists of a silicon strain gauge pressure transducer measuring pressure in 
a plenum which sends an electrical signal to a computer which drives (through an 
amplifier) a speaker connected to the plenum. In contrast, power flows both ways in 
a bilaterally coupled feedback system. For example, if the transducer, computer and 
speaker of the previous example were replaced by a spring loaded, movable plenum 
wall, a bilaterally coupled feedback mechanism would be obtained. These two types 
of implementations are also commonly distinguished as active and passive 2 . 

The idea of stabilizing a compressor using unilateral feedback was proposed by 
Epstein et al [11]. An initial theoretical investigation of surge stabilization was con- 
ducted by Chen [6] who examined the stabilization of compression system using a 
movable plenum wall and a variable throttle area at the plenum exit using linear 
stability theory. Chen also investigated the use of nonlinear sliding mode control 
with a plenum exit control valve and concluded that the approach was promising. 
A reference available in Russian [22] also discusses application of unilateral feedback 
stabilization of compression systems. 

Two experimental investigations ( Pinsley et al [40] and Ffowcs Williams and 
Huang [12]) demonstrate that a unilateral feedback control can be used to allow a 
compressor to operate steadily at naturally unstable points on the compressor map. 
Both control schemes were based upon a linearization of a lumped element model 


2 The terms active and passive already have useful precise definitions in the electrical circuit 
literature (See for example Wyatt et al [52].) which do not necessarily coincide with this common 
usage. Although unilaterally coupled feedback systems are typically active and bilaterally system 
are typically passive (in the electrical circuit sense) this correspondence need not hold in general. 
Thus, to avoid confusion, this terminology will not be used here. 


26 


similar to the type described in the previous section. Ffowcs Williams and Huang 
used a measurement of plenum pressure to provide a feedback signal which was used to 
drive a speaker mounted on the plenum wall, effectively changing the plenum volume 
in response to perturbations in plenum pressure. Pinsley also used measurements of 
plenum pressure to provide a feedback signal, but used it to drive a servo- actuated 
valve at the plenum exit. In both cases, the feedback control law (compensator) was 
implemented using a linear, frequency dependent, analog network. 

Both investigators were successful at moving the surge line appreciably in some 
regions of the performance map but not in others; in particular, little effect was 
found at higher wheel speeds. Both researchers indicate that their inability to obtain 
control at higher speeds was attributable to diminished effectiveness of their actuated 
variable, although no detailed analysis was performed. Neither researcher considered 
the dynamics of the feedback components except to study the effect of changing the 
compensator gain and phase at a single frequency (the open loop system natural 

frequency). 

There have been basically two types of bilateral implementations: close-coupled 
resistances and tuned absorbers. Dussourd [8] showed experimentally that adding a 
flow resistance immediately downstream of the compressor exit extended the surge 
free mass flow range. He found that, for compressors whose pressure-rise versus 
flow characteristics were not too steeply positively sloped, substantial increases in 
surge free operating range could be achieved with a relatively small pressure drop 
penalty. Further, the stability limit of the overall system could be well predicted 
based upon the simple criterion that the slope of the combined pressure-rise versus 
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flow characteristic of the compressor and valve be zero at instability. (This would 
be expected from the simple lumped element model for systems with a reasonably 
large value of B parameter. ) It should be noted that this result was apparently well 
known, at least empirically, as long ago as the 1930’s with commercially available 
devices providing throttling at the compressor inlet as described by Kearton [23]. 

At least two tuned resonator approaches have been reported in the literature. 
Bodine [4] used a no moving parts Helmholtz resonator (a plenum chamber with a 
short neck) mounted at the compressor exit to extend the stable flow range but did 
not provide any detailed analysis. Gysling [20] used a mechanical resonator to extend 
the stable flow range. In Gysling’s purely aero- mechanical scheme, the compression 
system was coupled to a mass-spring-dashpot arrangement through a movable plenum 
wall. Functionally, plenum pressure was fed back to produce changes in plenum 

volume. 

Gysling modeled both the dynamics of the mass-spring-dash pot system and the 
compression system using simple lumped element models and showed that an ap- 
preciable surge line shift (20-30% in mass flow) could be obtained by optimizing the 
choice of system parameters. Gysling’s design was based upon a linearization of the 
system dynamic model which was then checked using a nonlinear simulation. Good 
agreement was shown between the experimental and theoretical range of stabilization. 
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1.3 Current Research Objectives 

The specific objectives identified for the research reported in this thesis are classified 
into two main categories: modeling and model validation related objectives and those 
concerned directly with feedback stabilization. In this section, the identified objec- 
tives in both of these categories are first motivated, and then are briefly summarized. 

1.3.1 Modeling and Model Validation 

Although an established approach for modeling compression system surge instabilities 
is available, two research goals involving modeling related work would be beneficial 
were identified as useful to pursue. The first was to develop control oriented models 
which would include a variety of attractive actuation and sensing schemes. Many of 
the currently available models are oriented toward studying the autonomous system 
dynamics and do not include any actuators and sensors. Development of models in 
a form readily amenable to performing control system design and analysis was thus 
a necessary step for the subsequent thesis work described here. Availability of such 
models also provides the benefit of making this interesting problem accessible to other 

researchers. 

The second modeling related research goal was to experimentally assess in de- 
tail, the widely used approach of modeling the compressor as a quasi-steady actuator 
disk coupled to a one-dimensional incompressible duct. Although indirect evidence, 
based upon studies of overall system dynamics, exists to support this modeling ap- 
proach, there have been few unambiguous studies of the compressor alone. Since it is 
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the compressor behavior which drives the surge instability, accurate modeling of the 
compressor is an item of key importance. Thus, the results of such an experimental 
validation would provide an important contribution. 


1.3.2 Feedback Stabilization 

Work to date on controlling compression system instabilities has focused on analyz- 
ing and or implementing particular control schemes, and little of the available, and 
apparently applicable, control theory has been used to study this problem. Several 

research objectives were accordingly identified. 

The first was to establish and quantify the fundamental theoretical limits to sta- 
bilizing a compression system. The second objective was to evaluate a variety of 
potentially attractive control schemes and to identify the most promising generic 
stabilization approach to pursue. The third objective was to obtain a detailed un- 
derstanding of the most promising concepts and to experimentally demonstrate the 

viability of this option. 

1.3.3 Summary 

The basic objectives of this research are briefly summarized as follows: 

• Assess experimentally the quasi-steady compressor model. 

• Develop control oriented system dynamic models of the compression system for 
the more attractive actuation and sensing options 
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• Characterize and establish quantitative limits to feedback stabilization of com- 
pression systems 

• Determine the most promising generic approach to stabilizing compression sys- 
tems. 

• Develop a detailed understanding of the most promising stabilization approach 
and experimentally demonstrate its viability 

1.4 Overview of Thesis 

To accomplish these research objectives a combined analytical and experimental ap- 
proach was taken. This section provides a brief overview of the overall thesis orgam- 
zation. 

• Chapter 1 - Introduction and Background 

The general problem of compressor instability and motivation for stabilizing 
compression systems is introduced. Potential benefits and risks are discussed. 
Previous work is reviewed to establish the starting point for this research. Re- 
search objectives are outlined. 

• Chapter 2 - Experimental Facility 

This chapter provides an overall description of the experimental facility utilized 
for experiments described in the thesis. 

• Chapter 3 - Modeling of Actuation and Sensing 

This chapter develops the models necessary for performing control system design 
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or analysis for the compression system using the various actuation and sensing 
schemes which were considered. Inclusion of fuel modulation and mass injection 
amongst the set of potential actuators requires some extensions to the simple 
model presented initially and these are made here and explained. For example, 
the isentropic plenum assumption must be dropped introducing an additional 
states. Models of actuation and sensing previously developed must also be 
modified to be consistent with the new assumptions and this is also done here. 
The models are put into a form readily amenable to control system design and 
analysis. 

• Chapter 4 - Compressor Model Validation 

Modeling of the unsteady behavior of the compressor itself is the most funda- 
mental aspect of an overall compression system model. To directly assess our 
ability to model this behavior, the frequency response of a representative cen- 
trifugal compressor was determined. This chapter presents both the experimen- 
tal technique and data analysis methods which were developed for this purpose 
and the comparison of experimental results with theoretical predictions. 

• Chapter 5 - Limitations to Stabilization of Compression Systems 

This chapter establishes and quantifies some fundamental limitations to stabi- 
lization of compression systems. First it is established that achieving nominal 
linear stability is not generally problematic; with most choices of actuator and 
sensor the system is controllable and observable so the poles may be placed 
arbitrarily in the left half plane. This point had not previously been clearly es- 
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tablished. With this issue clarified, attention is directed at what actually limits 
our ability to control such systems. This is the topic of the second part of the 
chapter where more fundamental limits to stabilization imposed by practical 
considerations such as bounded actuation and model uncertainty are character- 
ized and quantified. 

• Chapter 6 - Evaluation of Alternative Stabilization Strategies 

In this chapter a side by side quantitative assessment of a representative assort- 
ment of potential pairings of actuators and sensors is performed. The results 
presented here serve two main purposes. First, they provide a rational basis 
for selecting the more viable options for more intense scrutiny. Second, they 
establish, in a concrete manner, the real impact that choice of actuator and 
sensor location has. This fact, particularly regarding sensor selection had not 
been recognized previously. 

• Chapter 7 - Close-Coupled Control Theory and Design 

The use of an actuator and sensor close-coupled to the compressor can have sig- 
nificant advantages as could be qualitatively anticipated from the model physics 
and quantitatively established in Chapter 6. In this chapter, the rational de- 
sign of a particular implementation (a valve close-coupled to the compressor exit 
with mass flow measurement) is analyzed in detail and experimentally evaluated 
for the first time. A nonlinear, Lyapunov stability analysis of this close-coupled 
control strategy is also provided in Appendix E. This analysis is useful when 
considering the effect of large disturbances acting on the system. 
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Chapter 8 - Summary and Conclusions 

In this final chapter the overall research effort is summarized, a final set of 

r 

conclusions conclusions are provided, and suggestions for further research are 
outlined. 


34 



Figure 1-1: Extended range of surge free operation provided by feedback stabilization. 
(Data of Gysling [20], and Pinsley, [38]) 
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Chapter 2 


Experimental Facility 


A turbocharger consisting of a single stage centrifugal compressor driven by a radial 
inflow turbine was used for the experimental phases of the present work. An overview 
of the experimental facility is provided by Figure 2-1 which shows schematically the 
arrangement of the various subsystems, which are listed below. 

Compression System The compression system consists of the compressor itself, a 
plenum chamber and connecting ductwork. 

Actuation System The actuation system is made up of a valve which controls flow 
through the compression system, along with a high bandwidth D-C servo motor, 
servo-amplifier and motor position controller. 

Measurement System The measurement system consists of the various transduc- 
ers (pressure, temperature, speed, etc.) along with their power supplies, and 
signal conditioning hardware. 
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Data Acquisition and Control System The data acquisition and control system 
consists of a digital computer, A/D (analog to digital converter), D/A (digital 
to analog converter), user interface, mass storage, and software. 

Auxiliary Support System The auxiliary support system includes the drive tur- 
bine, high pressure air supply (used for the turbine drive), oil supply, cooling 
water and cooling air supply for the servo motor. 

The basic facility was originally constructed by Pinsley, and is described in detail in 
[39]. The compression system and auxiliary support systems have not been signif- 
icantly modified and Pinsley [39] can be referred to for full details. The actuation 
system, measurement system, and data acquisition system components have been ei- 
ther extensively modified or replaced for the purpose of performing the experiments 
described here. A full description will therefore be provided in the following sections 

for these systems. 

2.1 Compression System 

The compression system consists of the following components: inlet duct, compressor, 
compressor exit duct, plenum chamber, throttle valve and metering duct. 

2.1.1 Flow Path 

The arrangement of these components and air flow path is shown in Figure 2 2. 
Ambient air enters the system through a bellmouth at the entrance of the inlet duct. 
The air passes through the inlet duct into the compressor which discharges into the 
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inlet duct diameter 

5.2 cm 

annular inlet area 

12.5 cm 2 

hub to tip radius ratio 

0.37 

exit tip diameter 

5.5 cm | 

exit duct diameter 

4.4 cm 


Table 2.1: Compressor Geometry 


compressor exit duct. The air then passes through the compressor exit duct into 
the plenum chamber. The air exits from the plenum chamber through the throttling 
valve into the metering duct. In the metering duct the flow passes through a series 
of perforated plates and finally discharges through a sharp-edged ASME orifice into 
the room. The servo- controlled valve could be located either at the plenum entrance 
or exit as described subsequently. 

2.1.2 Description of Compressor 

A Holset Model HID turbocharger was used for all tests. In this turbocharger the 
compressor and drive turbine are mounted on a common shaft. The compressor is 
a centrifugal type, with a vaneless diffuser, a volute and no inlet guide vanes. The 
impeller has 6 blades and 6 splitter vanes. An exploded view of the compressor 
(reprinted from the Holset Manual) is shown in Figure 2-3. Some specifications of 
the compressor geometry, taken from Pinsley [39], are given in Table 2.1. The com- 
pressor has a peak rated pressure ratio of 3.1 at 140,000 rev/min at a flow rate of 
approximately 0.10 kg/s. The tests described here were done at approximately 90,000 
rev/min, with a pressure ratio of approximately 1.75. An overall performance map 
from [39] is shown in Figure 2-4. 
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2.1.3 Flow Path Components 

The inlet duct is constructed from 2” PVC Schedule 40 piping. Different inlet lengths 
were used and the length of the inlet used for each particular test is specified sub- 
sequently when these tests are described. The plenum is constructed of 10 304 

Stainless Steel Schedule 10S pipe. The plenum volume is adjustable and the actual 
volume used is provided in the description of the individual tests. A l|” Jenkins Ball 
Valve located at the plenum exit was adjusted manually to throttle the flow. The 
metering duct consists of a 1.07 m section of 4 Schedule 10, Type 304 L pipe, baffled 
with 3 perforated plates terminating in a 1.75 inch (4.4 cm) ASME standard sharp 

edged orifice plate. 

2.2 Actuation System 

The actuation system provides the means to impose control actions upon the com- 
pression system. The actuation system consists of the following components: 1) a 
custom made rotary valve, 2) Pacific Scientific Model 4VM81-220-1 Low Inertia, Per- 
manent Magnet D.C. servo motor, 3) Copley Controls Corporation Model 240 D.C. 
servo amplifier 4) Galil Model DMC - 400 - 10 Digital Motion Control Board, and, 
5) BEI Motion Systems Model BEI-MX-213-18-1024 optical shaft encoder. The an- 
gular control valve was constructed previously by G. Guenette and J.E. Pinsley and 
a detailed description is provided in [39]. For completeness, the valve is shown in 
Figure 2-5. The remaining components form a position control loop for the valve. 
This position control loop operates cyclically at a rate of 2000 cycles/sec as follows: 
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1. The digital motion control board reads its input buffer to obtain the commanded 
valve position. 

2. The digital motion control board decodes the shaft encoder signal to obtain the 
shaft position 

3. The digital motion control board generates an error signal by subtracting the 
actual shaft position from the commanded value 

4. The error signal is input to a digital filter which implements a PID (Proportional- 
Integral-Plus-Derivative) control algorithm. 

5. The output of the PID controller is converted to a 0 to 1 volt analog signal 
which is input to the D.C. servo amplifier 

6. The D.C. servo amplifier outputs a pulse width modulated current to the servo 
motor. The time mean value of the amplifier output current is proportional to 
its voltage input. 

7. The current supplied to the D.C. motor produces an accelerating torque which 
rotates the valve towards the desired position. 

The small signal closed loop bandwidth (the frequency at which the ratio of the actual 
valve position to that commanded becomes less than unity with a small sinusoidal 
drive signal) of the valve position control loop is approximately 100 Hz. 
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2.3 Measurement system 


The measurement system is comprised of the transducers, transducer excitation, and 
signal conditioning devices. Both high and low frequency response measurements 
were made to determine the time resolved and steady-state behavior, respectively, of 
the system. Measurements were made of the following quantities: 

1. inlet duct mass flow rate, temperature and wall-static pressure 

2. compressor rotational frequency, exit total pressure and exit static temperature 

3. control valve angular position and total to static pressure-drop (differential pres- 
sure) 

4. plenum wall-static pressure 

5. orifice upstream temperature and differential pressure 

The location of these sensors is shown schematically in Figure 2-6. Each channel of 
data was sampled at a rate of 1 KHz (1000 samples /sec/channel) by an analog to 
digital converter. A full description of the A/D is is given in a subsequent section. To 
prevent aliasing all time resolved measurements were low-pass filtered using 8-pole 
Butterworth Filters with a -3 dB cutoff frequency of 500 Hz. 
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2.3.1 Inlet Duct Instrumentation 

Time Resolved Mass Flow Rate 

The time resolved (instantaneous) mass flow rate was obtained using a DISA Type 
55D01 Constant Temperature Anemometer with a TSI-1210-T1.5 hot-wire sensor 
located on the duct center line, 120 mm (approximately 2 duct diameters) downstream 
of the bellmouth inlet. The anemometer was directly calibrated statically against 
the ASME orifice at the start of each test. To obtain good resolution of the hot 
wire signal, it was offset using a TSI Model 1057 Signal Conditioner (the filters 
on the Model 1057 were not used) which provided a 0-1 volt output over the flow 
range of interest. This signal was then low-pass filtered at 500 Hz and amplified 
by a factor of 10 using a Frequency Devices 901F 8-pole Butterworth filter. The 
conditioned anemometer output, which now ranged from 0-10 volts over the flow 
range of interest, was terminated in the A/D converter. At the start of each test 
the mass flow was measured using both the ASME orifice and the anemometer at 
ten steady operating points over the flow range of interest and the data were fit with 
a fourth order polynomial. The calibration curve thus obtained provided mass flow 
rate directly as a function of conditioned anemometer output voltage. 

Time-Averaged Inlet Temperature 

The time averaged temperature of the air entering the system was measured using a a 
type-K unshielded thermocouple inserted 65 mm downstream of the bellmouth inlet, 
radially located 6 mm from the duct wall. The thermocouple output was conditioned 
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using an Analog Devices 3B47 Isolated Linearized Thermocouple Input Module which 
provided a linear high level 100 degree-K/ volt output to be read by the A/D. 

Time Averaged Inlet Wall-Static Pressure 

Time averaged wall-static pressure was measured using a Druck-PDCR-820, 1 psi range, 
transducer located 65 mm downstream of the bellmouth inlet. The low level output of 
the Druck transducer was amplified using a Analog Devices 2B31K high performance 
strain gauge conditioner which also provided the required bridge excitation. The 
conditioned output provided a linear high level output of approximately 600 Pa/volt. 

2.3.2 Compressor Instrumentation 

Time Resolved Rotational Frequency 

A magnetic pickup in conjunction with a magnetized nut on the compressor shaft 
provided a sinusoidal signal at the rotor frequency. This sinusoidal signal was then 
fed to an Analog Devices Model 451 F/V (frequency to voltage) Converter. The F/V 
converter provided a linear high level voltage output of approximately 250 Hz/volt. 
This signal was then low pass filtered at 500 Hz with a Frequency Devices 901-F 8-pole 
Butterworth Filter. For the compressor frequency response measurements, greater 
resolution of the rotational frequency was required. For this purpose, the F/V output 
signal was split into two separate channels. One channel was used to determine the 
time-averaged rotational frequency the other was offset and then further amplified. 
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Time-Resolved Total Pressure 


Time-resolved total pressure at the compressor exit was measured using a 50 psi, 
Kulite XCE-093-50D pressure transducer connected to a modified United Sensor KB- 6 
Kiel Head Probe. The Kiel Probe was inserted 65 mm downstream of the compressor 
exit with the head of the probe at the duct center line. The modified probe arrange- 
ment is shown in Figure 2-7. As shown in Figure 2-7, the Kiel probe was modified by 
removing 3 inches (76mm) from the probe stem. The Kulite transducer was cemented 
with epoxy into a short section of 1/8” thin- walled stainless steel tubing with the sens- 
ing face of the transducer flush with the end of the stainless tube. The tube with 
the transducer is mated to the Kiel probe using a Scanivalve quick release connector. 
(This arrangement allowed for future re-use of the Kulite transducer.) The low level 
signal from the Kulite transducer was amplified using an Analog Devices 2B30K am- 
plifier which also provided excitation for the transducer. The resulting output was a 
linear high level signal of approximately 14000 Pa/volt. 

Time-averaged Static Temperature 

Time averaged compressor exit static temperature was measured with a type-K un- 
shielded thermocouple inserted 65 mm downstream of the compressor exit at a radial 
location 6 mm from the the duct wall. The thermocouple output was conditioned 
using an Analog Devices 3B47 Isolated Linearized Thermocouple Input Module. This 
module provided a linear high level 100 degree-K/volt output to be read by the A/D. 
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2.3.3 Control Valve Instrumentation 

Time-resolved Valve Position 

The angular position of the valve was measured using a BEI-MX-213-18-1024 opti- 
cal shaft encoder. The shaft encoder provided a resolution of 4096 quadrature counts 
per revolution. The encoder signal was decoded by the Galil Model DMC-400-10 

motion control board. 

Time-averaged Total to Static Pressure-Drop 

Time averaged total to static pressure drop across the control valve was measured 
using a Druck 20 psi PDCR-120/35-WL differential pressure transducer. The high 
pressure port of this transducer was connected to the Kiel probe described m Sec- 
tion 2.3.2 and the low pressure port is connected to a wall static tap in the plenum 
chamber into which the valve discharges. The time-resolved compressor exit total 
pressure and time-averaged control valve pressure drop were never both measured for 
the same test allowing the Kiel probe to be shared for the two purposes. The low 
level signal from the Druck transducer was amplified using an Analog Devices 2B30K 
amplifier which also provided excitation for the transducer. The resulting output was 
a linear high level signal of approximately 14000 Pa/volt. 
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2.3.4 Plenum Chamber Instrumentation 

Time-resolved Wall Static Pressure 

The time resolved wall static pressure in the plenum chamber was measured using a 
30 psi Druck PDCR-820-0479 transducer. The low level signal from the Druck trans- 
ducer was amplified using an Analog Devices 2B30K ‘amplifier which also provided 
excitation for the transducer. The resulting output was a linear high level signal of 
approximately 20000 Pa/ volt. 

2.3.5 Metering Duct Instrumentation 

Time-averaged Orifice Differential Pressure 

The time averaged differential pressure across the ASME orifice was measured using a 

2.5 psi Druck PDCR-811 transducer. The pressure tap for the transducer was located 
in the wall of the upstream flange. (Details are given in [39]). The low level signal 
from the Druck transducer was amplified using an Analog Devices 2B30K amplifier 
which also provided excitation for the transducer. The resulting output was a linear 
high level signal of approximately 900 Pa/volt. 

Time-averaged Orifice Temperature 

The temperature of the air entering the orifice was measured with a type-K un 
shielded thermocouple inserted upstream of the orifice. The thermocouple output 
was conditioned using an Analog Devices 3B47 Isolated Linearized Thermocouple In- 
put Module. This module provided a linear high level 100 degree- K/ volt output to 
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be read by the A/D. 


2.4 Digital Data Acquisition and Control System 

The digital data acquisition and control system has two basic functions: 1) Closing 
the compression system stabilization feedback loop with a digital control, and 2) 
storing data for further analysis. Due to limitations of the available hardware, it was 
necessary to integrate these two functionally independent tasks into a single system. 
The hardware system consists of: 1) Data Translations Model DT2801 A/D (Analog 
to Digital) converter, 2) 20MHz Intel 80386 microprocessor based Hewlett Packard 
Model HP-Vectra-RS/20 Microcomputer with a 120Mbyte Hard-disk, and, 3) Burr- 
Brown Model PCI-20093 W-l Analog Output Board (D/A converter). 

2.4.1 System Operation 

The cyclical operation of the system can be described functionally as follows: 

1. Measurements from the compression system are sampled by the A/D converter 
and put into a digital format 

2. The sampled signals are converted to engineering units. 

3. A new valve position command is computed based upon the most recent and 
(possibly previous) measurements. 

4. The valve position command is sent over the communication bus to the valve 
position controller (Galil digital motion control board). 
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5. Selected channels are output by the D/A converter to provide real time moni- 
toring on an oscilloscope. 

6. Data for selected channels is stored in the volatile memory (RAM) of the mi- 
coprocessor. 

At completion of the test, data stored in volatile memory is transferred to the hard 
disk for permanent storage. 

2.4.2 Loop Timing 

Details of the loop timing are given to help provide and understanding of the data re- 
duction procedures employed. The overall data acquisition and control cycle actually 
involves a number of independent, concurrent, processes. These processes are: l)The 
software loop, 2) the A/D conversion process and, 3) the actuator position control 

process. 

The actuator position control process is executed by the Galil digital motion 
control board as described in Section 2.2. The A/D conversion process involves the 
Data Translation A/D board and the DMA (Direct Memory Access) controller of 
the Intel 80386 Microprocessor. The A/D board executes an endless loop, in which 
it successively samples the first, second, third, up to the n th measurement channel, 
and then repeats. The board contains an internal clock and at each “tick” of this 
clock the next channel is sampled. After each channel is sampled, the DMA controller 
transfers the digitized value directly into the memory of the 80386 Microprocessor. 
There is a single, fixed memory location in the 80386 assigned to each A/D channel 
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and this memory location is overwritten each time the A/D obtains a new digitized 
sample for that channel. Thus, if there are a total of n channels of data which are 
measured, and the clock ticks at a uniform rate of / ticks/sec, each channel will be 
sampled every nff seconds and, correspondingly, each assigned memory location in 
the 80386 will also be overwritten every n/f seconds. The amount of time, n/f will 
be referred to as the A/D cycle time or sampling period, T a . With this system, if n 
channels of data are read by the A/D in each scan, then the first and last channel in 

a data scan are separated by a time interval of T, — 1//. 

For control and data analysis purposes it is preferable to sample all the channels 
simultaneously. The process described above was adopted by necessity to accommo- 
date various limitations imposed by the available hardware. 

The software loop is a program which is executed by the 80386 microprocessor 
itself. It is responsible for processing the data obtained by the A/D, converting it to 
engineering units, implementing the digital control algorithm, storing data to an in 
memory buffer, sending commands to the actuator position control, and communi- 
cating with the user via keyboard input. This program operates as an endless loop 
(until it is terminated by keyboard input from the user.) The cyclic operation of the 
software loop is as follows: 

1. The program enters a “tight loop”, checking if all of the A/D channels have 
been updated (overwritten) that is, it waits until the completion of the next 

A/D cycle. 

2. As soon as it has been determined that an A/D cycle is completed, all the 
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updated A/D channels are “read” i.e. they are copied to another memory 
location where they will not be overwritten. This set of values is called the 

current A/D scan. 

3. The valve position command is computed. 

4. The program enters a tight loop, polling the Galil DMC (digital motion control) 
board to determine if it is ready to communicate. The DMC board communi- 
cates at uniform time intervals of 0.5msec, corresponding to its internal 2KHz 
cycle time. This clock is not synchronized with the clocks on either the 80386 or 
the Data Translations A/D board. At the end of each DMC cycle, it sends the 
current motor shaft position and receives the current position command. Since 
the DMC uses the same register in the 80386’s I/O address space for sending 
and receiving data, the shaft position must be read before the command can be 

sent. 

5. Once the DMC board is ready to communicate, the current shaft position is read 
and the new position command is sent. Both the incoming and outgoing data 
must be sent one byte at a time with “hand-shakes” occurring after each byte. 
Thus, the 80386 is unavailable for other tasks until all the data is transferred. 

6. Once the data transfer is completed, the current A/D scan is converted to 
engineering units, and the digital control law calculation is performed. 

For control and data analysis purposes it is necessary to have the cycle time for 
the software loop the same as that of the A/D. Since the steps described above are 
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executed serially, it can be seen that the total amount of time required for the software 
loop to execute can be considered to be the sum of three time periods: 1) the waiting 
time for the A/D cycle to complete, 2) the waiting time for the communication with 
the DMC to complete and, 3) the computational time for implementing the digital 
control and other auxiliary calculations, such as conversion to engineering units. To 
insure proper timing of the software loop, the A/D cycle time thus had to be longer 
than the sum of the DMC waiting time and the computational time. This put an 
upper bound on the A/D cycle time of 1msec providing a sampling rate of lKHz per 

channel. 
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Figure 2-1: Overview of experimental facility 
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Figure 2-5: Control valve used for experiments 
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wall static pressure 



Figure 2-6: Sensor locations for experimental studies 
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Figure 2-7: Modified Kiel Probe used to measure compressor exit total pressure 
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Chapter 3 


System Modeling 

In this chapter the system dynamics of a compression system equipped with a number 
of actuators and sensors are modeled in a form suitable for control design and analysis. 


3.1 Description of System 

The compression system to be studied is shown schematically in Figure 3-1. The 
system consists of an inlet duct, compressor, plenum chamber and throttle valve, and 
a number of actuators and sensors. This configuration is typical of that found m a 
gas turbine engine where the plenum chamber would represent the combustor and the 
throttle represents the turbine nozzles. All of the actuators and sensors shown here 
for implementing a stabilizing control would not, however, be found on a production 

engine. 

Ambient, low pressure gas is drawn into the inlet duct where it flows to the com- 
pressor. The compressor then raises the pressure of the gas, discharges the pressurized 
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gas into the plenum chamber, where combustion occurs, raising the gas temperature. 
Hot, pressurized gas exits the plenum chamber, through the turbine nozzles, and is 
exhausted through the turbine to the ambient surroundings. 

The system to be analyzed is equipped with five actuators which provide alter- 
native means to influence the behavior of the system and four sensors to provide 
feedback signals. This set of actuators and sensors were selected as being representa- 
tive of a diverse number of potentially attractive options. In the inlet duct upstream 
of the compressor is an injector which can be used to inject additional gas into the 
compressor inlet. The flow through the injector is controlled by modulating the in- 
jector supply pressure. A close-coupled valve is located immediately downstream of 
the compressor. The pressure drop across the valve is controlled by modulating the 
open area of the valve. A plenum bleed valve located in the plenum wall discharges 
high pressure gas from the plenum to ambient conditions. Flow through the valve is 
controlled by modulating the open area of the valve. A piston, which will be referred 
to as the movable plenum wall, is used to vary the volume of the plenum. The volume 
of the plenum is controlled by inserting or retracting this piston. The fifth means of 
influencing the system is to vary the rate of heat addition to the plenum. In a gas 
turbine this could be accomplished by varying the fuel flow. 

Four sensors are available to measure the behavior of the system. Total and static 
pressure sensors are located at the compressor face. A mass flux sensor (such as 
a hot wire or hot film) also located in the compressor inlet, is used to obtain the 
instantaneous mass flow rate through the compressor. Since one-dimensional axial 
flow will be assumed, the radial locations of the transducers in the duct need not be 
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specified. A pressure sensor in the plenum wall provides measurements of the static 
pressure in the plenum. Since uniform pressure in the plenum is assumed, it is not 
necessary to specify the location of the sensor. All of the sensors are assumed to 
provide perfect measurements of the sensed quantities, free of noise and any other 
inaccuracy. This considerable idealization can be relaxed later by further augmenting 

the model. 

3.2 Component Models 

The dynamics of the compression system described above will be approximated with 
the lumped parameter model depicted in Figure 3-2. The model is constructed using 
a number of components which will be first modeled individually and then inter- 
connected. All locations (flow stations ) referred to in the following derivation are 
identified in figure 3-2. The presentation here will be brief, highlighting the major 
assumptions and the basic physics of the model. 

3.2.1 Ducts 

The inlet duct is broken into two separate ducts; one upstream and one downstream 
of the injector. These will be referred to as the upstream and downstream ducts 

respectively. 
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Assumptions 


Flow in the ducts is assumed to be incompressible, inviscid, and one-dimensional 1 . 
The duct area is not assumed to be constant but must vary slowly enough to justify 
the one-dimensional flow assumption. 


Mathematical Model 

With these assumptions, the flow in the duct is completely described by specifying 
the mass flow, m through the duct and this will be used as a state variable. The time 
rate of change of the mass flow through the upstream duct, m«, is found from the 
axial momentum balance: 


dm u /dt = (Poi — Ptb)? 

*-'uj 

where p 0l and p 02 are total pressures at station- 1 and station-2 respectively; A c is the 
annular inlet area of the compressor, and C u is the equivalent length of the upstream 
duct computed using A c as the reference area 2 . Similarly for the downstream duct: 


drrid/dt = 


1 

ZJa c 


(Po, -P0,)> 


(3.2) 


i B y one-dimensional it is meant that at all points in the duct the fluid moves only in the axial 
direction, and at any given axial location the axial velocity is constant across the duct cr °“ sectl °£ 
2 The equivalent duct length for a duct of non-constant cross sectional area is defined to be he 
length of duct with a given constant area which would experience the same acce eration as the 
nomconstant area ductwhen the same total pressure difference was applied L Usrng ^ 
annular inlet area as the reference area, the equivalent length, £, of a duct of length L, g 
by . c = A c 1 /A{x) dx , where A(x) prescribes the variation of duct area as a function of axia 

distance along the duct center line. 
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The equivalent length of the downstream duct, £ d , includes the flow path through 
the compressor. 

3.2.2 Plenum 

Assumptions 

1 

All fluid properties in the plenum are assumed to be spatially uniform, that is, the 
plenum is well mixed. Fluid velocities in the plenum are assumed to be negligible. 
The effect of combustion is, idealized as a specified rate of heat addition to the gas 
contained in the plenum. Ideal gas behavior 3 , with constant specific heat is assumed. 

Mathematical Model 

With the above assumptions, the plenum behavior is fully described by specifying 
the mass, m p , internal energy, U p , and volume, V p of the gas in stored in the plenum, 
which may vary with time due to the movable wall. These quantities will thus be 

used as the state variables for the plenum. 

The time rate of change of the mass of gas, m p , in the plenum is found from the 

mass balance: 

dm p /dt = (m 6 — (m 7 + m 8 )), (3-3) 

where m 6 is the incoming mass flow, and m 7 and m 8 are the mass flows exiting 

through the throttle and bleed valve respectively. 

The time rate of change of internal energy, U p , of the gas in the plenum is found 

3 The pressure, temperature and density are related by p — pRT 
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from the energy balance: 


dU p /dt = riieCpTos - rh 7 C p T 07 - Th 8 C p T 0s + Q + p P A w u w , ( • ) 

where, C p is the specific heat at constant pressure, T 0o , T 0l , and T 0g are the absolute 
total temperatures of the entering and exiting streams, Q is the heat addition rate, Pp 
is the absolute static pressure in the plenum, A w is the projected area of the movable 

wall and u w is the velocity of the wall. 

The time rate of change of the volume of the plenum is given by: 

dVpjdt = A w u w ( 3 - 5 ) 


3.2.3 Compressor 

Assumptions 

The compressor is modeled as an actuator disk, that is, a region of infinitesimal axial 
extent producing a discontinuity in total pressure and temperature. Any fluid inertia 
associated with the flow passages of the actual compressor is included in the equivalent 
length, C d , of the duct attached to its upstream flange. The actuator disk assumption 
implies that no mass can accumulate in the compressor. Unless otherwise noted, the 
instantaneous total pressure-rise and total temperature ratio of the compressor are 
assumed to be the same as the steady state values for the same inlet conditions. This 
will be referred to as the quasi-steady assumption. 
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Mathematical Model 


The total to total pressure rise, {po s -po 4 ), the compressor is given by. 

p 0i _ po< = 1/2 p a u T 2 'H c (m 4 /(p a u T A c ), u T /a), ( 3 - 6 ) 

where p a is the ambient density, u r is the tip speed of the compressor, m 4 is the mass 
flow through the compressor, a is the speed of sound at ambient conditions, and * c 
is the non-dimensional compressor pressure rise characteristic. is an empirically 
determined function which must be measured for the compressor of interest. Similarly, 
the compressor temperature ratio, 7o 5 /7o 4 , is given by. 

T 0 JT 0t = Q c (m 4 l(paUTAc),u T /a) ( 3 - 7 ) 

where, 0 C is the total temperature ratio characteristic of the compressor. This is also 
an empirically determined function. 

In some cases, we will not make the quasi-steady assumption. Instead, we will 
assume that the instantaneous, non-dimensional pressure rise of the compressor, V>(0, 
is given by: j>{t) = where ^ci(rh 4 /(p a u T A c ),u T /a) is the theoretical, ideal 

pressure rise of the compressor, and *,(*) is the instantaneous non-dimensional total 
pressure loss due to viscous effects and other non-idealities. The losses will be assumed 

to evolve in time according to: 

djJdt = - (¥ ci (rh 4 /(p a u T A c ),u T /a ) - * c (rh 4 /(p a u T A c ),u T /a ) - </>,(*)) , ( 3 . 8 ) 

ti 
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where ti is an empirically determined time constant. This will be referred to as the 
unsteady loss model. 

3.2.4 Throttle 

Assumptions 

The throttle flow is assumed to be quasi-steady, one dimensional, compressible, and 
isentropic. 

Mathematical Model 

The throttle is modeled as a quasi-steady one-dimensional isentropic converging noz- 
zle. The mass flow through the throttle, m 7 is given by: 

m 7 = p a UTAc^t{Pp/PatTp/T a ,At/ A c ), (3-9) 

where, p p and T p are, respectively, the absolute temperature and pressure in the 
plenum, p a and T a are the ambient absolute pressure and temperature and A t is 
the exit area of the throttle. The function, is the non-dimensional throttle flow 
characteristic which is obtained assuming one- dimensional, steady, compressible flow. 
The nozzle model is applicable for both choked and non-choked conditions. 
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3.2.5 Plenum Bleed Valve 


Assumptions 


The plenum bleed valve was also modeled as a quasi-steady, isentropic, converging 


nozzle in exactly the same manner as described for the throttle. 


Mathematical Model 

The mass flow rate through the bleed valve is given by: 

rh 8 = p a UTA c $b(Pp/Pa,T p /T a ,A b /A c ), (3.10) 

where is the bleed valve area, and the function $ b is the non-dimensional bleed 
valve flow characteristic. (It is actually the same function as ^t-) 

3.2.6 Close-Coupled Valve 

Assumptions 

The close-coupled valve, as its name implies, is located immediately at the exit of the 
compressor. It is assumed that mass flow rates entering and leaving the close-coupled 
valve are at all times equal to each other and to the flow through the compressor. 
That is, there is no accumulation of mass in the close-coupled valve or between it and 
the compressor. Pressure loss across the close-coupled valve is assumed to be small 
enough for the flow through the valve to be accurately modeled as incompressible. 
The valve is assumed to discharges into the plenum with negligible static pressure 
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recovery. 


Mathematical Model 

The total pressure drop across the valve, (po 5 -Po 6 ) 1S S iven ^ 

(Po 5 — Poe) = l/lp^T 2 y v {m h l{p a u T A c ),OL c ), (3.11) 

where, m 5 is the flow through the valve, a c is the valve fraction open (the area of 
the valve divided by the area of the valve when wide open), and is the non- 
dimensional pressure drop characteristic of the close-coupled valve. This function is 
either estimated analytically or measured empirically. 

3.2.7 Injector 

The section of duct between station-2 and station-3 along with the injection tube will 
be referred to as the injector. 

Assumptions 

It is assumed that the injector behaves quasi- steadily, that is, that the steady state 
characteristics continue to hold even when the upstream and downstream mass flows 
or supply pressure are unsteady. It is assumed that the mass of fluid contained in 
the injector duct remains constant, and that the axial velocity is spatially uniform 
at station-3, (i.e., the flow is fully mixed out at the exit of the injector). The flow 
between station-2 and station-3 is not necessarily one-dimensional. All of the fluid 
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inertia of the injector is included in C, u . 


Mathematical Model 

The injector must be characterized as a two port device. Specifically, the difference 
jjj pressure between station-2 and the injector supply, (P02 POj)> given by 

P(h - P0t = l/2p a u r 2 ^i u (Th 2 /(paUTA c ),rh 3 /(paU T A c )), (3.12) 

where $, u is defined as the non-dimensional injector upstream pressure characteristic, 
and m 2 and m 3 are the upstream and downstream mass flows, respectively. Similarly, 
the difference in pressure between station-3 and the injector supply, (po 3 — Po,)» i s 
given by: 

p 0i - p 0i = 1/2 p a u T 2 i d (m 2 /(p a u T A c );m 3 /{p a u T A c )), (3.13) 

where 'P, is defined as the non-dimensional injector downstream pressure character- 
istic. 

Analytical expressions which are used for the steady characteristics and 
are derived in appendix A. 

3.3 System Model 

The component models described in the previous section must be interconnected to 
form an overall system model. Applying continuity, matching pressures at component 
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interfaces, and using the spatially uniform plenum assumption gives: 


dm u /dt 
dthdl dt 
dm p /dt 
dU P ldt 
dVp/dt 


1/2 p a UT 2 ( ( Pa -POi) _v.\ 

Cut Ac \ l/2p a uj 2 ") 

\J2pmL (V + (p* Z lal + *'-y v - & ' ~ 

£ d /A c V ’ d l/2/>oUr 2 l/2p a ur / 

p a ujA c ( rh d /(p a u T A c ) - - $!>) , 

(a„c p t„0 c - miWA + + pp- 4 -"”) 


(3.14) 


For brevity, the arguments of the various characteristic functions have been sup- 
pressed. The system model, Equations 3.14, is put into non-dimensional form by 
normalizing the various dimensional quantities by the reference values given in Ta- 
ble 3.1 to form the non-dimensional quantities listed in Table 3.2. For example, non 
dimensional time, r, is defined as the ratio of the physical time, t, to the reference 
time, (•/cJ r )/(a p ^A c ), from Table 3.1 to give: 


t 
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quantity 

reference 

value 

time 

a/C.V'p 


pressure(differential) 

l/2p a u T 2 

pressure( absolute) 

Pa 

mass flow rate 

palirAc 

temperature 

T 

± a !i 

heat addition 

p a UTA c CpT a 

area 

h A c 

plenum mass 

C3 

plenum energy 

PaV p C v T a 

plenum volume 

Vp 

tip speed 

a 

length 

Ci 

11 

hi 

e 

+ 

a. 

wall displacement rate 

ujA c 


Table 3.1: Reference quantities for non-dimensionalizing the system model 


It is also useful to utilize the B parameter 4 which is defined by: 


B = 


UT Vp 
2d p \| ( C c A c ) 


(3.15) 


The resulting non-dimensional equations are then given by: 


dfa/dr 

d<f> d /dT 

dfrip/dr 


-B 

CjCc 

B 

C d /C c 
Mt 2 
2Bt % 


(fa + ^i v (fa,<f>d)) 

4>d) + fa + M t) - *v(fa, a c) - fa) 


(3.16) 


4 One could also define a stability parameter, (Bip pea k)/<t>peak, where <j> pta k is the 
at which the peak compressor pressure rise occurs and *..* is the peak pressure rise. Fo glasses* 
compressors where the slope of the compressor characteristic roughly scales with . this 

stability parameter provides a useful means to compare the stability characteristics of compression 
systems containing compressors whose pressure-rise differs greatly. 
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non-dimensional 

symbol 

quantity 


time 

r 

pressure (differential) 

xP 

pressure (absolute) 

7 r 

mass flow rate 


temperature (absolute) 

9 \ 

heat addition 

Q 

area 

A 

1 plenum mass 

rh p 

plenum energy 

U p 

plenum volume 

% 

tip speed 

Mt 

wall displacement rate 

£ 

length 

L 

Greitzer B parameter 

B 


Table 3.2: Non-dimensional symbols for the system model 


dU p /dr 

dV v jdr 


MzSL (<f> u Q c (<j>d, Mr) — ($<(*> 7 0p) + ^«>( 7r p> ^p))^p + Q + x p() 

2B9 p ' 

BM w Mt e 


To obtain closure, some auxiliary relations are required to allow all the quantities 
appearing on the right hand side of equation 3.16 to be determined in terms of the 
state and input variables. Using the definitions of the non-dimensional variables and 
routine algebra gives the required relations as: 


7T p 

= %!% 

(3.17) 

V’p 

= 210,1%- 1)/(7-Vt 5 ) 

(3.18) 

dp 

= l Op/ rhp 

(3.19) 
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3.3.1 Model Reduction 


The model order is reduced if either the injector or movable wall is not used. Without 
the injector, the flow through the upstream and downstream ducts is the same and 
the model reduces to: 


d<f>/ dr 
drh v jdT 
dUp/dr 
dVp/dr 


( 3 - 20 ) 

— ^p) + ^p( 7r p> + 0 + **) 

2nv p 

BM w Mt c 

W v 4 ’ 


where is now the common mass flow through both the up and downstream duct. 
With a fixed plenum volume, the system equations reduce to. 


dKlir = -T^r + < 3 ' 21) 

L'u! *^c 

dfa/dr = + ^ + 

L'diL'c 

drh p /dT = —-j^r{<f>d — {$t('XpiOp) + ®b('’ r p'>0p))) 

2t D U p 

dUp/dr = (<l> u Qc(<t>d, Mt) - ($<(*>, <? P ) + $b(*p,9p)Wp + Q) • 

( 3 . 22 ) 


With no injection and a fixed plenum volume, the system equations become: 

d(t>/dr = ( 3 - 23 ) 
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dm p /dr = + 

Zntlp 

dUp/dr = (^Qc(^i^t) _ (^t( 7r p>^p) + $&( 7r p>0p))^p + Q) 

2 13 u p 

(3.24) 

Removing any of the other actuators does not require a change in the number of states. 
If the close-coupled valve is not used, the valve pressure drop characteristic, is held 
fixed at zero. When no bleed valve is present, the plenum bleed characteristic, $ 6 , 
is set to zero. With no unsteady heat addition Q is held constant, although not 

necessarily zero. 


3.4 Simplified System Model 

The system model can be further simplified by assuming that the thermodynamic 
process occurring in the plenum is isentropic, i.e., p/p 7 = constant, where 7 is the 
ratio of specific heats. With this assumption, and considering only the close-coupled 
valve and plenum, the following second order (two state) model is obtained: 

dcfr/dr = -^p + Vv-Vc) ( 3 - 25 ) 

diftp/dr = Jjid’d — ($< + $<>))• (3.26) 

This model, given by Equations 3.25 and 3.26, will be referred to as the basic com- 
pression system model. 
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3.5 Linearized System Model 


The preceding nonlinear models are valid for perturbations of arbitrary amplitude. 
For many purposes it is useful to have simpler models which are valid for small per- 
turbations about an equilibrium operating point. Throughout, when such models are 
required they are obtained by the usual Taylor series expansion about an equilibrium 
point (see for example [24]). We will subsequently have need to refer to the linearized 
basic compression system model which is given by: 

d4>cldr 

dx/jp/dr 


m r B 


1 


-B 


i 

Bm T 


to 

Wp 


+ 


Bv c 0 
0 










. (3.27) 


where, m c 





and, Vb = 

C 



are evaluated at the 


equilibrium operating point of interest. 
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Q Heat Addition Plenum 



Figure 3-1: Compression system equipped with actuators and sensors 


Plenum 



Figure 3-2: Flow stations for for modeling compression system 
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Chapter 4 


Validation of Compressor Model 

4.1 Introduction 

Central to the modeling of the compression system is the treatment of the compressor 
itself. As discussed in previous chapters, most compression system models assume 
that the compressor can be represented by a possibly nonlinear fluid-resistor (quasi- 
steady actuator disk) in series with some amount of fluid-inertia. In some cases, this 
model is elaborated slightly by using a first order lag to relate the compressor instan- 
taneous pressure rise and the quasi-steady value at the same flow rate. The use of a 
low order, lumped element, (finite number of ordinary differential equations) repre- 
sentation is common to both these approaches, however, representing a considerable 
simplification of the detailed behavior of the flow field existing in an actual compres- 
sor. Although much indirect support for this approach is provided by the success of 
the overall models which incorporate these assumptions, little direct validation of this 
modeling approach has been done previously. To provide such a direct validation, a set 
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of experiments have been performed in which the transfer function between sinusoidal 
pressure perturbations across the compressor (input) and the mass flow perturbations 
through the compressor (output) was measured. The measured frequency response 

could be directly compared with that predicted by a lumped model. 

In the following sections the theoretical frequency response of the compressor is 
first obtained for the compressor with and without a first order lag. Next, the exper- 
imental methodology used to obtain experimental measurements of the compressor 
transfer function is given. Finally, the experimental measurements and theoretical 
response are compared and the results interpreted. 


4.2 Theoretical Compressor Frequency Response 

4.2.1 Background- Admittance, Impedance, and Phasors 

Consider the compressor shown schematically in Figure 4-1. Let us assume that the 
the non-dimensional mass flow rate at inlet and exit, & and <f> 2 are at all time equal 
and refer to this common flow as <f>. Suppose that the difference (exit minus inlet) 
in non-dimensional total pressure, *, across the compressor is caused to fluctuate 
sinusoidally about some mean value rp, with amplitude, ip m , so: 

xp(t) = ? + 


where: 

ip(t) = ipm cos(u )t + #l) 
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For a linear time invariant system, which is assumed here, the resulting perturbations 
in mass flow rate, J(t), are also sinusoidal at the same frequency with amplitude, <f> m 

and phase angle 0 2 : 

4>(t) = <f>m COS (wf + 0 2 ) 

It is convenient (and quite standard, see for example Chua [7] ) to represent these real 
valued, time domain, sinusoidal signals as phasors. That is, we associate the unique 
complex number (phasor), tf> = with the pressure signal, </>(*) and the unique 

complex number (f> = with the mass flow rate signal (Note that bold 

face will be used to distinguish phasor quantities.) The corresponding time domain 
signals can be uniquely recovered using: t />(*) = Tle{^} and, = lU&e?'*}, 
where, 7 le{}, indicates the real part of a complex number. Using the phasor notation, 
we define (in analogy to electrical system terminology) the complex valued, frequency 
dependent admittance, Y(ju), of the compressor as the ratio of the flow phasor, 0 , 
and the pressure phasor, rj> 

Y (ju>) = 

The sign convention adopted here means that the direction of positive power flow is 
into the element. Thus, the magnitude of Y{ju) is given by: \Y(ju)\ = KH’m which 
is the amplitude ratio of the output flow perturbation to input pressure perturbation. 
Similarly, the angle of Y(ju), which is given by: l(Y(ju)) = 0 2 - 9 U represents the 
phase shift between the mass flow perturbations and pressure perturbations. It will 
also be useful to define the complex valued, frequency dependent impedance, Z(ju), 
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of the compressor as the inverse of the admittance as given by: 

Z{ju) = 1 /Y(ju). ( 4 ‘ 2 ) 

The frequency dependent behavior of either Y(ju) or Z(ju) will be referred to as the 
frequency response of the compressor. 

Analytical expressions for the compressor frequency response will now be obtained 
for both the quasi-steady and the unsteady loss models of the compressor as described 

in Section 3.2.3. 

4.2.2 Quasi-steady Model 

From Section 3.2.3 we have for the linearized quasi-steady compressor: 


d(t>/dt = utI(2C c )(^4> ~ 0)- 


(4.3) 


For now, we maintain the use of the physical time, t. Assuming a sinusoidal input 
pressure perturbation and flow response, substituting into Equation 4.3 and perform- 
ing routine manipulations yields: 


Y(j“) = 



( 4 . 4 ) 
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where the compressor resistance, R e is defined as R c = The impedance can 

then be obtained by applying Equation 4.2 to Equation 4.4, giving. 

Z(ju>) = R c + j2C c u/uT ( 4 - 5 ) 

Thus, for the quasi-steady model, the real part of the impedance is (within an alge- 
braic sign) simply the compressor slope, and is independent of the forcing frequency. 
As a function of frequency, the imaginary part of the impedance would plot as a 
straight line with a slope in direct proportion to the equivalent duct length (fluid 
inertia). The impedance representation thus isolates the two important physical ef- 
fects, flow dependence of the compressor pressure rise and fluid inertia, into real and 
imaginary parts, respectively, of the compressor frequency response. This provides 
a means for directly checking the validity of the quasi-steady assumption. Non- 
quasisteady behavior effecting the compressor pressure rise will manifest itself as a 
frequency dependence on Ke{Z(ju)}, and and dlm{Z(ju)} /du. 

4.2.3 Unsteady Loss Model 

Using model described in Section 3.2.3 the linearized model for the compressor with 
unsteady losses can be obtained as: 

dt/dt = u T /{2C c )(^ i 4> -'Fl-'P) ( 4 - 6 ) 

dfjdt = -$*-«>. < 4 - 7) 
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where r L is defined as the ratio of the loss time lag to the compressor flow through 
time. We will only be concerned here with radial bladed machines for which the 
ideal pressure rise is a constant, independent of flow, so that, = 0. Once 

again defining R c = the admittance and impedance can be obtained as in the 

preceding section. The admittance is then given by. 


Y{ju) = 1 /Rc— 


f 1 + 


(-2C c t l Vc2 + l) + 

t WZXX7 ^ > RcU T j 


(4.8) 


R c ut 2 A { 


This reduces to the quasi-steady model, given by Equation 4.4 as the flow-through 
time becomes small with respect to the drive frequency, i.e. t l u -+ 0. The impedance 
is found from Z(ju>) = 1 /Y(ju?) and is given by. 


R c + j((2CJur - 


ZO") 




(f>UT A 


l+( 


<$>utA c 


(4.9) 


As previously noted, with the unsteady behavior modeled, the real part of the impedance 
becomes frequency dependent. In the limit, as the frequency approaches zero, the 
real part of the impedance approaches the slope of the compressor pressure rise char- 
acteristic but at high frequency it approaches the slope of the ideal pressure rise 
characteristic, which is zero in the present case. The derivative with respect to fre- 
quency (slope) of the imaginary part of the impedance, is frequency dependent. For 
the unsteady model, even at low frequencies, the slope of the imaginary part of the 
impedance is dependent on the time mean flow coefficient as well as geometry. 
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4.2.4 Frequency Normalization 


The expressions for the theoretical frequency response of the compressor derived in 
the preceding sections have been in terms of dimensional angular frequency (e.g. 
radians/sec). It is useful in interpreting data to normalize the driving frequency to 
some relevant physical time scale. For this problem there are two different important 
time scales. The first, is the exponential growth (decay) time constant inherent in 
the linearized, quasi-steady model given by: 

t, = \2 CJ(R c u t )\, ( 4 - 10 ) 

which is seen to be a function of the compressor slope and fluid inertia (equivalent 
length) of the the compressor duct. This is analogous to the L/R time constant of 
an electrical circuit consisting of a an inductor with inductance (L) in series with a 
resistor with resistance (R). Based on this time constant, we define w fl , the quasi- 
steady normalized frequency by: 


«,=«/(! /*,)• (4 ' U) 

The other relevant time scale is the flow through time of the compressor. The flow 
through time is defined here as the compressor volume, V„ divided by the volume 
flow rate of fluid, which can be expressed in terms of the flow coefficient and tip speed 

as: 

K = VJ{<j> c u T A c ), ( 4 - 12 ) 
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from which the reduced frequency u u is defined as 


Uu =u/(l/tu). ( 4 ‘ 13 ) 

The flow through time, as defined by Equation 4.12, is the amount of time it would 
take a fluid particle in a constant density, axisymmetric flow (plug flow) to pass 
from the inlet to the exit of a conduit whose volume matched that of the compres- 
sor. Note that the cross sectional area of the conduit is irrelevant for this calcula 
tion. It is remarked that the flow through time of a centrifugal compressor with a 
volute is a somewhat poorly defined quantity in that the flow through time is de- 
pendent on the particular path taken through the compressor, and the density is not 
in fact perfectly constant. Thus, the definition of flow through time adopted here 
is by no means unique, but this representation should provide the correct order of 
magnitude and capture the major dependence on mean flow and geometric scaling. 
The volume of the Holset HID compressor used for the experiments was estimated 
at V c = 2.36 x 10- 4 m 3 , based upon a rough approximation of the annular volume 
bounded by the centrifugal impeller hub and the casing. 

4.3 Experimental Methodology 

4.3.1 Summary 

Experimental frequency response measurements were obtained for the Holset HID 
centrifugal compressor described in Chapter 2 using the short inlet configuration. 
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The frequency response was measured at four different time mean operating points 
along a constant corrected speed line with M r = 0.76, as shown m Figure 4-2. The 
flow coefficients at these time mean operating points, which range from 0.1 to 0.2, were 
chosen to include both negatively and positively sloped portions of the compressor 
characteristic. For the latter, the downstream valve acted to stabilize the system, 
thus allowing the compressor to operate without surging even on the positively sloped 


portions of its characteristic. 

To obtain frequency response data, the compressor was forced sinusoidally by 
modulating the area of a valve located immediately downstream of the compressor 
exit. The system was driven at 25, logarithmically spaced, values of frequency, ranging 
between 1 Hz and 100 Hz. This provided data on the compressor behavior from the 
low to high frequency asymptotic limits based upon the quasisteady time scale given 
by Equation 4.10. For each drive frequency, the system was operated for several 
minutes prior to taking data which allowed any transient behavior associated with 

the start up of the forcing to die away. 

Time resolved measurements were made of the total pressure difference imposed 


across the compressor (inlet to exit) and the resulting flow rate fluctuations through 
the duct. The compressor exit minus atmospheric inlet pressure difference was mea- 
sured using the high response total pressure probe described in Chapter 2. Flow rate 
measurements were derived from a single hot wire anemometer located at the duct 


center line, several duct diameters upstream of the compressor. 


The locations of the 


instruments are shown in Figure 4-3. 

The signals were Fourier analyzed to obtain their amplitude and phase at the drive 
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frequency, from which the desired frequency response was then directly calculated. 

4.3.2 Detailed Test Procedure 

The basic configuration used for the frequency response measurements is shown 
schematically in Figure 4-3 with the various instruments and transducers as described 
in in Chapter 2. The frequency response data for the compressor was obtained at 
four different time mean operating points along a line of constant compressor speed 
(M T = 0.76) as shown in Figure 4-2. The time mean flow was set to the desired value 
by adjusting the nominal area of the control valve. (The throttle at the plenum exit 
was wide open for these tests) A series of test runs were then performed in which the 
system was driven at a particular fixed forcing frequency by the control valve. 

A collection of such runs at a fixed time mean operating point will be referred to 
as a test. Each test consisted of 25 runs in which the frequency was varied between 
1 Hz and 100 Hz which gave data over a wide range of reduced frequencies for the 
quasi-steady model; from the low frequency to high frequency asymptotic limits. A 
logarithmic frequency spacing was used. The drive frequency was set using the analog 
Wavetek signal generator which was then sampled by the computer and used as an 
input to the digital controller which, in turn, drove the servo-motor position control 

loop. 

The length of time and sampling rate of the runs was varied according to the drive 
frequency (at least 10 samples per period) to provide sufficient spectral range and 
total number of cycles (at least 50 cycles) to obtain sufficient spectral resolution. The 
cutoff frequency of the low pass anti-aliasing filters was set at one half the sampling 
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frequency, to avoid aliasing. 

The amplitude of the drive signal was adjusted at each frequency to the minimal 
value required to obtain a clearly distinguishable flow perturbation as determined 
by direct observation with a digital storage oscilloscope. This method was adopted 
to minimize the effects of non-linearity, which increase with signal amplitude, while 
maintaining a large signal to noise ratio. Post-test, spectral analysis showed that 
the resulting sinusoidal amplitude of the flow coefficient perturbation at the drive 
frequency varied between 0.005 and 0.038 over all the tests conducted. Amplitude of 
the forcing pressure coefficient perturbations at the drive frequency ranged from 0.01 

to 0.1. 

To account for effects of rotor speed fluctuation, it was also necessary to ex- 
perimentally obtain estimates of the rate of change of the compressor time mean 
pressure-rise characteristic with speed i.e., This data was obtained at the four 

time mean flow coefficients at which the frequency response was measured. At each 
flow coefficient, the compressor steady state pressure rise was determined at ten dif- 
ferent rotational frequencies equally spaced to cover a range of ±15 Hz nominally 
centered at the time mean rotational frequency as used for the frequency response 
tests. (This range of flows was chosen to reflect the observed maximum fluctuations 
in rotational frequency) At each rotational frequency the throttle was adjusted to 
maintain the desired time mean flow coefficient. Over the small range of speed vari- 
ations, the resulting data showed an essentially linear variation of pressure-rise with 
rotor speed, as would be expected. The slope of a straight line fit to this data was 
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used accordingly to provide an estimate of 

4.4 Signal Processing and Data Analysis 

The overall goal of the signal processing and data analysis was to obtain an accurate 
measurement of the frequency response of the compressor from the raw, time series 
(digitally sampled) data. Specifically, for each drive frequency at which the compres- 
sor was tested it was desired to obtain both the magnitude ratio and phase difference 
between the input pressure signal and output fluctuation in the compressor mass flow. 
The methods applied for this purpose will now be described. 

Drive Frequency Determination 

The first step in the data analysis was to determine the actual drive frequency, /<*• 
Because the drive frequency signal was generated externally by the Wavetek signal 
generator and its frequency was not known exactly. The drive frequency was deter- 
mined by using the DFT (Discrete Fourier Transform) of the finite length sequence, 
x [n], obtained by sampling the drive signal over the duration of a run, where the DFT 
is defined as in Oppenheim [37], for a length N sequence by: 

X[k\ = I] x[n]e~^ N)kn . 

n=0 

For computational efficiency, the DFT was actually computed using a FFT (Fast 
Fourier Transform) algorithm rather than directly from the definition above. The 
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drive frequency was then given by the spectral component of the DFT with the 
greatest magnitude. As discussed in Oppenheim [37, Example 11.5], the frequency 
corresponding to an isolated spectral peak can be more accurately determined by 
zero padding the sequence (i.e. extending the length by appending a sequence of 
zeros) prior to computing the DFT. An iterative procedure was implemented for this 
purpose which would double the sequence length by zero padding until successive 
estimates of the peak magnitude agreed within a prescribed tolerance. (A value of 
1% was used for the analysis reported here.) The frequency corresponding to this 
peak was then taken as the drive frequency for that run. Based upon the minimal 
spectral resolution obtained with this procedure the resulting relative error in the 
frequency determination is conservatively estimated at less than ±1% . 

Magnitude and Phase Determination 

Once the drive frequency, f d , was determined, the magnitude and phase of each 
signal was computed by evaluating the DTFT (Discrete Time Fourier Transform) at 
the drive frequency. For the length N , real sequence x[n], the DTFT of x[n], denoted, 
X(e) UJ ) evaluated at the frequency tv is given as in Oppenheim by: 

X(en = E x{n)e-i“ n . 

71 = 0 ’ 

If the physical drive frequency, f d , has dimensions of cycles per unit time and the sam- 
pling period is T, then u in the above expression is defined by u; = 2tt f d T. Denoting 
the real and imaginary part of X(e^) by a and b , respectively, the amplitude, A, and 
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phase angle, 9, of the signal at the drive frequency is then found using: A = yj a 2 + fe 2 
and 9 = arctan(b/a). 

A/D Skew Correction 

As discussed in Chapter 2 there is a known, fixed, inter-channel time delay (A/D 
skew) between the various signals. For the situation of interest here, a sinusoidal 
signal at a known frequency, f d , this effect was easily correctable. For a delay time, 
t d , the phase lag due to the delay is given by 69 = 2 rf d t d . The phase lag was 
calculated in this way for each channel and then subtracted off to restore all channels 

to a common time base. 


4.5 Experimental Results 

4.5.1 Quasi-steady Model 

The simple quasi-steady model, as expressed by Equation 4.4, requires two parame- 
ters, compressor slope ^ and equivalent duct length, C c . These parameters were 
identified, i.e. their numerical values were determined, by varying them to obtain a 
good fit between the calculated and experimentally measured frequency response. 

Since the equivalent compressor length (fluid inertia) is a purely geometric pa- 
rameter, which does not vary with the time mean flow, a common value was used to 
fit the frequency response data for all four operating points. The compressor slope 
in general is expected to vary with the operating point and so different values were 
allowed for each operating point. The values obtained were C c = 0.75 m for all flows 


90 


an( j _ i _ 5 ; 0.80, 0.38, —2.1 corresponding to <f> = 0.10, 0.12, 0.15 and0.20 re- 
spectively. Independent estimates for these parameter values, based purely on known 
compressor geometry and steady state data, are consistent with these identified val- 


ues, as detailed subsequently in this section. 

The frequency response plots comparing the measured admittance with the the- 
oretical response using these parameter values are shown in Figure 4-4. The “Bode 
Plot” format (log magnitude versus log frequency aqd phase angle versus log fre- 
quency) used here was selected as the most relevant for control design purposes as it 
gives a direct indication of the required compensation and associated gam and phase 
margins. Theory and measurement are seen to agree well over the higher frequency 
range but not at lower frequencies. 

The poor agreement at low frequencies can be attributed to the fluctuations in ro- 
tor speed which occur at the lower frequencies. This is demonstrated by the improved 
agreement at low frequency shown in Figure 4-5 where the data has been corrected, 
as described below, to account for the measured rotor speed fluctuations. At higher 
forcing frequencies, the inertia of the rotor tends to maintain a constant rotational 
speed, so the correction has little effect. 

To obtain this plot, the measured rotor speed fluctuations were accounted for by 
assuming that the compressor responds linearly (so that superposition applies) and 
quasi-steadily to rotor speed fluctuations. Using the measured speed fluctuation pha- 
sor, M r , and sensitivity, (obtained by the procedure described in Section 4.3.2) 

the pressure rise fluctuation phasor due to speed fluctuation alone, V>M r , waS cal ' 
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culated using, Mj = ^M r . The effect of this was considered as an additional 
forcing term acting on the compressor duct which was then superimposed, using 
phasor addition, to obtain the net driving pressure. With the speed fluctuations ac- 
counted for, overall agreement between the simple quasi-steady model and measured 
data is quite good. Unless otherwise noted, in the remaining data presented in this 
chapter, speed fluctuations are accounted for in this manner. 

A number of qualitative features shown by the data in Figure 4-5 warrant further 
comment. The magnitude of the admittance is seen to be essentially monotomcally 
decreasing with frequency, in agreement with theory. This indicates that, over the 
tested range of frequencies, there are no unexpected resonances. Thus, if any higher 
order dynamics have been neglected by the simple model over this frequency range, 
they are well damped and unlikely to cause closed loop stability problems. For the flow 
coefficients where the compressor slope is positive, (<j> = 0.10,0.12,0.15) the phase an- 
gle characteristics are seen to monotonically increase while for the operating point on 
the negatively sloped portion of the characteristic, <f> = 0.20, the phase angle decreases 
monotonically. This difference in behavior can be understood by considering Equa- 
tion 4.4 in which the phase of the admittance is given by arctan( 7 ^g), showing that 
the algebraic sign of the phase depends upon the slope of the compressor characteris- 
tic. Alternatively, this change is indicative of the migration of the system’s pole from 
the right to the left half of the complex plane as the compressor slope changes sign 
with increasing flow. Finally it is noteworthy that the “break frequency”, where the 
relatively constant portion of experimentally determined magnitude of the frequency 
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response begins to decrease (“roll-off”), occurs at a normalized frequency, of unity, 
for all the flows, in agreement with the theoretical curves. Physically, the break fre- 
quency occurs at the transition between quasi-steady behavior with negligible fluid 
inertia effects at low frequency to inertia dominated behavior at high frequency. The 
agreement between experiment and theory indicates that the relevant time scale is in 
fact set by the fluid inertia and compressor slope as defined by Equation 4.10. 

As previously noted, the parameter values for this model were also estimated in- 


dependently to provide a check on the reasonableness of the values identified from the 
frequency response data. An estimate of the equivalent length based purely on geom- 
etry gave a value of C c = 0.63(m) compared with the identified value of 0.75(m). Dif- 
ferentiating a fourth order polynomial fit of the steady state compressor speed line to 

obtain estimates of the compressor slope gives values of ^ = 1-5, 0.86, -0.05, -4.1 

in order of increasing flow coefficient, compared to the corresponding identified val- 
ues 0 f _ x.5, 0.80, 0.38, -2.1. Thus, considering the relative crudeness of the 
independent estimates, the identified and estimated parameters appear to agree well, 
lending additional support to the conclusion that the important physics are properly 
captured by the model. Note that without this confirmation, the compressor might 
still be adequately modeled by a first order system with the same functional form of 
Equation 4.4, but be governed by an entirely different set of physics. The difference 
between the identified and estimated values indicates that, for detailed control design 
purposes, it may be worthwhile to perform the system identification as done here 


rather than to rely on a priori estimates. 
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4.5.2 Unsteady Model 

For basic control design purposes, the quasi-steady model appears to provide an ad- 
equate representation of the compressor as shown in the preceding section. However, 
in many compression system dynamic studies, for example, Greitzer [17], Fink [14] 
and Bons [5] the existence of some unsteady lag in the compressor response has been 

assumed. 

The experimental frequency response data obtained here was therefore examined 
to determine if there were direct evidence of such lagging behavior. To do this, 
the theoretical quasi-steady impedance, Equation 4.5, was first compared with the 
observed data. The previously identified values of the slope and equivalent length 
were used to compute the theoretical response. For this comparison, the impedance 
rather than the admittance was used, as it provides a more sensitive and diagnostic 
measure of the unsteady effects, as discussed in Section 4.2.2. 

The result is shown in Figure 4-6 where the real and imaginary parts of the 
impedance are plotted against the reduced frequency, w„, based on flow through time 
as defined by Equation 4.13. As discussed in Section 4.2.4, u; u would be expected to 
provide the relevant frequency scaling for unsteady aerodynamic effects. 

Two observations are made regarding Figure 4-6. First, the real part of the 
impedance, which is predicted by the quasi-steady theory to be a constant, numer- 
ically equal to the slope of the compressor pressure rise characteristic, but opposite 
in algebraic sign, in fact varies with frequency. In particular, for the three oper- 
ating points on the positively sloped portion of the characteristic, the real part of 
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the impedance goes from a negative value at low frequency to a positive value at 
high frequency. The second observation is that the imaginary part of the normalized 
impedance, Im{Z(ju)/4>}, which is predicted to be proportional to the frequency 
(with the constant of proportionality independent of operating point) actually varies 
with the time mean flow coefficient. 

To determine whether a relatively simple model of the unsteady compressor re- 
sponse could account for the observed deviations from quasi-steady behavior, the 
experimental data was compared to the theoretical impedance of the unsteady model 
described by Equation 4.9. Along with the the compressor slope and equivalent length 
of the quasi-steady model, the unsteady model contains a single additional parame- 
ter, r L , the ratio of lag time and flow through time. Leaving the original parameters 
unchanged, the time lag parameter, r L , was set to a value of 1.4 which gave the best 

overall agreement between theory and experiment. 

The results are shown in Figure 4-7, which compares the experimental impedance 
with the unsteady model predictions. It can be seen that the unsteady model predic- 
tions are consistent with both the previously observed trends, that is: the frequency 
dependence of the real part of the impedance and the mean flow dependence of the 
slope of the imaginary part of the impedance. Although the agreement is not perfect, 
particularly for the real part of the impedance at the highest frequencies and flow 
rates, it can be concluded that, overall, the observed behavior appears consistent 
with the hypothesized unsteady loss model. Considering the probable complexity of 
the loss development process, it is not surprising that such a crude model, which 
lumps all of the loss development into a single time lag, is not adequate for accurate 
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prediction as reduced frequency approaches unity. The unsteady loss model, based 
on flow through time accounts well for the dependence of the imaginary part of the 
impedance on the time mean flow rate. 

This dependence on mean flow had also been observed by Abdel-Hamid [11 who 
accounted for it as a flow dependent equivalent compressor length. Since equivalent 
length is a purely geometric quantity this explanation is not physically plausible. 
In contrast, the simple unsteady loss model appears to offer a physically reasonable 
explanation for the observed trend. 

4.5.3 Relevance of Unsteady Behavior for Control Design 

Finally, from a control design perspective, it is interesting to examine the Bode Plots 
of the admittance computed with the unsteady loss model. This is done in Fig- 
ure 4-8, which compares the experimentally measured admittance with the theoreti- 
cal response using the unsteady loss model and the previously identified parameters. 
Comparison with the quasi-steady model results given in Figure 4-5 shows that the 
unsteady model provides a marginally better match to the experimental data, partic- 
ularly the phase. This would reduce the required phase margin in the control design. 
Thus, for detailed design and fine tuning, use of the- unsteady loss model may be 

warranted by the additional fidelity which it provides. 

The more significant aspect of the observed unsteady behavior is that it suggests a 
physical upper bound on the required control bandwidth which does not exist for the 
quasi-steady model. Specifically, the data shows that the real part of the impedance 
becomes positive at high frequencies. This implies that the system cannot support 
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sustained oscillations above some finite high frequency limit. 


4.6 Summary and Conclusions 

The results of the compressor model validation study are summarized as follows. 

• The experiments show that this class of compressor can be well represented 
using a lumped parameter, linear model for small sinusoidal disturbances over 
the frequency range of interest. 

• The quasi-steady theory provides a good approximation which is useful for 
control design purposes. 

• For the particular drive arrangement used for these experiments fluctuations in 
rotor speed occurred at low frequencies. Effects of the speed fluctuations can 
be accounted for by superimposing the quasi-steady response to measured rotor 
speed fluctuations. 

• A simple unsteady loss model gives good agreement with data and can be used 
to account for departures from quasi-steady behavior at intermediate reduced 

frequencies. 

• The experimental data shows that when operating on the positively sloped 
portion of its characteristic, at reduced frequencies (based on flow through time) 
much less than unity the compressor impedance (as defined by Equation 4.2) is 
negative. As the reduced frequency approaches unity, the compressor impedance 
becomes positive. 
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• The observed high frequency change in the algebraic sign of the compressor 
impedance, which is not predicted by the quasi-steady model, indicates that 
the quasi-steady model will be conservative in terms of predicting the required 


control system bandwidth. 
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NOTE: NOT TO SCALE 


Figure 4-3: Equipment configuration for frequency response measurement 
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Chapter 5 


Limits to Compression System 
Stabilization 


Previous experimental and analytical work ([12, 20, 40] has shown that feedback can 
be used to extend the stable (surge free) operating range of a compression system. 
These studies have also indicated that, for the particular schemes investigated, there 
were limits to the range of system parameters (flow coefficient, B parameter, com- 
pressor slope etc.) over which the control remained effective. Motivated by these 
findings, an investigation has been carried out to determine what fundamentally lim- 
its the ability to stabilize a compression system. Relevant control theory is developed 
and applied to illuminate the nature and extent of these limitations. 
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5.1 Nominal Stabilization 


The most basic requirement of the feedback control used to extend the operating 
range of the compression system is that the resulting closed loop system be stable to 
infinitesimal perturbations, that is, linearly stable. In designing such a control, the 
true system is idealized with a model. If this model, when coupled with the feedback 
control, is linearly stable, the system is said to be nominally stabilized. Because of the 
approximations involved in the modeling process, nominal stability does not ensure 
that the actual system with the same control will also be stable, but, from the point 
of view of performing a rational design, nominal stability is a minimal requirement. 

Additional measures must be taken to ensure that the implementation shares the 
stability properties of the nominal model. This is the stability robustness problem 
of control design and will be addressed subsequently. In this section the more basic 
requirement of achieving nominal stability for the basic compression system model 
will be addressed. 

5.1.1 Theoretical Background 

We will consider linear, time invariant, SISO (single input single output) system 
models of the form: 


x = Ax + btx (5-1) 

y = cx, 
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where, the vector x denotes the system state, u is the scalar control input, y is the 
scalar output and A,b,c, are appropriately dimensioned constant matrices. It is well 
known (see for example [24, Theorem 5.2]) that the system, given by Equation 5.1, 
can be stabilized using a linear time invariant, finite dimensional, feedback control if 
and only if it is both stabilizable and detectable (The terminology here follows [24] but 
is fairly standard.) If the stronger conditions of controllability and observability are 
met, then not only can the system be stabilized but, in addition, the closed loop poles 
of the system can in fact be placed arbitrarily in the left half plane. Complex poles 
must, however, be placed as complex conjugate pairs. Such stabilizing controllers can 
readily be synthesized using standard design techniques and numerical algorithms. 

For our purposes therefore, the question of whether a compression system can 
be nominally stabilized is reduced to determining under what conditions its linear 
system model is stabilizable and detectable. These properties, in turn, depend upon 
the particular choice of actuators and sensors, which determine b and c respectively 
in Equation 5.1, as well as the homogeneous system dynamics expressed by the system 

A matrix. 

5.1.2 Application of Theory to Compression Systems 

The theoretical results outlined in the preceding section will now be applied to the 
question of nominally stabilizing a compression system using a single actuator and 
and a signal sensor. Note that, u and y in Equation 5.1 are scalars. Correspondingly, 
b and c are column and row vectors respectively. Consider the linearization of the 
basic compression system model with two states given by Equation 3.25 with an as 
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yet unspecified sensor and actuator. The system is described by: 
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(5.3) 


which is of the form of Equation 5.1. The (as yet unspecified) input and output ma- 
trices, b, and c, are determined by the particular choice of actuator and sensor. The 
compression system can be stabilized if it is controllable and observable 1 . Applying 
standard rank tests on the controllability and observability matrices shows that this 
second order system is controllable if and only if 


rank( 


b Ab 


) = 2 


and it is observable if and only if 


ran 


k( 


cA 


) = 2 . 


Before considering special cases, i.e. particular choices of b and c, it is noteworthy 
that if the open loop poles are complex, i.e. have non-zero real part, then these 
conditions will be met by any actuator/sensor pair. This can be seen by noting that 


i As previously discussed in Section^. 1.1 it need only be stabilizable and detectable. However 
the more conservative requirement of controllability and observability is somewhat more straight 
forward to check and will therefore be used here. 
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if, 

rank( b Ab ) < 2 > 

then, b and Ab are linearly dependent, which implies that b is an eigenvector of A. 
But since A and b are real, this is impossible unless the eigenvalue is also purely real. 
A similar argument holds for the observability matrix. Thus, for any operating point 
and set of characteristics for which the open loop system A has complex poles, a 
stabilizing control can be designed for any choice of actuator and sensor. Limitations 
on the allowable range of parameters over which a particular choice of actuator and 
sensor can be used to nominally stabilize a compression system will, therefore, only 
occur for parameter values which result in purely real open loop poles. 

From the rank tests it is seen that limitations occur if b or c are aligned with 
an eigenvector of A or A T respectively. Fortunately, for most specific choices of 
actuator and sensor, this cannot occur for any physically meaningful value of the 

system parameters. 

As a practical example, consider a variable plenum exit throttle actuator with 
plenum pressure sensing. This corresponds to Pinsley’s [401 experiment. In this case 

we have: r -i 


0 

(5.4) 

1 


~B da J 


0 1 

(5.5) 
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For this two dimensional case, if, 


det( 


b Ab 


)/ 0 , 


then, 


rank( 


b Ab 


) = 2 - 


Similarly, 


det( 


)^ 0 , 


cA 


implies that, 


rank( 


C 

CA 


) = 2 - 


From which we obtain the condition for controllability. 




( 5 . 6 ) 


which will be satisfied for any finite value of B so long as the mass flow through the 
throttle valve (with pressure drop held constant) is a strictly increasing function of 
the valve area, which will always be the case. Similarly, the condition for observability 
can be obtained : 

I + 0, (5.7) 

which again holds for any finite value of B. Thus it has been shown that, for all 
physically reasonable values of the system parameters, the compression system can be 
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nominally stabilized using a variable area throttle at the plenum exit as the actuator 
and a plenum pressure sensor. It should be noted that limitations in terms of the 
maximum B-parameter which can be stabilized exist if the control law is restricted 

i 

to a simple proportional control (see for example Pinsley [40]). In order to nominally 
stabilize this system for any value of B-parameter a dynamic compensator (see [24, 
Chapter 5]) will be required. 

Similar calculations can be done to show that many other reasonable choices of 
actuator and sensor also yield systems which are both controllable and observable. 

This calculation does point out, however, that as B becomes infinitely large the 
system would in the limit no longer be controllable or observable. It should then be 
anticipated that large but finite values of B will be difficult to control with this choice 

of actuator and sensor. 

The main point which has been established is that the nominal stabilization of the 
basic compression system does not generally set a fundamental limit to control. This 
is an important theoretical result because it points out that, as a means of comparison 
between various actuation and sensing strategies, the ability to nominally stabilize 
the system provides little discrimination. 


5.2 Bounded Actuation 

In the analysis presented in the previous section, it is assumed that there is no bound 
on the allowable magnitude of the control signal u. This section considers the more 
realistic situation where the control is bounded. 
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5.2.1 Theoretical Background 

Suppose that a stabilizing feedback control is connected to the linear system described 
by Equation 5.1. Then, with no bound on the control signal, any initial state will 
be driven back to the origin by the controller. This is not the case if the control 
magnitude is bounded, ie., u mia < « < ^max- With the control bounded and the 
plant open loop unstable, the system cannot be restored to equilibrium for all initial 

states. 

The set of initial states which can be restored to equilibrium subject to the con- 
straint on the control magnitude includes only a limited region in the overall state 
space. This region (which will be more precisely defined subsequently) will be referred 
to as the recovery region 2 . For a control system to be effective, in practical situations, 
the recovery region cannot be too small. 

To make the previous notions more precise, consider the linear discrete time 
system 3 : 

x n+ j = A,jX n + b(ju n , (5-8) 


which describes the behavior at the sample points of the continuous time system, 
Equation 5.1 driven by a discrete time controller 4 . Let the control sequence be 


2 The concepts of recovery region and degree of controllability were developed ^^pendendy m 
the course of this thesis research and were only later found to coincide with Schmitendorf s [42] 

work. Schmitendorf’s terminology will be used here. t . , pvpr th analvsis 

3 A similar analysis can be performed for the continuous time control case, however, the analysis 

iS input to the plant la held constant over th. sample interval. 
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bounded by the constraint, 


u min — u » — u niax * — 3... (5-9) 

Following Schmitendorf 5 [42], 7 l{nT), the recovery region at time nT is defined as the 
set of all initial states, x(0), which can be driven back to the origin using a control 
sequence u u u 2 , •••,«„ satisfying the constraint equation 5.9 in the time interval, 

0 < t < nT where T is the sampling period of the discrete time control. The degree 
of controllability g is then defined as the shortest distance between the origin and the 
boundary of the recovery region, i.e., 

g = inf { ||x|| : x G 7Z c (nT)}, (5-10) 

where, 7l c {nT) is the complement of the recovery region, 7 Z(nT), and inf{} is defined 

as the greatest lower bound on the indicated set. 

For a system with two states, for example the basic compression system model, 
the degree of controllability is simply the radius of the largest circle centered at the 
origin (equilibrium point) which can be contained in the recovery region. In general, 
it is the radius of the largest hypershere centered at the origin which can be contained 
inside of the recovery region. Every initial state inside this circle, or hypersphere for 

the general case, can be returned to the origin. 

For open loop unstable systems, the degree of controllability approaches a finite 

5 Schmitendorf defined these terms terms for the continuous control case and his terminology has 
been adapted to the discrete time case treated here. 
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limit as the number of steps, n increases to infinity. That is, due to the bounded 
actuation, some initial states cannot be driven back to the origin no matter how long 
a time we allow. If for some choice of actuator and set of system parameters the 
magnitude of g(nT) as n - oo is small, with respect to the expected magnitude of 
initial conditions, (assuming an appropriately scaled set of state variables is used) an 
actual system using this actuator cannot be expected to be effective. In such a case, 
the slightest perturbation in the initial system state could place the system into a 

region which could not be returned to equilibrium. 

The magnitude of the degree of controllability thus provides a figure of merit for 
assessing the relative effectiveness of various actuators. Since the analysis assumes a 
linear plant model, it is most appropriately employed as a means for distinguishing 
ineffective actuators, as opposed to determining the most effective actuators. This 
can be understood by considering that for an ineffective actuator the degree of con- 
trollability will be small, and even small initial perturbations, for which the linearized 
dynamics are accurate, cannot be returned to equilibrium. If the degree of control- 
lability is too small it is, therefore, unlikely that the system will work in practice. 
On the other hand, it is important to note that one cannot assume that an actuator 
with an extremely large degree of controllability would actually achieve this radius as 

nonlinear effects may become important. 

Finally, it is noted that the recovery region is the largest possible region that can 

be returned to the origin without violating the control constraints. It is thus the best 
that can be done with any control law, linear or nonlinear, and any choice of sensor. 
The degree of controllability thus provides a means of comparing the effectiveness of 
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various actuators which is independent of the choice of control law or sensor. 

5.2.2 Numerical Results for Compression System 

A numerical algorithm to calculate the boundary of the recovery region and the degree 
of controllability was developed and implemented. The method taken closely follows 
the continuous time approach described in [42], but has been adapted for a discrete 
time control (see Appendix B). To illustrate the basic concepts for a practical case of 
interest, recovery regions, 1l(nT), calculated for various values of nT for an unstable 
compression system with a plenum bleed valve actuator, are shown in Figure 5-1. The 
maximum fluctuations in bleed flow are constrained to be less than ±1% of the mean 
mass flow. In the planar case shown here, the degree of controllability at time nT is 
defined as the radius of the largest circle centered at the origin which can be contained 
inside of the recovery region at time nT. It can be seen that, for increasing values of 
the recovery time, the size of these circles and, therefore, the degree of controllability, 
approaches a limiting value as expected for an open loop unstable system. 

It is apparent that the degree of controllability is dependent upon the scaling of 
the state variables and input. For a clear physical interpretation, the degree of con- 
trollability should be based on an appropriately scaled set of state variables. That 
is, a perturbation of unit magnitude should have the same physical significance re- 
gardless of its direction in state space. To this end, scaled state variables are chosen 
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where ~fc T P , and J b are the values at the equilibrium operating point, and ?Wx 
is the maximum bleed flow. The resulting scaled recovery region at time nT = 2, ie. 
two Helmholtz periods, and the corresponding circle whose radius defines the degree 

of controllability, are illustrated in Figure 5-2. 

With this scaling, the degree of controllability for the plenum bleed valve is plotted 

as a function of system B parameter and compressor slope in Figure 5-3 for nT = 10. 
Additional computations using larger values of nT showed that no significant increases 
in the degree of controllability would be obtained by further increasing the recovery 
time. To help interpret the results, if the degree of controllability were equal to 
0.01, then, with this scaling the maximum perturbation in initial plenum pressure or 
compressor flow must be less than 1% of the equilibrium value. Figure 5-3 shows that 
the degree of controllability decreases asymptotically towards zero as the B parameter 
and or compressor slope increases. 

The analysis quantifies a trend which can be anticipated from the following phys- 


as: 

Z\ = 

X2 = 

with the input scaled as: 

u = 
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ical considerations. As the compressor slope becomes increasingly positive, the force 
accelerating the compressor duct fluid for a given flow perturbation increases, and 
consequently so does the growth rate of the instability. As B increases, the relative 
compliance of the plenum increases which decreases the effectiveness of the bleed 
valve in influencing the compressor. Since the control is constrained, both of these 
trends require the allowable region of initial conditions to decrease, he., the degree 
of controllability is reduced. This calculation shows that even though the system 
can be nominally stabilized, the allowable set of initial conditions when the control 
is bounded becomes so restrictive that for practical purposes the system cannot be 

stabilized. 

It is revealing to compare this to the case of an actuator close coupled to the 
compressor exit. The actuator (which could be a valve) is idealized here as a device 
which can produce a desired pressure rise (or drop) which is simply added to that 
of the compressor to obtain the overall pressure rise of the compressor and actuator. 
The resulting degree of controllability, assuming that the actuator could produce 
a maximum pressure perturbation of ±1% of the time mean compressor pressure 
rise (and otherwise using the same parameter values as for Figure 5-3) is shown in 

Figure 5-4. 

In contrast to the plenum bleed actuator case, the degree of controllability for 
the close-coupled actuator does not asymptotically approach zero. On this basis, the 
plenum bleed valve can be rejected as being the less favorable option. It would be 
anticipated that bounded actuation will similarly severely limit the effectiveness of 
many control schemes which might appear viable if only nominal stabilization were 
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considered. 


5.3 Limitations Imposed by Model Uncertainty 

The approximations used in developing the simple lumped compression system model 
are valid over only a limited range of frequencies. Above this range of frequencies we 
would expect that the behavior of the system would depart significantly from that 
of the model. It is well known that these unmodeled high frequency dynamics can 
potentially result in instability when a feedback control loop is closed around the 
system, even though the nominal model with high frequency dynamics neglected is 
theoretically stable. 

To prevent such instability, the feedback control must be properly designed. This 
is one aspect of what is known as the stability robustness problem and there is a large 
body of literature treating this subject. The basic approach to ensuring stability in 
the face of uncertain high frequency dynamics, which is common to many of these 
methods, is to impose limitations on the closed loop system bandwidth. 

For an open loop unstable system, which is the main concern here, it is not 
clear to what extent this prescription can be followed. That is, can the closed loop 
bandwidth in fact be made arbitrarily small for such systems? Intuitively, one would 
expect problems to arise as the time scales of the instability approach the time scales 
at which the model becomes inaccurate. 

In this section it is shown that this is in fact the case. The problem will be quan- 
tified as a design tradeoff between stability robustness to high frequency model error 
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and nominal performance. This is a fundamental control problem which apparently 
has not been previously addressed in this fashion and so the analysis is of some gen- 
eral interest. Focus will, however, be primarily maintained on the implications for 
the compression system, to which end specific numerical results will be given for an 

illustrative case. 

The presentation is organized as follows. Section 5.3.1 provides a brief summary of 
relevant background concepts and theory. Section 5.3.2 apphes the theory to formu- 
late the desired quantitative design tradeoff and Section 5.3.3 gives some numerical 
results for the compression system application. 


5.3.1 Theoretical Background 

Closed Loop Transfer Functions 

A frequency domain (Laplace Transform) description of the linearized system is used 
for the analysis. A block diagram of the system to be analyzed is shown in Figure 5- 
5. The notation used here follows [15) and is fairly standard. Laplace transformed 
quantities will be distinguished by the use of s, or ju, if evaluated on the imaginary 
axis, as their argument. Referring to Figure 5-5, P(s) is the plant transfer function 
representing the compression system and F(s) is the transfer function of the feedback 
system used to stabilize the loop. The Loop Transfer Function, I(s), is defined by: 

L(s) = P(s)F(s). (5 ' U) 
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The Sensitivity Transfer Function , 5(a), and Complementary Sensitivity Transfer 
Function , T(s), are defined respectively by: 


5(s) 


1 

1 + L{s) 


T(s) 


1 + L(s ) 


(5.12) 

(5.13) 


The two are related by the identity: 


S(s) = 1 T(s) 


(5.14) 


By expressing the closed loop response in terms of the Sensitivity and Comple- 
mentary Sensitivity Transfer Functions, the role played by these functions in dictating 
the overall performance in the face of inevitable disturbances and sensor noise is made 
clear. Referring to the block diagram, Figure 5-5, the response of the plant output, 
y (s ), to disturbances at the plant input, d,-(s), disturbances at the plant output, d 0 (s), 
sensor noise, n(s), and reference command, r(s), is given by: 


y(s ) = S(s)d 0 {s) + P{s)S(s)di(s ) + T{s)n(s) + T{s)r(s) (5.15) 


The response of the controller output, u(s) to these inputs is given by: 


u(s) = -S(s)F{s)d 0 (s) - T(s)di{s) - S{s)F(s)n(s) + S{s)F(s)r{s) (5.16) 


121 


Stability Robustness 


Let L(s) be the nominal loop Transfer Function, given by the approximate model, 
and L'(s) be that of the actual system. It will be assumed that the nominal system 
is closed loop stable and that the true and nominal plant models are related by: 

L'{s) = L{s)(l + e{s)) . ( 5 - 17 ) 


where, e{s) is a stable transfer function whose magnitude can be estimated but is 
otherwise unknown. This is known as a multiplicative representation of unstructured 


uncertainty, see for example [50]. The relative error between the nominal and true 
model is given by: 


L'(s) - L(s) 

L(s) 


|£(s)|. 


(5.18) 


Using the Nyquist Stability Criterion, it can then be shown, see for example [50], that 
the actual closed loop system will also be stable so long as the nominal Complementary 
Sensitivity Transfer Function evaluated along the entire ju axis satisfies: 


\T{joj)\ < l/|e(jw)| 


(5.19) 


Thus, the Complementary Sensitivity Transfer Function gives an indication of sta- 
bility robustness to unstructured multiplicative uncertainty as well as governing the 
nominal closed loop response to disturbances and noise. 
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Frequency Domain Design Specification 

Ideally, the closed loop system would exactly follow the input command in the face 
of the various disturbances and noise sources. From Equations 5.15 and 5.16 it can 
be seen that, for the system output, y(s), to follow the command and to keep the 
response to disturbances small, we should have \T(s)\ « 1. On the other hand, to keep 
the response of the output to measurement noise small we should have \T(s)\ < 1. 
From Equation 5.19 it is seen that for stability robustness we also require \T(s)\ < 1 
along the portions of the ju axis where the modeling error is large. 

If the commands and disturbances are limited to low frequencies and the sensor 
noise and model uncertainty are largely confined to high frequencies, then reasonably 
good performance can be achieved if \T(s)\ can be held close to unity at low frequen- 
cies and then made to decrease (roll off) above some high frequency cutoff. Such a 
specification is illustrated in Figure 5-6. It is demonstrated in the following section 
that, unfortunately, specifications of this sort cannot necessarily be satisfied for an 
open loop unstable plant. 

Frequency Domain Integral Constraints 

Freudenberg and Looze, [15], extending the approach pioneered by Bode, [3], have 
shown that the achievable closed loop transfer functions, S(s) and T{s), for plants 
with right half plane (unstable) poles, zeros, and time delays cannot be shaped arbi- 
trarily. They have expressed these limitations in terms of integrals of the logarithm of 
the magnitudes of the S(s) and T(s) evaluated along the ju axis. These will be called 
frequency domain integral constraints here. As a direct consequence, reductions in 
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the magnitude of S(s) or T(s) over one range of frequencies can only be achieved at 
the expense of increased magnitude at other frequencies. 

5.3.2 Theoretical Stabilization Limitation 

The general limitations implied by Freudenberg and Loose's frequency domain inte- 
gral constraints can be used to obtain quantitative estimates of the price pa.d for 
high frequency model uncertainty for the basic compression system. Suppose that we 
wish to stabilize the compression system at an open loop unstable operating point. 
To obtain such a bound, let us assume that, above some high frequency cutoff, u c , the 
magnitude of the relative error between the nominal and actual frequency response 

is bounded above by £ m • That is. 

k(i w )l < £ m f° r M > Uc e m > 1- (5.20) 

Further, to satisfy the stability robustness condition, Equation 5.19, the closed loop 
is constrained to satisfy: 

|T(;cv)| < l/|e„l for M > Wc ( 5 - 21 ) 

Let p = a + jb denote the most unstable of the two open loop poles of the nominal 
compression system model. (That is the one with the largest real part.) If at the 
operating point of interest the system has two complex conjugate poles either one can 
be selected and called p. Let «, i = l,n. denote the right half plane zeros, if any, of 
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L(s). 

With these assumptions, Freudenberg and Looze’s results can be extended to show 
that the complementary sensitivity has a peak magnitude on the ju axis, denoted by 
lirCHIU, which is lower bounded by: 

||T(;u,)|U > (5-22) 


where 

8{u c ) - arctan + arctan ^ ( 5 - 23 ) 


and, 


nn> ( z ' + p) 

lh=i 


1 


; n z > 0 

; n z = 0 


(5.24) 


Because k 2 > 1 and by assumption £ m > 1, Equation 5.22 shows that the peak 
complementary sensitivity must be greater than unity. Since \T{ju) < 1| for \u\ > u c 
the peak must occur at a frequency below w c . The true peak value is always at least 

as large as the estimate given by Equation 5.22. 

It is noted that when the system has right half plane zeros, k z is greater than unity, 
so that the situation is always made worse by the presence of these zeros. Whether or 
not there are any right half plane zeros depends on the choice of actuator and sensor. 
Inspection of the denominator of k 2 shows that the peak complementary sensitivity 
will become extremely large when there is a near pole zero cancellation in the right 
half plane. Thus, the actual location as well as the mere presence of right half plane 


zeros and poles is important. 
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Due to the relationship between T(s) and S{s) given by Equation 5.14, the sen- 
sitivity transfer function must also have a peak value greater than unity when ever 
||T(ju>)||oo >■ 1. This is shown as follows. Since, 

||S(ju;)||oo = P — ^CMIloo 

> ||T(i«)lloo - 1 

it can be seen that ||S(j«)l|oo wil1 be substantially greater than unity whenever 
||T(j'w)||oo > 1- A geometrical construction showing the function 6{w e ) as an an- 
gle in the complex plane is shown in Figure 5-7. From Equation 5.22 it is apparent 

that, for ||T(ja>)||oo to be small, we must have 0(w e ) « 

From the geometrical construction shown in Figure 5-7, it can be seen that, for 
a fixed « e , as either the growth rate of the instability (real pole case) or the natural 
frequency of the instability (complex conjugate case) or both become large, with 
respect to u ; c , 0(o> c ), will become small, implying a large ||T(jw)||oo-. Equation 5.22 can 
thus be used to quantify the intuitively anticipated trend that acceptable performance 
can only be achieved when the instability growth rate is small with respect to the 
high frequency limitation on the range of model validity. 

5.3.3 Numerical Results for a Compression System 

To appreciate the severity of the limitation imposed by the bound given by Eq 
tion 5.22 it is useful to compute numerical values for the basic compression system 
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model using some typical parameter values. The results of performing such a cal- 
culation are shown in Figure 5-8 for the case of a dose-coupled valve actuator with 
measurement of compressor mass flow rate. For this pairing of actuator and sensor 
there are no zeros in the right half plane so that k, = 1 in Equation 5.22. In this 
figure, the lower bound on ||T(jw)|U is plotted “ ». f “ nction ° f compressor slope 
and u> c /u k , the ratio of the high frequency cutoff to the system undamped natural 
frequency (Helmholtz Frequency). The maximum relative modeling error is held fixed 
at £ m =2.0 and the parameter B is held fixed at B = 2. 

For the lower values of compressor slope, m c , the open loop poles are complex 
conjugates and, as shown in Figure 5-8, the peak complementary sensitivity can be 
made reasonably close to unity even for relatively small values of u e /u h . That is, the 
model need only be good up to the Helmholtz frequency. For the steeper compressor 
slopes or larger values of B, the figure shows that we must have uju h greater than 
8 to reduce the lower bound on ||TtfuO||oo to a reasonable value. In this case, the 
unstable open loop poles are real, the important time scale has become the growth 
rate of the fastest unstable pole, and the model must be accurate for these fast time 

scales. 

For other choices of sensors and actuators, there are also right half plane zeros so 
that k 2 > 1 and, from Equation 5.22, the limitation would be expected to be more 

severe. 

High frequency modeling uncertainty, therefore, imposes an effective limit on the 
range of compressor slopes and B which can be robustly stabilized while maintaining 
acceptable performance. A means for estimating this bound has been derived and 
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numerical results utilizing this procedure have been presented to demonstrate the 
severity of the limitation for an illustrative case. The existence of right half plane 
zero increases the severity of these limitations and properly choosing actuators and 
sensors to eliminate such zeros is an important design consideration. 

5.4 Summary 

In this chapter, important theoretical limitations to stabilization of compression sys- 
tems have been examined in detail. For many choices of actuators the compression 
system can be nominally stabilized. In the more realistic situation where actuators 
are bounded and high frequency modeling errors are considered the ability to stabilize 
the system may be severely restricted over the parameter ranges of practical interest. 
These limitations are more severe for some actuator and sensor pairs than others and, 
therefore, a proper choice of actuator and sensor is essential. This finding is further 
amplified and explored in detail in Chapter 6 where a systematic evaluation of various 
actuator sensor pairs is performed. 
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Figure 5-1: Recovery regions for plenum bleed actuator as function of recovery time 
(illustrative case) 
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Figure 5-2: Recovery region for plenum bleed actuator showing circle whose radius 
determines degree of controllability (illustrative case) 




Figure 5-3: Degree of Controllability as function of system parameters for plenum 
bleed actuator 
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Figure 5-4: Degree of Controllability as function of system parameters for close 
coupled actuator 



Figure 5-5: Block Diagram for defining closed loop transfer functions 
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Figure 5-6: Typical frequency domain performance specification for the complemen 
tary sensitivity transfer function 




Figure 5-7: Geometrical Construction for the weighting function 0(w c ) a)real pole 
b)complex conjugate pole 
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Figure 5-8: Lower Bound on peak complementary sensitivity for close-coupled valve 
with mass flow measurement 
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Chapter 6 


Evaluation of Alternative Control 
Strategies 


6.1 Introduction 


Almost all research to date on the use of feedback stabilization of compression sys- 
tems has been aimed at proof of concept of the technique, and only a few of the many 
possible actuators and sensors have been considered. In this chapter, we therefore 
present the first systematic definition of the influence of sensor and actuator selection 
on increasing the range of stabilized compressor performance 1 . The results show that 
proper choice of sensor as well as actuator crucially affects the ability to stabilize 
these systems and that, overall, those actuators which are most closely coupled to the 
compressor (as opposed to the plenum or throttle) are most effective. In addition, 


!The work of this chapter was originally 
International Gas Turbine and Aeroengine 


presented in a similar form by Simon et al [46] at the 
Congress and Exposition, Cologne, Germany June 1-4, 


1992. 
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the source of the disturbances driving the system (for example, unsteady compressor 
pressure rise or unsteady combustor heat release) has a strong influence on control 
effectiveness. This chapter both delineates general methodologies for the evaluation 
of compressor stabilization strategies and quantifies the performance of several ap- 
proaches which might be implemented in gas turbine engines. 

A compressor feedback stabilization system to prevent surge conceptually consists 
of sensors to detect fluid disturbances within the compression system, actuators to 
introduce desired perturbations, and a suitable control law connecting the two. The- 
oretically, only a single sensor and actuator are required, with many choices available 
for their type and location. Sensors may measure pressure, mass flow, velocity, or 
temperature within the compressor duct, upstream or downstream of the compressor, 
in the plenum, or at the throttle. Similarly, there are many methods to introduce 
unsteady fluid perturbations: varying throttle area, moving a plenum or duct wall, 
introducing or bleeding off mass flow, varying the heat addition in the plenum (when 
it is a combustor), as well as introducing a variable throttle between the compressor 
and plenum, to name a few. A representative list of sensing and actuation options is 

shown in Table 1. 

The central point of the chapter is to show that selection of sensor and actuator 
type and location is a critical factor in determining the effectiveness and practical- 
ity of a stabilization system. Because of the practical interest in this question, the 
chapter presents a methodology to compare different implementation alternatives, as 
well as carries out this comparison for a number of candidate strategies. Each sen- 
sor/actuator pair when coupled to the compressor, forms a different physical system, 
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Actuation 

Sensing 

Injection in compressor duct 

Inlet mass flow 

Valve close-coupled to compressor 

Plenum pressure 

Plenum bleed valve 

Compressor face Ptotal 

Plenum heat addition 

Compressor face Pstatic 

Variable plenum volume 

Plenum temperature 

Variable inlet geometry 


Fast inlet guide vanes 


Tangential inlet injection 


Auxiliary compressor stage 


Plenum mass injection 


Inlet duct bleed 



Table 6.1: Sensing and Actuation Options 

with differing dynamic behavior, physical limitations, and overall performance. The 
steps included in the evaluation of control strategies are thus as follows: 


• Models of specific systems are analyzed to elucidate the relative performance 
sensitivities to non-dimensional system parameters. 


• Numerical calculations are carried out to quantify the limits to control of various 
sensor/actuator pairs when connected by a simple proportional control law. 

• Optimal control theory is used to evaluate actuator effectiveness, given complete 
knowledge of the system state. 

It is noted that the methodology followed here was developed to obtain quantitative 

i 

comparisons in a manner which was thought to be most accessible to the to the tur- 
bomachinery community. The general problem of actuator /sensor selection has also 
been considered by a number of researchers for other applications (see, for exam- 
ple, Norris and Skelton, [33]; Schmitendorf, [42]; Muller and Weber, [31]) and these 
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approaches should certainly be considered in future investigations of this type. 

The goals of this work are to make clear the large influence that implementation 
(sensor/actuator selection) has on the effectiveness of a compressor stabilization con- 
trol scheme and to elucidate the basic physical limitations to control. We thus do 
not consider all possible actuators and sensors but chose a representative selection of 
schemes which are potentially applicable for a large variety of compression systems. 


6.2 Analytical Comparison 

We use the simplified lumped parameter model of compression system dynamics, 
developed in Chapter 3 and linearize it about a particular operating point. The be- 
haviors of systems with different sensor/actuator pairing are most clearly revealed 
from their transfer functions, which are defined as the ratio of the Laplace trans- 
formed system output (sensor signal) to input (actuator motion). Open loop transfer 
functions for three different actuators and four different sensors illustrated in Fig- 
ure 3-1 are presented in Table 6.2. The transfer functions for all the sensor/actuator 
pairs in Table 6.2 have the same denominator polynomial denoted as D(s) which is 
given by: 

£) W = j! + (i^r - B ^) 5 + ( 1 ”^)’ (61) 

but the numerator polynomials differ so that the inherent differences are captured by 
the latter. For generality, the expressions given in Table 6.2 use the equivalent throttle 
slope, m Te , of a throttle in parallel with a bleed valve, and equivalent compressor slope, 
m Ce of a compressor in series with a close-coupled valve. If no bleed valve is used, 
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Close- Coupled Valve 

Plenum Bleed Valve 
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Compressor 
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Compressor Face 
Static Pressure 
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w> 

W) 

D{s) 


Table 6.2: Open Loop Transfer Functions 


m T = m Tt . Similarly, if no close-coupled valve is used, m Ce = m c . In this section, 
we will use the simplest control law, a proportional relationship between input and 
output, to focus on the effects of actuator and sensor selection, but in a subsequent 
section, we will examine the impact of the form of the control law. 

6.2.1 Stability Modification Using Proportional Control 

System stability under the influence of a proportional control law is determined from 
the roots of the closed loop characteristic equation (this is derived in many texts; see, 
for example, Ogata, [34]) 

GoW + KG„U) = 0 < 6 - 2 ) 

In Equation 6.2, Gw(a) and GoW are the numerator and denominator polynomials, 
respectively, o{ the transfer function given in Table 6.2 and the gam, K, is a real 
constant of proportionality in the control law. The system will be stable if and only 
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if all the roots of the characteristic equation have strictly negative real parts. For the 
present case, where the characteristic equation is a second order polynomial, this is 
equivalent to requiring all the coefficients of the closed loop characteristic equation 

to be strictly positive. 

Equation 6.2 is the sum of two terms, G N (s) and Gd{- s)- When the control is off 
(K set to zero), stability is determined from the roots of C?o(s) on ^Y anc ^ thus, is 
the same whatever the sensor/actuator pair. As the gain K is increased from zero, 
system stability becomes increasingly modified by the G N (s) term and, because the 
various sensor/actuator pairs have different numerator polynomials, Gjv(s), the effect 
of feedback also varies. This is best illustrated by several specific examples. The first 
two examples show the effect of sensor type, whereas comparison between behavior 
in the second and third illustrate the impact of the actuator. 


Example 1: Close-Coupled Control Valve With Mass Flow Measurement 

Using the appropriate transfer function from Table 6.2, substituting into the charac- 
teristic equation, Equation 6.2, and rearranging results in the closed loop character- 
istic equation: 


s 2 + 


Bm 


— Bm c + KBv, 


T e 


: ) s+ ( 1_ 


m r 


+ 


m 


T* 


Kvc 

m Te 


= 0 


(6.3) 


For stability we require: 


( — Bm Ct + KBv c ] > 0 

\Bm Te ) 


(6.4) 
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and 


fi3 + £!;]>o (6.5) 

V ™T t m Tt J 

All the parameters in these two inequalities are positive numbers except for the com- 
pressor slope, m Ce . If m Ce < 0, as it typically is for high flow rates, the system will be 
stable with no feedback, that is with I< set equal to zero. As the flow rate becomes 
lower, m Ct will become less negative, reaching zero at the peak of the compressor 
characteristic and then moving to a large enough positive value, so that the system 
will be unstable without feedback. For sufficiently large values of the gain I<, how- 
ever, both the inequalities expressed in Equations 6.4 and 6.5 can be simultaneously 
satisfied and the system can always be stabilized. 

t 

Example 2: Close-Coupled Valve With Plenum Pressure Measurement 

In this case, as given in Table 6.2, the numerator polynomial of the transfer function 
contains only a constant term. The closed loop characteristic equation is. 

— + W )=0 ( 6 . 6 ) 

\Bm T<; ) \ m T t J 

The system can now only be stabilized if the equivalent compressor slope is small 
enough so that the term - Bm c ^ is positive; in other words, the equivalent 

compressor slope obeys the inequality m Ce < l/(B 2 m Te ). The ability of proportional 
feedback to stabilize this system is thus limited to a certain range of parameters. For 
many applications, B is unity or larger and the throttle slope, m Te , is on the order of 
ten to one hundred, so the useful range can be quite small. 
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Example 3: Plenum Bleed Valve With Measurement of Plenum Pressure 

From Table 6.2, the numerator polynomial for this case is 


G N (s) = -^{s~Bm Ct ) 


(6.7) 


and the closed loop characteristic equation is given by: 


s 2 + 


1 n Kv h 

— Bm Ce 5~ 

Bm r & 


l T c 


s + fl + Kvbm c ^\ — 0. 

V m T e / 


( 6 . 8 ) 


Whether the gain K is chosen to be positive or negative, it will have the desired 
effect on only one of the two coefficients of the characteristic equation. As a result, 
stabilization is limited to cases where m c . < 1 /B. The limitation is associated 
with the sign change between the leading and the constant coefficient of G N (s), as 
given by Equation 6.7, which implies that G„(s) has a zero in the right half of the 
complex plane. Systems whose transfer functions have numerators with zeros in the 
right half of the complex plane are called non-minimum phase systems. As discussed 
in Section 5.3.1, the ability to control non-minimum phase systems is known to be 
subject to certain fundamental limitations (see also Freudenberg and Loose, [15)) 
and this is just one manifestation of the generally poor behavior encountered in 
such systems. The remaining sensor/actnator pairs whose transfer functions appear 
in Table 6.2 have also been analyzed. For stability, all coefficients in the closed loop 
characteristic equation must be positive; from this requirement, the capability of each 
scheme to extend the flow range to high values of compressor slope can be determined. 
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Table 6.3: Limitations on Compressor Flow Range Increase With Proportional Con- 


trol 


The behavior of the different pairs is shown in Table 6.3, which summarizes the 
limitations of each pair. 

In some instances, for example, sensing compressor mass flow and actuating with a 
close-coupled valve, the range of parameters over which stabilization may be achieved 
is unlimited, although large values of gain may be required with a large compressor 
slope or large value of the B parameter. In these cases, the asymptotic behavior of 
the required gain is therefore given for large B and compressor slope, in order to 
show the trends to be expected in these regimes. Excessive gain must be avoided 
in practical situations as will be discussed below. As shown in Table 6.4, however, 
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the gain increases either linearly or quadratically with compressor slope and B and 
there will thus be a practical limitations on maximum slope or maximum B at which 
the system can be stabilized. An important point is that the limitations on nominal 
stabilization can, in most cases, be relaxed or removed by using a dynamic control 
law. Thus, the limitations expressed in this section reflect the combined properties of 
actuator, sensor, and control law taken together, and not necessarily of the individual 
elements. Within the restriction to a fixed control law maintained here, however, the 
comparison of different sensor/actuator pairs is both valid and useful. 


6.3 Practical Limits To Control 

The analytical results so far indicate that the ability to stabilize the system with 
proportional control depends strongly on proper pairing of actuator and sensor, as 
well as on the values of the system parameters, particularly the compressor slope 
m c and B. In the above examples, however, we have considered ideal linear systems 
in which only the nominal system dynamics system are considered. To address the 
issue of implementation, it is also necessary to consider bandwidth limitations and 
actuator constraints (for example, servo dynamics and mechanical stops) which are 
encountered in any physical realization. As a matter of definition, by actuator we 
refer here to the entire actuation system including the flow train element (e.g. the 
valve), the motor that drives it, and any included feedback elements. The bandwidth 
limitations may be imposed by the sensors, processor, actuator, or some combination 
of the three. System bandwidth must also be constrained to maintain stability in 


144 



the presence of unmodeled high frequency dynamics (see Section 5.3). Unless the 
bandwidth of the actuator is much greater than that of the compression system, 
there is a non-negligible time lag between the command output of the control law 
and the response of the flow train element. The lags introduced by the actuator 
generally result in reduced control effectiveness although, to some degree, they can 
be compensated for by use of a control law more sophisticated than proportional 
control. Another constraint on control effectiveness is the need to avoid actuator 
saturation 2 . For example, valve areas can only be modulated between 0 and 100% 
(i.e. the valve must be somewhere between full open and full closed). This implies 
that control gains must be kept small enough to avoid saturation for the expected 
range of perturbations. 

6.3.1 Sensor and Actuator Pairs and Fluid Model 

In this section, the limitations imposed by gain and bandwidth constraints are quanti- 
fied for specific cases. Five actuators and four sensors are studied as representative of 
a diverse set of implementation options. The selected actuators were: 1) injection in 
the compressor duct; 2) close-coupled control valve; 3) plenum bleed valve; 4) plenum 
heat addition; and 5) a movable plenum wall. The selected sensors were: 1) compres- 
sor duct mass flow; 2) plenum pressure; 3) compressor face static pressure; and 4) 
compressor face total pressure. These actuators and sensors are shown schematically 
in Figure 3-1. At the level of idealization used here, the close-coupled valve could 
be either at the compressor inlet or exit without changing the results. The linear 
2 An actuator is said to be saturated when its output reaches the limit of its allowable range. 
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lumped parameter models of the previous section were extended as required, because 
the various actuators imply additional system states. The differential equations that 
describe the system dynamics are obtained by performing balances on momentum in 
the ducts and mass and energy in the plenum as detailed in Chapter 3. The general 
description is nonlinear, but it is small perturbations which are of primary interest 
here, and the governing differential equations are linearized to yield sets of equations 

of the form: 

x = Ax + bu. (6-9) 

The nonlinear output equations are also linearized to obtain the resulting vector of 
output perturbations y defined as 


y = CX + du 


( 6 . 10 ) 


In Equations 6.9 and 6.10 , the linearized state variables x, the inputs tz, and the 
outputs y, are perturbations from the corresponding equilibrium values, and A, b, 
c, and d are appropriately dimensioned constant matrices. The state variables have 
been normalized as detailed in Table 6.4, so that unity magnitude for any of these 
perturbation variables has approximately the same physical significance. 

All twenty pairings of the five actuators and four sensors have been evaluated with 
a proportional control law. Such a comparison provides two useful results. One is the 
identification of sensor/actuator pairs which may be stabilized over a significant range 
of system parameters using the simplest possible control law. In addition, for those 
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Perturbation Variable 

Normalized By 

Flow 

Time mean flow through compressor 

Pressure 

Time mean compressor pressure rise 

[ Thermal input 

Mean compressor work 

Moving wall work input 

Mean compressor work 

Plenum bleed valve area 

r Area to fully close valve 

Close- coupled valve area 

Area to fully open or close 


(whichever is smaller) 


Table 6.4: Normalization Factors 


pairs with significant stabilization, the required gain gives a measure of the combined 
effectiveness of this choice of sensing and actuating locations. As was discussed, it is 
useful to eliminate those systems which could mathematically be stabilized but stand 
little chance of succeeding in an actual implementation. To this end, two constraints 
were imposed. First, the allowable magnitude of the normalized proportional gain was 
limited to be not more than twenty. For example, at the maximum allowable gain, a 
live percent change in compressor mass flow would yield a one hundred percent change 
in the allowable plenum bleed valve area; that is, the valve would be fully opened 
or closed. Second, the bandwidth of the feedback loop was limited by modeling a 
two-pole, low pass Butterworth filter in the feedback path. This filter can be given 
various physical interpretations such as probe dynamics, amplifier dynamics, actuator 
dynamics or unmodeled dynamics in the compression system itself. Whatever the 
interpretation, the insertion of the filter insures that the feedback path has finite 
bandwidth, a constraint which will always exist in practice. The study was carried 
out with the cutoff frequency of this filter maintained at ten times the Helmholtz 
frequency of the system formed by the plenum and compressor ducts and sensitivity 
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to this assumption will be examined subsequently. 

The figure of merit used to assess the sensor/ actuator pairs was to examine the 
stability boundaries in a compressor slope versus B parameter plane. Preliminary 
studies and consideration of the underlying system dynamics showed that these two 
parameters have a dominant effect on system stability. It is more relevant, for ex- 
ample, to quantify the amount of stabilization that can be achieved in terms of the 
compressor slope, which enters into the stability criterion in an explicit manner, rather 
than the change in mass flow at staff. Since the positive compressor slope typically 
increases as flow decreases, all else being equal, the greater the maximum stabilizable 
slope is, the greater the range extension will be. The relative extent of the stabilized 
region in this compressor slope B-parameter plane thus provides an appropriate and 
useful basis for comparison. The stability boundaries were computed by performing 
an incremental search over the three-dimensional (slope, B parameter gam) parame- 
ter space. For each fixed B, the value of gain which maximized the slope at instability 
(within the allowable range), as well as the corresponding slope, was found. The sta- 
bility boundaries in the B parameter versus slope plane thus represent the maximum 
slope which could be stabilized using any normalized gain with a magnitude less than 

twenty. 

6.3.2 Results of the Control Scheme Evaluations 

The results of the calculations are summarized in Figure 6-1 which shows the stability 
boundaries for the twenty sensor/actuator pairs. The figure is broken into four plots, 
one for each sensor. Within each plot, the five curves indicate the different actuators. 
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The region below and to the left of any given line is the region in which stabilization 
can be achieved. In the upper left hand plot, for example, all the region to the left 
of the dash-dot line represents the range of compressor slope and B in which the 
combination of compressor mass flow sensor and close-coupled valve is capable of 
suppressing the instability. 

Several general conclusions can be drawn from the results in Figure 6-1: 

1. The overall trend is that control becomes more difficult as the compressor slope 
and B parameter increase, with the maximum stable slope decreasing with 

increasing B. 

2. Only the actuators located in the compressor duct, which act upon the com- 
pressor duct momentum (injector and close-coupled control valve), are capable 
of stabilization at steep slopes over the full range of B. 

3. Plenum heat addition gives little or no stability enhancement 

4. In general, there is no best sensor independent of the actuator. 

For reference, the B parameter that might be associated with large axial gas turbine 
engines is very roughly 0.5 (e.g., Mani [27]). For other engine geometries and indus- 
trial compressors with large downstream volumes, B parameters of 2 or greater might 

be encountered. 

A more specific conclusion is given by the comparison of the results with the mass 
flow sensor to the other sensor locations. As B reaches a value of roughly unity, the 
ability of all the pairs to stabilize the system becomes quite small, except for the 
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close-coupled valve and the injector which use mass flow sensing. This points out 
clearly that not only is actuator position important, but sensor position is as well. 

The conclusion about the effect of actuator position is one that is in general accord 
with intuitive ideas of system behavior, but that having to do with sensors is generally 
less familiar. It is therefore worthwhile to give some physical motivation for the impact 
of sensor position. The different dynamics brought about by the various sensors can 
be understood with reference to the non-dimensional characteristic equation for the 
unsteady system behavior with no feedback: 

(6U) 

For stability the coefficients of the second and third terms must be positive, so that 

- ( 6 ‘ 12 ) 
Jj TTiTc. 

m Cc < rn Te ( 6 - 13 ) 

As described by Greitzer, [19], the mechanism of instability can either be static, 
corresponding to the inequality in Equation 6.13 being violated or dynamic, corre- 
sponding to the violation of that in Equation 6.12. Whichever of these events occurs 
the underlying cause is a positive slope of the compressor pressure rise charactenstic. 
Therefore, let us examine how the destabilizing effect of positive slope is ameliorated 
when different sensing schemes are used. As a case of practical interest, we consider 
the close-coupled valve with two different sensors, one measuring compressor mass 
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flow and the other measuring the total pressure at the compressor inlet face. 

In the first of these, the valve position and, hence, the valve pressure drop, is 
proportional to the sensed perturbations in compressor mass flow. The pressure 
perturbations across the compressor and across the valve are both proportional to 
the mass flow perturbations. Because the valve is just downstream of the compressor, 
the two act in series, creating an effective compressor slope which is the sum of the 
(positive) slope across the compressor and the (negative) slope across the valve. It is 

this combined characteristic which the system sees . 

Suppose the constant of proportionality between sensed mass flow and valve open 
area is K, and the rate of change of valve pressure drop per increment in open area 
is v c . For a given mass flow perturbation, the effective slope of the compressor will 

change from m c , , with no feedback, to K - fedback is Se " Sing 

the compressor mass flow and feeding the signal back to the valve actuator thus works 

directly on the cause of the instability, the positive compressor slope. Insertion of the 
slope (m c . - Kv,) in Equation 6.1 1 in fact gives Equation 6.3, the previously derived 

closed-loop characteristic equation for this situation. 

A different situation prevails for the inlet total pressure sensor. The inlet total 
pressure perturbation is related to the derivative of the inlet mass flow, i.e., the fluid 
acceleration in the compressor duct, through the unsteady Bernoulli Equation 

„ f 1 du J r (6.14) 

Po 0 - Poi = P J Q dx v 


151 


Non-dimensionalizing and linearizing gives 


-4>i = — d </>/ dr 


( 6 . 15 ) 


where fi is the total pressure perturbation at the compressor face in non-dimensional 
form. If the valve area perturbation is proportional to the compressor inlet total 
pressure, the result is to create a pressure change across the valve proportional to the 
acceleration, which is the derivative of the state variable cj>. In other words, the effect 
of feedback in this case is to alter the overall pressure difference from duct inlet to 
exit for a given fluid acceleration rate. As far as the system is concerned, this is seen 
as a change in the equivalent duct length, £, of the compressor. That is, the longer 
the duct, the larger will be the instantaneous pressure change required to produce 
a given fluid acceleration. An increase in effective length of the duct does several 
things to the system, the most important being that it changes the effective value of 
the B-parameter, which scales as 1 /y/Z. This can be seen in the resulting closed loop 
characteristic equation which is given by: 


(1 + Kv c )s* + f ^ Buie, + 2 Kvc4> c B + 


— I<v c 


Bm-i 


J V 


+ 


2 Kv c <f> c 

m T 


( 6 . 16 ) 

in which the critical value of B for instability is increased by y/1 + Kv c , with (1+Kv c ) 
being just the factor by which the effective length is increased. In the case of total 
pressure sensing, however, only the dynamic instability is influenced, because the 
increase in effective mass does nothing to alter the static stability. The instability 
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limit is thus still m Ce > ^r e - 

Somewhat similar arguments can be made for other sensor/actuator pairs; the 
main point is that the use of different sensors for the feedback gives the system quite 
different dynamical behavior. 

6.3.3 Effect of Control System Bandwidth 

The preceding study of the sensor/actuator pairs was carried out with the bandwidth 
of the feedback loop limited by a two-pole Butterworth inserted in the feedback path. 
The cutoff frequency, u c , of this filter was held fixed at ten times the compression 
system Helmholtz frequency u> h and we can also examine the influence of this param- 
eter, in other words, of controller bandwidth, on the stabilization process. Figure 6-2 
shows the changes in instability onset that occur with different controller bandwidths 
for the close-coupled control valve, with feedback on mass flow, and a B-parameter 

of 2. 

In the figure, the horizontal axis is the controller gain, and the vertical axis is the 
slope of the compressor characteristic at instability. Curves are shown for values of 
uJiOh from 1.0 to 100, representing extremes of this ratio, and it is seen that large 
changes in the maximum slope that can be stabilized result from these changes in 
bandwidth of the controller. 

Several trends are exhibited in Figure 6-2. First, for a controller bandwidth of 
u J Uh = 1, use of the control is actually destabilizing. The more gain the less steep 
the compressor slope that can be achieved. In this regime, the control system is 
adversely affected by the slow actuator dynamics. Second, as u c /io h becomes larger, 
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increasing the controller gain increases the compressor slope that can be attained. 
Third, for a given level of bandwidth, increasing gain increases stability only up to a 
point; beyond this point the stabilization decreases as gain increases. The maximum 
slope which can be stabilized, however, always increases with bandwidth. 

The necessity to go to higher bandwidth at higher slope arises because, as the 
slope increases, the system dynamics become faster relative to the undamped natural 
frequency (Helmholtz frequency). The relevant time scale becomes the fluid acceler- 
ation time in the compressor duct and thus the actuator must be fast relative to this 
time scale rather than that set by the Helmholtz frequency. This point is an impor- 
tant one since it is difficult, in many cases, to engineer actuators with bandwidths 
many times that of the Helmholtz frequency. The analysis suggests that actuator 
bandwidth will be a prime determinant of practical system performance and that 
research and development on this topic is useful to pursue. 


6.4 Comparison Based On Actuator Response To 
Disturbances 

From the study of individual sensor/actuator pairing using a proportional control, 
the close-coupled valve and injector emerge as the most promising actuators. A 
further question to address, however, is whether the type of compensation (control 
law) or choice of sensors would affect this conclusion. In this section, linear, optimal 
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stochastic, control theory is used to provide a definite answer to this question . 

6.4.1 Background 

To motivate the approach taken, some basic aspects of linear dynamic system stabi- 
lization are first presented. For small disturbances, the compression system, actuators 
and disturbances are described by the linear differential equation set: 


x = Ax + bn + Lw 


(6.17) 


where x is a vector of system states, u is a scalar representing a particular actuated 
variable, w is a p-dimensional vector of external disturbances, and A, b, and L are 
constant matrices of appropriate dimensions. This is just Equation 6.9 modified to 
include the effect of external disturbances. It is known (see, for example, Kwakernaak 
and Sivan, [24]) that systems of this type can be stabilized for all conditions using the 
control law u = -Ax, where K is now a one-by-five constant gain matrix, provided 


only that the pair {A,B} is controllable. Such a control law is known as full state 
feedback. For a system with n states, this would require n properly placed sensors. For 
the situation of interest here, over the range of parameters which have been analyzed, 


the requirement of controllability is met. Thus, with enough properly placed sensors, 
nominal stabilization of the idealized linear system is theoretically possible. (See also 

controllability. The degree of controllability explici y f stochastic approach explicitly 

r r „r r tzszz* .» - RM s m e» 

square) activity rather than fixed bounds. 
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Chapter 5 which considers the case of a single sensor.) 


6.4.2 Methodology 

All actuators are not equally effective. All else being equal, an ineffective actuator 
will require more motion to maintain stability in the presence of disturbances than 
an effective one. In the analysis, therefore, the actuators are compared based upon 
their minimal required RMS (root mean square) response to a persistent broadband 
disturbance, while maintaining system stability. This comparison is independent of 
choice of sensor, because it is assumed that the state of the system is known at all 
times. Further, the comparison is based upon the minimal possible RMS amplitude, 
and hence there is no question as to whether a particular actuator would perform bet- 
ter if another control law were used. In this sense, the comparison is also independent 

of the control law. 

Specific details of the computations performed are included in Appendix D and 
a more general treatment can be found in Kwakernaak and Sivan [24]. The central 
concept utilized is that, if the disturbance can be described as a stationary, Gaussian, 
white noise process, a particular gain matrix K can be found which will minimize 
the root mean square value of the actuated variable u. The gain matrix K and the 
RMS value of u will depend upon the matrix L which determines how the disturbance 
enters the system and upon the statistical properties of the disturbance. 

For this analysis, the disturbances driving the system must be chosen. Three 
disturbances were studied: 1) a compressor with an unsteady pressure rise; 2) an 
unsteady heat release in the plenum chamber; and 3) an unsteady outflow through the 
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throttle. In order to compare disturbances which have effectively the same strength, 
in each case, the input disturbance amplitude (white noise intensity) was normalized 
to that required to produce a one percent RMS fluctuation in compressor mass flow at 
a fixed, stable operating point of = 0.20. The results of this analysis are illustrated 
in Figure 6-3. In the figure, the maximum allowable slope at which the system can 
be stabilized is plotted as a function of B for live different actuators, based on the 
restriction that the RMS actuation is no greater than 25% of the full actuator range, 
for the most deleterious type of disturbance. The choice of this level of maximum 
RMS actuation is somewhat arbitrary, but it is taken to be one that can be ach.eved 
in practice. In addition, computations have been carried out at other levels, and the 
results show similar trends. 

For all the actuators, the maximum allowable slope that can be attained decreases 
as the B- parameter increases. However, the effect of varying B on the maximum 
allowable slope varies markedly between different actuators. As with proportional 
control, those actuators which are most closely coupled to the compressor ( injection 
and the close-coupled valve) are the most effective. The figure shows that, for B 
greater than unity, the maximum allowable compressor slope is quite limited. (For 
reference, the characteristic of the turbocharger used in the experiments, as described 
in Chapter 2, has a maximum slope of approximately six, i.e., m c - 6). In addition, 
except for the close coupled valve, whose performance becomes independent of B , the 
maximum allowable slope decreases monotonically with B parameter. Heat addition, 
which might seem attractive because of ease of implementation through fuel injection, 
shows little potential for stabilization at values of B larger than unity. 
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Another result is that the behavior of the plenum bleed, the injector, and the 
moving wall are roughly comparable, although the injector has some advantage for 
larger values of B (greater than two, say). At these high values of B parameter, 
however, the close coupled valve has a clear advantage, over all of the other schemes 
examined, in stabilizing the system. 

6.4.3 Effect of Disturbance Type 

The influence of the disturbance type on the minimal required actuator activity is 
shown in Figure 6-4 for the close coupled valve, the actuator shown to be most 
effective in enhancing flow range for large values of B. The horizontal axis is the 
slope of the compressor characteristic, and the vertical axis is the normalized RMS 
actuation level. Three curves are plotted, showing the RMS level corresponding to 
three disturbances of equivalent strength but different type. The difference between 
the three curves indicates the impact of disturbance type on the ability to control 
the system. Compressor pressure rise disturbances, such as might arise from local 
unsteady flow in the impeller or diffuser, create a situation that is more difficult 
to control than disturbances due to combustor heat release or throttle mass flow 
fluctuation. 

The implication of Figure 6-4 is that the nature of the disturbances driving the 
system is an important factor in setting the requirements for stabilization. There is 
little known at present about the detailed disturbance structure within turbomachin- 
ery and engines and characterization of these unsteady phenomena is thus a research 
item of considerable practical concern. 
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6.5 Summary 


An evaluation of strategies for the control of compression system surge has been 
carried out as a first step towards developing rational design procedures for surge 
stabilization. A basic result is that the close-coupled actuators and sensors which 
measure and act upon the momentum of the fluid in the compressor duct are the most 
effective for geometries and compressor slopes of interest for gas turbine applications. 
Although this result was qualitatively anticipated based upon the system physics, 
the analysis has quantified the severity of these trends, showing them to be extremely 
important over the parameter range of interest. The following specific conclusions 

can be made: 

• Proper choice of actuator and sensor is an important part of the overall design 
of a surge stabilization system. 

• Mass flow measurement with either a close-coupled valve or an injector for 
actuation are the most promising approaches of those evaluated. 

• Fuel modulation is not a promising candidate for practical ranges of system 
parameters. 

• Characterization of compression system disturbance sources is important for 
determining the requirements for stabilizing control schemes. 

• As either the compressor slope or the B parameter increase over the anticipated 
range of interest, control becomes substantially more difficult. 
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• Actuator bandwidth can be an important constraint in many practical imple- 
mentations. 
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Compression System B Parameter Compression System B Parameter 


Figure 6-1: Influence of sensor and actuator selection on maximum stabilized slope 
for a bandwidth and gain limited system 
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Slope at Stability Limit 



Figure 6-2: Influence of normalized control system bandwidth (p£) or 
gain required and slope at instability onset;close-coupled control valve 


the feedback 
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Compression System B Parameter 


Figure 6-3: Maximum slope at which the RMS actuation will be no greater than 25% 
of the maximum actuator authority 



Slope of Compressor Characteristic 

Figure 6-4: Influence of compression system disturbance type on actuator motion; 
close-coupled valve, B = 2 
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Chapter 7 


Close-Coupled Control 

The analysis of Chapter 6 indicates that close-coupled control schemes, i.e., measure- 
ments and actuation in close proximity to the compressor, appear to be the most 
promising generic strategy for this class of systems. Further investigation of this 
specific approach is therefore warranted. This chapter presents analytical and exper- 
imental results for a particular implementation in which a valve close-coupled to the 
compressor serves as the actuator and a hot wire anemometer in the compressor duct 
is used as the sensor. 

The chapter is organized into two main sections. The first considers some the- 
oretical aspects of close-coupled stabilization and shows that there are benefits to 
stabilizing the system using a feedback controlled valve as opposed to a simple flow 
restriction. In the second section, the results of close-coupled stabilization experi- 
ments are presented and compared with theory. 
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7.1 Close-coupled Control Theory 


7.1.1 Introduction 

As discussed in Chapter 6 the compression system will be linearly stable at any operat- 
ing point where the slope of the compressor pressure-rise characteristic (pressure- rise 
as a function of flow) is negative. The basic mechanism for close-coupled stabilization 
using a close-coupled valve is conceptually straightforward, the valve in series with the 
compressor effectively forms an equivalent compressor whose pressure-rise character- 
istic is is negative at the desired operating point so that the overall system is stable. 
Stabilization can be accomplished using either a valve whose area is modulated in 
response to measured flow perturbations or with a fixed flow restriction, i.e., with or 
without feedback, respectively. These two approaches are illustrated schematically 
in Figure 7-1, which shows, for both cases, the compressor characteristic, valve char- 
acteristic, and combined characteristic which is negatively sloped at the operating 
point. In the case of the feedback controlled valve, a family of dashed curves is shown 
representing the valve pressure-drop characteristics corresponding to a range of fixed 
valve areas. The feedback control adjusts the instantaneous valve area according to 
the measured flow rate resulting in the solid curve which is identified in the figure as 
the valve characteristic with feedback. 

7.1.2 Pressure-drop Penalty With and Without Feedback 

As suggested by Figure 7-1, the pressure loss across the valve required for stabilization 
can be reduced using feedback and this possibility will now be investigated. Consider 
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the basic closed-loop characteristic equation for a compression system with a close- 
coupled valve actuator and compressor duct mass flow measurement (see Chapter 6): 

s 2 + (1 l(Bm T ) - B(m c -m * e ) - BI<(s)v c )s + (l-(m c - m Vc )/m T + K(s)vJm T ) = 0 

(7.1) 

In Equation 7.1, B is the Greitzer B parameter, K(s) is the transfer function of the 
compensator (not necessarily the constant gain of Chapter 6), m c is the slope of 
the compressor pressure rise characteristic, m Vc and m T are the slopes of the close 
coupled valve and throttle pressure drop characteristics and u c is the valve sensitivity 
(change in non-dimensional pressure drop per unit change in valve closure). With the 
exception of m c all of these coefficients will always be positive. 

I 

Close-coupled Stabilization Without Feedback 

For stability all the coefficients of the characteristic equation, Equation 7.1, must 
be positive. This could be accomplished with no feedback at all (i.e., A (s) = 0) 
if m Vc > m c that is, if the slope of the pressure-drop versus flow characteristic of 
the close-coupled valve is steeper than the positively sloped pressure-rise versus flow 
characteristic of the compressor. However, such a stabilization mechanism has an 
associated cost because the slope of the valve pressure-drop characteristic is directly 
related to the pressure loss across the valve. For an incompressible flow: 

= C{a c )<f>v 2 , 
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where C{a c ) is a fixed constant for any particular valve opening which depends upon 
the particular geometry and piping arrangement for the valve. The slope and pressure 
drop across the valve are thus related by: 

m Vc = 2'& v /<f> v (^*3) 

which is independent of the valve coefficient C(a c ) and. cannot be modified by chang- 
ing the valve design. The use of a close-coupled valve not augmented by feedback 
control therefore has associated with it an unavoidable pressure loss penalty given by: 

> <A c m c /2 ( 7 - 4 ) 


which is obtained from Equation 7.3 and the stability requirement, m Vc > m c . 

Close-coupled Stabilization With Feedback 

Now consider the situation when the valve area is modulated in response to the 
measured flow using a proportional feedback. In this’case, from Equation 7.1, it is 
sufficient for stability that, 


m c — (m Vc + Kv c ) > 0, ( 7 -5) 

so that the valve slope, m Vc , required for stabilization may be reduced by increasing 
the proportional gain, I<. The pressure loss across the valve required for stabilization 
as given by Equation 7.3, is accordingly reduced. The use of feedback is thus seen 
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to reduce the required pressure drop required for stabilization. The sum, m Vc + I<v c 
appearing in Equation 7.5 can be interpreted as the slope of the effective pressure 
drop characteristic of the feedback controlled valve. Thus, as depicted schematically 
in Figure 7-1, the effective pressure-drop characteristic of the feedback controlled 
valve can have a steep slope at the desired operating point with very little steady 

state pressure-drop. 

Pressure Drop Penalty With Feedback 

Ideally, the required steady state pressure drop across the control valve could be made 
arbitrarily small by sufficiently increasing the gain. For several reasons, in practice, 
some drop across the valve will be required. First, in order for the valve to be 
able to produce pressure perturbations which act to accelerate as well as decelerate 
the duct flow it must be able to both open and close slightly from its equilibrium 
operating point. This requires some amount of valve closure and corresponding steady 
state pressure loss. The actual amount required depends upon both the proportional 
gain and on the level of disturbances which must be counteracted, and is therefore 
application dependent. 

The second factor which may necessitate some steady-state pressure-drop across 
the control valve is that there are inevitably additional dynamic elements in any 
real system. The effects of these dynamics are not accounted for in Equation 7.1, 
which indicates that any slope can be stabilized by making the gain sufficiently large. 
In reality, it is not possible to make the gain arbitrarily large without eventually 
experiencing instability. On the other hand, if the gain is too small the system 
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will also be unstable. As a result, there will generally be some maximum value of 
the equivalent compressor slope, m Cc = m c - m Vc which can be stabilized. This 
is discussed in Chapter 6 and such maximum values of equivalent slope are plotted 
in Figure 6-2 in which additional dynamics, represented by a low-pass filter, are 
included in the calculation. Let the maximum equivalent compressor slope which can 
be stabilized with the proportional control in the presence of the additional dynamics 
be denoted, m Ce . If m c > m Ce the system can still be stabilized but only by increasing 
m Vc so that m c - m Vc < rri Ce which in turn requires some steady-state pressure drop 

across the valve. 

7.1.3 Optimal scheduling 

As discussed above, at any given operating point, with or without feedback, there 
will be some minimum pressure-drop required for stabilization. If the compressor 
is to be stabilized over a range of flows it will, therefore, be desirable to schedule 
the mean valve opening with operating point so as to provide only the minimum 
pressure-drop required for stability. This will be referred to as the optimal scheduled 
valve area. Since the use of feedback reduces the valve pressure drop requirement, the 
optimal scheduled valve area will be different depending on whether or not feedback 
is used. In the following sections, numerical calculation of the optimal schedules will 
be considered both with and without feedback. 
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Optimal Schedule Without Feedback 

We now obtain an expression for the minimal valve pressure drop required to stabilize 
the system with no feedback. With no feedback, the stability requirements are, from 

Equation 7.1: 


1 

Bm T 


— Bm Ct 


> 0 


m Cc 

TYl t 


> 0 , 


(7.6) 

(7.7) 


where, by definition, m c . = m c - m v «. We assume that it is Equation 7.6 which is 
the more restrictive. Equation 7.6 can be rearranged to obtain the maximum stable 

value of m Cc as, 


For an incompressible flow through the throttle m T = 2fa/fa, where p is the plenum 
pressure and fa is the flow through the throttle. At equilibrium fa = ~ where 

^ c , and i/> v are the compressor exit pressure and valve pressure drop respectively and 
<j) c = fa. Making these substitutions Equation 7.8 can be written as: 


m Ce 


< 


<t>c 

2B 2 ('$ c - fa) 


(7.9) 


From Equation 7.9, the maximum stable slope with no pressure-drop across the valve 
is then given by 


m c < 


fa 

2 B 2 V C ' 


( 7 . 10 ) 
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If the compressor slope is small enough to satisfy Equation 7.10 the optimal valve 
area is therefore fully open. For steeper compressor slope some loss across the valve is 
required. Assuming incompressible flow through the valve so that m Vc = 2 ip v / <{> c , and 
using Equation 7.9 we find the stability that at the stability boundary, m Vc satisfies 

the quadratic equation: 

m v 2 - (m c + 2'$ c /phic)m Vc + (2^ c m c /<f> c - 1/B 2 ) = 0. (7- 11 ) 

We are interested in cases where some pressure drop across the valve is required 
from stabilization, so from Equation 7.10, 2'H c m c /<j) c - l/B 2 > 0. Using this fact it 
is relatively straightforward to show that Equation 7.11 has two real, positive roots. 
Since we are interested in the minimal stabilizing pressure drop, we choose the smaller 
of the two roots, which we will denote by m Vc . The minimum pressure drop is then 

calculated using, 

t/U = rn Vc (f> c /2. (7.12) 

Using this value for the valve pressure drop, the required valve area can then be 
calculated, with the specific relationship to be used 4ependent upon the particular 
valve geometry. 

Optimal Scheduling With Feedback 

It will not generally be possible to obtain a closed form solution to determine optimal 
scheduling with feedback and so a numerical procedure must be implemented. The 
following approach was taken here to determine the minimal stabilizing pressure drop 
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and corresponding valve area at a given operating point. 

1. Using a polynomial curve fit for the compressor characteristic, calculate the 
compressor pressure rise and compressor slope at the desired flow coefficient. 

2. Enter a loop in which the valve pressure-drop is incremented starting with zero 
pressure drop. Inside this loop the following steps are taken. 

(a) The plenum pressure and resulting slope of the throttle characteristic are 
determined. 

(b) The linear system, A,b,c matrices describing the compression system dy- 
namics are calculated. 

(c) The compression system model is augmented with a linear system model 
for the actuator dynamics. 

(d) The resulting system is tested to determine if there is any stabilizing gain. 
The algorithm described by [28] was implemented as a Matlab routine for 

this purpose. 

(e) If there is no stabilizing gain the valve pressure drop is incremented and 
the loop is repeated. 


7.1.4 An Illustrative Example 

preceding theory will now be illustrated with an example in which the optimal 


The 


scheduled valve areas are calculated with and without feedback using the compressor 
characteristic for the Holset Model HID turbocharger described in Chapter 2. The 
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positively sloped portion of the characteristic was experimentally measured by using 
a close-coupled valve without feedback to stabilize the system. 

For the example, a system B parameter of 2 was used and the actuator dynamics 
were represented by a second order, low-pass Butterworth filter, as in Chapter 6. The 
-3dB cutoff frequency of this filter was set to ten times the Helmholtz frequency of the 
compression system. The compressor tip Mach Number was 0.76. The flow through 
the close-coupled valve was assumed to be incompressible and the total pressure loss 
through the valve, which dumped directly into the plenum, was taken as 1.98 times 
the dynamic pressure based upon the minimal valve area. This value was based upon 
empirical data obtained for the valve used in the experimental facility. 

Optimal Scheduled Characteristics 

The resulting optimal scheduled performance with and without feedback for this 
illustrative example are shown in Figure 7-2. The figure shows net pressure-rise for 
the compressor /valve combination as a function of flow coefficient for three cases: 1) 
Following the optimal valve area schedule without feedback; 2) following the optimal 
valve area schedule with feedback; and 3) with the valve fully open (compressor alone). 

The vertical distance from the compressor alone characteristic to either of the 
optimal characteristics represents the minimal required valve pressure drop for stabi- 
lization. The figure shows the basic benefit of using feedback; the requisite pressure 
drop is decreased and the resulting delivery pressure of the stabilized system is there- 
fore higher. For this particular case it can be seen that, with feedback, essentially 
zero pressure-drop across the valve is required for flow coefficients from approximately 
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0.135 ( the natural surge point) down to 0.08. In this range the compressor slope is 
less than the maximum stabilizable slope and all the stabilization is accomplished 
using feedback. (Some additional pressure drop would be required depending on the 
disturbance level as discussed in Section 7.1.2.) At lower flow coefficients the com- 
pressor slope increases and some pressure-drop must be taken across the feedback 
controlled valve in order to stabilize the system. The corresponding optimal valve 
area and gain as a function of flow for the feedback controlled valve are shown in 
Figure 7-3 and Figure 7-4 respectively. These figures show that the optimal valve 
area and gain are highly nonlinear functions of flow, and thus a purely linear control 
scheme cannot be used if the optimal performance is to be obtained. The system 
implementation could, however, be split into a linear and nonlinear portion. A linear 
proportional control would act on a high pass filtered flow measurement signal (ac 
coupled) which provide an indication of the rapid Bow perturbations, and a nonlinear 
control would act on a low-pass measurement of the flow to slowly adjust valve area 
and gain to follow the optimal schedules shown in Figure 7-3 and 7-4. It is noted 
that a large part of the nonlinearity in the optimal gain is due to the particular valve 
geometry used here in which the valve area varied linearly with valve angular posi- 
tion. This valve geometry results in a valve whose per unit change in pressure drop, 
per unit change in angular position, varies significantly with the nominal valve area. 
As a result, the proportional gain of the controller must vary to maintain a constant 
overall loop gain. By properly tailoring the valve geometry it may be possible to 
maintain a constant gain, thereby simplifying the implementation. 


174 


Effect of Flow Rate 


It is interesting to examine the effect of changes in mean flow rate on the optimal 
scheduled performance. To do this, the optimal scheduled performance was recalcu- 
lated using a new compressor characteristic obtained by shifting the previous com- 
pressor characteristic so that corresponding points would occur at a higher flow. That 

I 

is, the new compressor characteristic, c \4>c ) is gi ven by : 

= ( 7 ‘ 13 ) 

where is the original compressor characteristic and A <j) c is the shift. The optimal 
characteristics for the original characteristic and the shifted characteristic are com- 
pared in Figure 7-5, for A<f> c = 0.15. Comparing the requisite valve pressure drop 
with and without feedback for the low and high flow cases, it can be seen that, all 
else being equal, the benefits of using feedback are greater for a compressor whose 
characteristic has steep positive slope at higher flow coefficients. This can be under- 
stood by considering that for the same compressor slope, the pressure drop required 
to stabilize the system without feedback increases directly with flow coefficient, as 
shown by Equation 7.4, while over a substantial range of flows a negligible amount of 
pressure drop is required for the feedback controlled valve. 
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7.2 Close-coupled Stabilization Experiments 

7.2.1 Description of Tests 

A discrete time (digital) proportional control was implemented which commanded the 
close-coupled valve area to be directly proportional to the mass flow at the compressor 
duct inlet (mass flow was derived from a constant temperature anemometer sensor). 
The control implementation is described in Chapter 2. The system was tested over 
a range of values of proportional gain in order to determine the optimal value to 
minimize the control valve pressure losses. The tests were conducted holding the tip 
Mach number, M T (or equivalently corrected speed) constant. A value of M T = 0.76 
was used for all tests. 

The same basic testing procedure was used for each gain setting. The operating 
point was first set by adjusting the area of the throttle valve located at the plenum 
exit, throttle The open area of the close-coupled valve was then gradually adjusted in 
order to stabilize the compressor with the minimum amount of pressure loss across the 
close coupled valve. After each small change in close-coupled valve area the turbine 
supply pressure would be adjusted to keep the compressor speed constant. Once the 
maximum stable valve opening was established, the time mean operating point data 
would be measured. The plenum exit throttle would then be closed slightly and the 
process would be repeated. This procedure was continued (that is operating points 
with gradually decreasing flow were measured) until the system could no longer be 
stabilized for any close-coupled valve area. With feedback, the system can not always 
be stabilized simply by reducing the area of the close-coupled valve. As the valve 
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is closed beyond a certain point, a high frequency instability develops, and when it 
is opened too much the system also goes unstable. At the limiting flow there is no 
longer any range between these two settings, that is any further closure or opening 
of the valve results in instability. 

7.2.2 Experimental Results 

The results of the close-coupled experiments are shown in Figure 7-6 which plots 
the total to static pressure rise characteristic of the compressor in combination with 
the close coupled valve for a range of gain settings.' Data points with the same 
gain are connected by straight lines. As described in the test procedure, each data 
point represents the maximum stable output pressure which could be obtained for 
a given gain and flow coefficient. For reference, the total to total characteristic of 
the compressor alone is also plotted in Figure 7-6. We are primarily interested in 
the optimum scheduled performance with and without feedback (see Section 7.1.2) 
as indicated in Figure 7-6 by the solid curves which have been fared in to form upper 
and lower bounds respectively for the experimental data. The figure shows that, as 
expected from the theory, the pressure-drop across the close-coupled valve required 
for stabilization can be decreased by using feedback,. thus providing higher output 
pressure over the stabilized range of flows. 
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7.2.3 Comparison with Theory 

It is of interest to compare the experimental results directly with theory. In order 
to make this comparison, an analytical model representing the actuator dynamics for 
the experiment is required. For these purposes, we include all the dynamics from the 
analog input to the computer to the actual position of the valve in one overall transfer 
function representing the actuator dynamics. The frequency response of this transfer 
function was measured with an HP3582A spectrum analyzer using a broadband drive 
signal. A fifth order, (five pole, four zero) transfer function was then fit to this data. 
This was found to be the lowest order model which could adequately represent the 
measured response. Both the experimentally measured frequency response and the 
response of the transfer function used to fit this data are shown in Figure 7-7, which 
shows that the actuator system has a natural frequency, as indicated by the peak in 
the magnitude characteristic, at approximately ten times the Helmholtz frequency. 
Using this model for the actuator dynamics and the model for the compression system, 
the theoretical optimal scheduled performance could be determined with and without 
feedback as described in Section 7.1.2. The compression system model was based 
upon the known geometry of the system and a polynomial fit to the measured steady 
state compressor characteristic. The quasi-steady compressor model (see 3 was used 

here. 

The resulting theoretical optimal scheduled performance and experimental data 
are shown in Figure 7-8. The curves indicating the theoretical optimal scheduled 
performance with and without feedback should bound the envelope of the experimen- 
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tal data. This is the case, indicating good agreement with theory, from the natural 
surge point at <f> c = 0.135 down to roughly <j) c = 0.085. Below this flow theory and 
experiment deviate both with and without feedback. 

The fact that the experimental data departs from the theory even without feed- 
back indicates that the discrepancy cannot be solely due to inaccurate modeling of the 
actuator or feedback system and we therefore turn attention to the open loop com- 
pression system model. It is thought that this discrepancy arises primarily as a result 
of unsteady compressible (organ pipe) behavior in the compressor duct, and in par- 

/ I 

ticular, the 0.3(m) section of duct between the compressor exit and the close-coupled 
valve. The model assumes that the flow in the compressor duct is one dimensional, 
and incompressible so that the compressor along with its upstream and downstream 
ducts may be represented using a lumped element element model. For this assump- 
tion to be valid, the acoustic time scale, that is the amount of time required for an 
acoustic wave to propagate down the length of the duct, should be small with respect 
to the growth rate of the instability. We will call the ratio of these two time scales, the 
reduced frequency; for model validity the reduced frequency should be much less than 
unity. The reduced frequency was calculated as a function of flow coefficient for both 
the inlet duct, (bellmouth inlet to compressor face) and exit duct (compressor exit 
flange to close-coupled valve) and the results are plotted in Figure 7-9. For both ducts 
the reduced frequency increases as the flow rate decreases. This is to be expected 
because at lower flows, the positive slope of the compressor characteristic becomes 
steeper and we would expect that the growth rate of the instability would increase 
and with it the reduced frequency. The reduced frequencies for the upstream duct are 
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higher because the upstream duct is longer. There is a prominent knee in both the 
curves (discontinuity in the slope), at approximately <j> c = 0.85, which corresponds 
to the flow at which the theory and experiment begin to deviate, and the reduced 
frequencies increase rapidly below this flow coefficient. The knee occurs as a result of 
a change in the nature of the instability which transitions from an undamped sinu- 
soid (complex conjugate poles in the right half complex plane) to a pure exponential 
growth (real poles in the right half plane) at <f> c = 0.085. as shown in Figure 7-10 
in which the real and imaginary part of the open loop compression system poles are 
plotted as a function of flow coefficient. The basis of the close-coupled theory is that 
there is no mass storage between the compressor and the close-coupled valve so that 
the characteristics of these components can be added together. Thus it is the valid- 
ity of the lumped assumption for the exit duct which is of primary concern, which in 
turn requires that the reduced frequency be much less than unity. For flow coefficients 
greater than 0.085, where agreement with theory was good, the figure shows that the 
reduced frequency based upon the downstream duct is less than 0.1., below this flow 
the reduced frequency increases rapidly as the system instability transitions to pure 
a pure exponential growth. It thus appears likely that the deviation between theory 
and experiment at lower flow coefficients is in large part due to unsteady compressible 
effects in the downstream duct. 

Because the flow sensor was located at the bellmouth inlet, while the close coupled 
valve was located downstream of the compressor, the effects of mass storage in the 
inlet duct also become quite important when the feedback control is used. Thus 
the fact that the reduced frequency of the upstream duct which approaches unity 
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at (f> c = 0.085 suggest that the performance of the system with feedback would be 
improved by collocating the actuator and sensor immediately at the compressor exit. 


7.3 Summary 

An analytical and experimental investigation has been made into the use of a partic- 
ular stabilization scheme in which measurements of compressor mass flow are fedback 
to adjust the area of a valve close-coupled to the compressor exit. The results are 
summarized as follows: 

• It has been shown that theoretically a compression system can be stabilized 
using a close-coupled valve either with or without feedback. The pressure-drop 
across the control valve required for stabilization can be reduced, however, by 
using feedback, thus providing higher delivery pressure over the stabilized range. 

• Ideally, with feedback, no pressure-drop is required across the close-coupled 
valve. Practically, some pressure-drop is required due to finite disturbances, 
and limitations on controller gain, and bandwidth. 

• A feedback controlled close-coupled valve has been implemented on a radially 
bladed, centrifugal compressor and it has been experimentally demonstrated 
that system can be stabilized using less pressure-drop with feedback over a 
wide range of flows. 

• Good agreement with theory was found when the theoretical growth rate of 
instability was consistent with the lumped parameter modeling assumptions for 
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compressor duct. 



pressure 


Fixed Valve 



Feedback Controled Valve 



Figure 7-1: Close-coupled stabilization with fixed valve and with feedback controlled 
valve (schematic) 
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Figure 7-2: Optimal scheduled performance with and without feedback for illustrative 
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Figure 7-3: Optimal scheduled valve area as a function of flow for illustrative example 
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Figure 7-5: Effect of flow rate on optimal scheduled performance similar compressor 
characteristics 
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Figure 7-6: Close-coupled stabilization experimental results 
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Figure 7-7: Bode plot of experimentally measured transfer function representing ac 
tuator dynamics showing fit to fifth order model 
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Figure 7-8: Comparison of close-coupled experimental results with theory 
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Figure 7-9: Reduced frequency as a function of flow rate for ducts in experimental 
apparatus 
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Figure 7-10: Open loop compression system pole location as a function of flow coef- 





Chapter 8 


Summary and Conclusions 


An in depth study of the use of feedback stabilization to prevent surge in compression 
systems has been carried out. While previous compressor stabilization studies have 
focused on particular implementations and experimental proof of concept, the work 
reported here has been directed towards obtaining a fundamental understanding of 
the problem. In this chapter, the approach taken is first briefly summarized and then 
specific conclusions are provided. 


8.1 Summary 

The sinusoidal frequency response of a radial bladed centrifugal compressor was exper- 
imentally measured over a wide range of equilibrium operating points. The resulting 
transfer functions were analyzed and used to validate a simple, lumped element dy- 
namic model of the compressor. This model includes the dynamic effects of fluid 
inertia but assumes that the pressure rise and flow through the compressor are alge- 
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braically related by the steady-state characteristic of the compressor. This is referred 
to as the quasi-steady assumption. 

Improvements in the model obtained by relaxing the quasi-steady assumption 
were also examined. In particular it was assumed that flow development process 
can be accounted for by a simple first order time lag between the instantaneous and 

steady-state pressure losses occurring in the flow passages. 

Much indirect support for the use of these types of simple component models for 
the compressor has previously been provided by the adequacy of overall compression 
system dynamic models which incorporated them. However, there has been very 
little unambiguous data and comparison with theory previously reported. The exper- 
imental results which have now been obtained show that m fact such simple models 
provide a good representation of the compressor; the element which is fundamentally 

responsible for the overall compression system instability. 

A second aspect of compression system modeling addressed in this thesis has 
been the development of overall system dynamic models which included a variety 
actuators and sensors of practical interest which had not been previously modeled. 
Development of these models in a form readily suitable for the analysis and design of 
control systems was of immediate utility in this thesis work for evaluating alternative 
control strategies and will also be beneficial to future researchers. 

The concept of stabilizing a compression system to prevent surge is a new one and 
its ultimate practicality remains to be determined. As an important first step towards 
providing this assessment, some fundamental limitations to compression system sta- 
bilization were identified and methods were developed for quantifying their severity. 
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This problem was approached by determining the relevant and useful control theory 
and then appropriately adapting and applying it to obtain results for this class of 
systems. This is thought to represent the first such examination of the fundamental 

limitation to feedback stabilization of compression systems. 

The study of fundamental limits to stabilization showed that the severity of such 
limitations, over the expected range of applications, could be greatly mitigated by 
proper actuator and sensor selection. The importance of proper actuator and sensor 
selection for this problem was further emphasized by the results of performing a first 
time, quantitative, side by side, comparison to assess the relative merits of a variety 
of potentially attractive actuation and sensing options. 

The overall sensor actuator evaluation study and consideration of the basic system 
dynamics indicated that close-coupled sensing and actuation has the best potential for 
being an effective, generic, approach to compression system stabilization. A detailed 
analysis and experimental investigation into this approach was, therefore, carried 
out. Basic design trade-offs were identified and quantified parametrically and the 
underlying mechanism for close-coupled stabilization was presented. 

A complete digital control and data acquisition system was designed and imple- 
mented on a representative compression system, in order to experimentally evaluate 
the close-coupled control concept. An 1-80386 microprocessor based system, capa- 
ble of implementing a linear, fourth order discrete time control algorithm at a lKHz 

sampling rate, was developed for this purpose. 

Using this system, experiments were conducted demonstrating the use of a close- 
coupled feedback controlled valve to extend the stable operating range of the com- 
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pression system. 


8.2 Conclusions 

On the basis of the analysis and experiments presented in this thesis the following 
specific conclusions have been reached: 

• Experimental measurements show that the centrifugal compressor which was 
tested has a clearly defined sinusoidal transfer function which relates the re- 
sponse of the flow through the compressor to a differential pressure fluctuation 
applied across the compressor duct (ie., an admittance). 

- Comparison of model predictions with the experimentally determined fre- 
quency response shows that for control and design analysis purposes the 
compressor can be assumed to be operating quasi-steadily, with all dy- 
namic effects due to the fluid inertia in the compressor flow passages. 

- The experimentally identified parameter values associated with the quasi- 
steady model are physically reasonable, indicating that both the model 
structure and assumed physics are consistent. 

- A simple unsteady loss model (the instantaneous pressure losses occurring 
in the compressor are related to the steady state losses by a first order lag) 
is useful for explaining observed departures from quasi-steady behavior 
which occur as the reduced frequency (based on flow through time) be- 
comes non-negligible. In particular, the model provides a physically plau- 
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sible means of accounting for the apparent increase in equivalent length (a 
purely geometric parameter) with decrease .in flow coefficient observed in 
a prior study [1]. 

- The experimental data shows that, as reduced frequencies approach unity, 
there is a substantial dependence of the real part of the compressor impedance 
(inverse admittance) with frequency. In particular, for a compressor oper- 
ating on a positively sloped portion of its characteristic (a condition where 
open loop instability can occur) the real part goes from being negative to 
positive as the frequency increases. With the sign conventions assumed 
here, this shows that the system goes from being active at low frequencies 
to passive at high frequencies. 

— The observed high frequency passivity, which is not predicted by the quasi- 
steady theory, is consistent with but not entirely accounted for by the 
simple unsteady loss model. 

- The fact that the compressor becomes passive at high frequency indicates 
that, above some frequency, it would no longer act as a destabilizing el- 
ement. Thus, this finding has the important implication that the quasi- 
steady model may be conservative in predicting control system bandwidth 

requirements. 

• There are many actuator/sensor pairs available for which the basic compres- 
sion system model is controllable and observable and the ability to nominally 
stabilize the system is generally not problematic. As a consequence, studies 
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which attempt to evaluate actuator sensor selection based upon their ability 
to nominally stabilize the system will provide little discrimination between the 
alternatives and are not recommended. 

• For many applications of interest, bounded actuation and the performance 
penalties required in order to robustly stabilize the system against unavoid- 
able high frequency modeling errors can be expected to severely restrict the 
ability to stabilize the system. 

• The limitations imposed by bounded actuation and high frequency modeling 
errors are much more severe for some actuator/sensor combinations than oth- 
ers, so that proper choice of actuators and sensors will be essential for many 
applications. 

• The non-dimensional slope, m c , of the compressor characteristic and the value 
of non-dimensional B parameter have been shown to be key parameters dic- 
tating the fundamental difficulty of the stabilization problem and the range of 
options available. The compressor slope determines if the compressor is active 
or passive. If it is active (positively sloped), its ability to supply the unsteady 
power required to sustain an instability scales with the magnitude of the slope. 
The B parameter determines the amount of coupling between the compressor 
and the downstream components. Thus, as values of either of these parameters 
increase stabilization becomes more difficult. 
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• For many potential applications, where B and compressor slope can reasonably 
be expected to be relatively large, (B > 1, m c > 2) only close-coupled control 
strategies which measure and act directly on th& fluid in the compressor duct 


will be viable. 

• Characterizing both the source and nature of the process disturbances acting 
on compression systems is important in determining the physical hardware re- 
quirements and expected range of stabilization of practical compression system 
stabilization schemes. , 

• The ability to stabilize a compression system using a servo controlled valve 
close-coupled to the compressor exit, with proportional feedback of compressor 
mass flow, has been experimentally demonstrated. The experiments verified 
theoretical predictions that the use of a high bandwidth feedback to modulate 
the valve allows the system to be stabilized with less steady-state pressure loss 
across the valve than would be required with a fixed valve (or one which is only 
slowly varied with operating point). 

8.3 Suggestions for Future Research 

There are a number of interesting questions and problems which were beyond the 
scope of this current research but should be be considered. These suggested areas 
of research are organized into two groups: those related to the specific problem of 
compression system stabilization and those which are more general problems in control 

theory. 
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8.3.1 Stabilization of Compression Systems 

• With the results and findings reported here providing a starting point, a basic 
systems engineering study to determine the potential risks and benefits asso- 
ciated with the stabilization of compression systems for particular application 
areas would now be quite useful. The results of such a study could help both in 
assessing the extent to which this technology should be pursued and to target 
the applications where it has the greatest potential benefit. 

• Either in conjunction With such a systems engineering study or alone, a demon- 
stration of feedback stabilization on an actual jet engine, industrial process 
compressor or gas turbine would answer many questions as to the practical 
applicability of this technology. 

• The work presented in this thesis strongly indicates that the availability of 
actuators with sufficiently high bandwidth and control authority for practical 
applications may be a potentially limiting factor. Development of such high 
performance actuators, either in conjunction with a particular application or as 
an independent effort is highly recommended. 

• A study to characterize both the source and nature of disturbances acting on 
typical compression systems, which has been shown here to be an important 
design factor, would provide extremely useful data for any of the studies rec- 
ommended above. 
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8.3.2 Control Theory 


The investigation of this particular control application has pointed out the following 
areas where it would appear beneficial to develop some general theoretical results. 

• There appears to be a great need for the development of control design and 
analysis techniques which are specifically suited to the problem of stabilizing 
open loop unstable systems. It appears that the majority of the controls lit- 
erature is oriented towards improving the performance of stable systems and 
either ignores or only indirectly addresses some of the difficulties inherent in 
stabilizing open loop unstable systems. 

Two problems which would be of particular interest to address are: 1) How to 
design a stabilizing control within the limitations of fixed hardware constraints, 
and 2) The related problem, given the objective of stabilizing a particular plant, 
what are the minimum hardware requirements in terms of experimentally veri- 
fiable performance specifications. 

• As discussed, actuator/sensor selection and evaluation problem has been inves- 
tigated by a number of authors. An interesting aspect of this problem which is 
recommended for future investigation is the one encountered in this thesis where 
the basic plant is not a fundamentally fixed entity with a pre-defined actuator 
sensor set to select amongst but instead, can be modified by the addition of 
arbitrary sensors and actuators. Under these conditions, can a physically real- 
izable choice of actuator and sensor which is in some sense optimal be directly 
synthesized? The problem is further complicated by the fact that addition of 
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certain actuators actually changes the plants dimension (ie, number of states 
in its state space description). For example adding an injector in the inlet duct 
of the compressor requires an additional state because the flow in the ducts 
upstream and downstream of the injector no longer need be the same. 
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Appendix A 


Modeling the Injector 


The injector is comprised of a short section of duct in which additional mass flow is 
injected coaxially as shown in Figure A-l. All flow stations referred to in the following 
analysis are as shown in Figure A-l. 

A.l Assumptions 

The following assumptions are made: 

1 . The densities of all fluid streams have a common density, />, which remains 
constant with time. 

2. Wall shear stress is negligible. 

3. The flow from plane- 1 to plane-4 is inviscid and irrotational. 

4. The static pressure across plane-4, P 4 , is spatially uniform. 

5. The static pressure, p 2 , and velocity, U2 across plane-2 are spatially uniform. 
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6. There is no loss in total pressure between plane-3 and plane-4. 

7. At plane-4 the flow velocity has two distinct velocities: the jet velocity, Uj, at 
the exit of the injector tube, and u m for the coaxial main flow surrounding the 
jet. 

8. The injector behaves quasi-steadily ie, we can consider the injector to follow its 
steady-state performance characteristics. 

t 

Let po, denote the total pressure at plane i. Define, 

Ap 0l = Po, - Po 3 
Apo 2 = Po 2 ~Po3- 

The ejector will be inserted between two sections of duct each having fluid inertia, 
thus, the upstream and downstream flow rates, rhi and m 2 will be state variables 
whose value are known at all time. We, therefore, want to express Apo, and Apo 2 as 
functions of nj[, m 2 , fluid properties, and the injector geometry. 

We first obtain the desired expression for Apo 2 as follows. At steady state, a 
momentum balance from plane-4 to plane-2 neglecting wall shear and assuming fully 
mixed conditions at plane-2 gives: 

p 4 A + m 3 Uj + rhiu m = p 2 A + m 2 u 2 , (A.l) 

where A is the cross sectional area of the duct. 
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Using the assumptions of uniform static pressure at plane-4 and no total pressure 
losses from plane-3 to plane-4 


Pa = Po 3 - 1 / 2 /> u j> 


(A.2) 


and by definition: 

P2 =P0 2 -1/2^2- ( A - 3 ) 

Making these substitutions gives: 


Ap 0 j = -1/2 pu) + rn^Uj/A + miU m /A + 1/2 pu\ - rn^jA. (A.4) 


Conservation of mass (continuity) gives: 


m 3 

= m 2 — mi 


m 2 — mi 


II 

^ , 

^rn 

mi 

1 

■SS 

1 

> 


ro 2 

«2 



where, Aj is the cross sectional area of the jet at plane-4. 
Substituting into Equation A.4 gives: 


Apo 2 = 


pA 2 


1 (m 2 - mi) 2 , (rh 2 -m^ 2 ^ m i 2 

'2 ( Ail A ) 2 {Ail A) (1 - {AjlA)) 


-m 2 


(A.5) 
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We now obtain the desired expression for Apoi as follows. With the assumptions 
made, we can apply Bernoulli’s equation from plane-1 to plane-4 to obtain. 


P4 = Po, - l/ 2 P U m- 


(A.6) 


Using Equation A.2, and continuity gives, after some manipulation, the desired ex- 


pression: 


mi 


A 1 ( (m 2 -m i) 2 - 

Ap °'~ 2pA 2 l {A JAY + (1 -HA J IA) + (A 1 IAfl 


(A.7) 


Equation A.7 and Equation A.5 are then made non-dimensional by normalizing 
pressure with respect to \/2pu T 2 , and mass flow with respect to pu T A c , where, u T 
is the compressor tip speed, and A c is the annular inlet area of the compressor. 
Defining, j3 = Aj/A, A* = A/A c , fa = mi/pu T A c , <£,„ = Ap 0 J(l/2pu T 2 ), and, 
= Apo 2 /(l/ 2 / 9 Uj 2 ) the desired non-dimensional expressions are obtained: 


* _ 1 ( ( fr -&) 2 , * a 

i \ u / . \ n I ’ 


(A*y 


(i -PY)' 


(A.8) 


and 


vlr = 

(A*) 2 


(fa - fa ) 2 (fa - fa ) 2 

■ w " + i 


h* fa* 
'•(!-/?) 2 


(A-9) 
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Figure A-l: Injector showing flow stations for analysis 
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Appendix B 

Degree of Controllability 


This appendix provides further details on the degree of controllability for discrete time 
systems. The development and terminology closely follows that for the continuous 
time case presented by Schmitendorf [42], 

B.l Recovery Region 

Consider the linear, shift-invariant, discrete time system described by: 


x n -^-i — Ax n T btin 


(B.l) 


with the scalar control bounded by the constraint: 


N < 1 i = 0, 1,2, 3,...’ 


(B.2) 
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It will be useful to define the set ft as the set of all control inputs satisfying Equa- 
tion B.2, that is, 

ft = {ui : |tt,-| < 1 t = 0, 1,2,3, ...} (B.3) 

Define 1Z( N ), the recovery region at step N as the set of all initial states, Xo which can 
be driven to the origin in N steps without violating the constraint B.2. An expression 
which can be used to generate the recovery region can be obtained as follows. Choose 
any control sequence, tt 0 , «x, . • • , which satisfies Equation B.2, and let x 0 be the 
as yet unknown initial system state. The system, Equation B.l then evolves as: 


Xi = Ax 0 + bu 0 

x 2 = A 2 x 0 + Abu 0 + bui 

x 3 = A 3 Xo + A 2 buo + Abtti + bu 2 


x N = A n Xq + A N ~ 1 bu 0 + A N 2 bu! + ■■■ + bu N _! 


So if Xjv = 0 we have: 

0 = A N Xo + A JV-1 b u 0 + A^“ 2 bui -f 1- bu^_i (B-4) 

t 

If A is invertible (For the case of interest here, where the system Equation B.l is 
obtained by discretizing a continuous time system , A is a matrix exponential and 
therefore always invertible.) then we can solve Equation B.4 for the initial state, 
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giving 

x 0 = - (A -1 tm 0 4- A -2 bui 4 1- A -iV bu^r-i) (B.5) 

So that, 

'R(N) = |xo : Xo = — (A -1 bu 0 + A _2 bu 1 4- H A ^bujv-i) , |w,| < l} (B-6) 

B.2 Attainable Set 

The attainable set at step N, , is defined as the set of all states which can be reached 
from an initial state, x 0 , in N steps, with the control bounded by the constraint, 
Equation B.2. Following a development similar to that just given for the recovery 
region, we have: 

cqv(xo) = : x = A N x 0 + A^-’bu,-.!, |u,| < l| . (B.7) 

We will make use of the fact that the attainable set is a closed, bounded, convex set 
[26]. For the subsequent development, we will also utilize the fact that 0 is an element 
of a N (x 0 ), that is the system can be driven to the origin in N steps from the initial 
condition x, if and only if the following condition holds. 

min |x[)Z 4 - max {^u,_ib / (A') _, z|| > 0. (B.8) 

||z||=il JJ 
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This is shown as follows. Define the support function on the attainable set by. 


h(x) = sup {A, x), (B.9) 

Aea^(Xo) 

for all *6R n for which the supremum is finite. 

Since the attainable set is convex, closed and bounded, it may be shown utilizing 
the separating hyperplane theorem (See [25, Theorem 4.12]) that 0 is an element of 
Ojv(xo) if and only if /i(x) > 0 for all x. 

Using the definition of the attainable set then gives the following necessary and 
sufficient condition for 0 to be in the attainable set: 

0 < sup | A'A^Xp + £ \'A n (B.10) 

u.e n ( ;=i ) 


Defining 


(AO^A 
Z “ II ( A / ) Ar A|| 


(B.ll) 


and substituting into Equation B.10 then gives the desired result. 


B.3 Degree of Controllability 

Define the degree of controllability at step N, as the largest ball centered at 

the origin 1 which is contained in 7 Z(N). 

We now find an expression for the degree of controllability. From Equation B.8, 
1 The ball of radius r centered at the origin, denoted, B r (0), is the set of points given by B r ( 0) — 

{x e R n : INI < 0}. 
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0 is an element of a n (xo) for every, Xo in a ball of radius r centered at the origin, if 


and only if: 

inf { min {x£,z + max{^ u i _ib , (A , ) - ’z}} > 0. (B-12) 

Xo€flr(o) || z ||_i tx.en fr[ 

Noting that, 

infKz:xo€£ r (0)} = -r||z|| (B.13) 

= —r (B.14) 

we have from Equation B.12, 


min | max | > ^t*,‘-ib / (A / ) 'zl > rl , (B.15) 

||z||=i \£i J J 

as a necessary and sufficient condition, for 0 6 a^Cxo) f° r every Xo €: B r { 0). Since, by 
definition, the degree of controllability is the largest value of r for which this condition 
holds we have the result: 

min (max!fn,. 1 b'(A')-zU = e (W). (B.16) 

||4=i l“i« n Is > >■ 


This expression can be further simplified to obtain: 


min {f>'(A')-’'z|} = f>W ( B - 17 ) 

IMI-dw j 


To obtain a numerical value for the degree of controllability, we must in general 
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use a numerical method to find the minimum indicated in Equation B.12. For the 
computations presented in Chapter 5, a system having two state variables was stud.ed. 
This requires a minimization of Equation B.17 around the unit circle. The search 
can be parameterized by a single variable over the closed interval, [0,2*], and the 
MATLAB routine “fmin” was used for this purpose. 
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Appendix C 


Frequency Domain Limitations 

This appendix provides a derivation for the lower bound on the Complementary 

Sensitivity function given by Equation 5.22. 

From Freudenberg and Looze, [15, Equation 2.17], for a system with no right half 

plane zeros or time delays we have: 

0 = / log \T(ju>)\dO p (u>) (C-l) 

J — CO 

Define the peak complementary sensitivity as: 

II^CMIlco = sup \T(ju)\. ( C - 2 ) 

<jJ 


Define 0 = {u : -u c < u < w c }, and let ft c denote the complement of ft. 
Suppose that for some £ > 1, \T{ju)\ < 1/e for all u 6 ft c . Let 0 = 1/e. From 
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Equation C.l we have: 


0 = J log \T(juj)\dO p (u) + J^ c log \T(ju>)\d9p(u) (C.3) 

< log(||T(ju;)|| 00 )^(fi) + \og(p)O p {n c (C-4) 

= log||T(jw)||oo Mn) +log^' (nC) ( C - 5 ) 

= log (\\T(ju)\\J’ {Q) p 6 > (QC) ) ■ (°- 6 ) 

Exponentiating both sides of Equation C.6 and rearranging provides the desired ex- 
pression: 

||T0‘ w )lloo>ei ^ ' (C ' 7) 


The derivation for the case with time delays and/or right half plane zeros is similar. 
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Appendix D 


Minimal RMS Actuator Activity 


Consider the linear time invariant system: 

x = Ax + bii + Lio, (D-l) 

where x is a vector of system states, u is the scalar control input signal, w is a scalar 
external disturbance signal, and A, b, and L are constant matrices of appropriate 
dimensions. Assume that the system is open loop unstable, ie, A has eigenvalues in 
the closed right half complex plane. Suppose that u is given by the linear feedback 
control law: 

u = Kx ( D - 2 ) 

where A is a real valued, constant matrix of appropriate dimension. Let us also 
assume that the closed loop system with feedback is stable, ie, the eigenvalues of 

A - b I< are all in the open left half complex plane. 

A measure of actuator effectiveness is given by the the minimum possible RMS 
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(root mean square) closed loop response of u to a Gaussian white noise disturbance 
w(t) with intensity V . 

To calculate this value we proceed as follows. Find the gain matrix, K , which 
minimizes the performance criteria 

Ji = J (x(t)'Rix(t) + pu 2 ) dt, (D.3) 

It is known (See [24, pages 260-264]that this same set of gains will minimize: 

J 2 = lim £{x(*)'Rix(f) + P« 2 }, ( D - 4 ) 

<—►00 

where E denotes the expectation operator. We are interested here in the asymptotic 
limit as p -> 00 . This limit will provide the solution with minimal RMS actuator 
motion and is known as the solution to the “expensive control problem. 

To find the required gains, we first utilize the results of [24, Theorem 3.11] which 
tells us the optimal closed loop pole locations. We then use a standard pole placement 
technique to find the the required gain matrix. In particular it is noted that the closed 
loop and open loop pole locations are identical for any open loop poles in the left half 
complex plane. Any unstable open loop poles are reflected about the ;w-axis into the 
left half of the complex plane to obtain the corresponding closed loop pole locations. 
Once the optimal gain, K , is known, we can solve for the RMS actuator activity 

using: 

JjkW} = y/KQK’, (D.5) 
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where, Q is the solution of the Lyapunov equation: . 

0 = (A - b K)Q + Q( A - b K)' + V. ( D - 6 ) 

(See for example [24, Equation 3-377 and 3-338] 
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Appendix E 

Lyapunov Analysis of 
Close-coupled Control 

E.l Introduction 

The analysis presented in the body of this thesis has been primarily based upon linear 
system theory. Although much valuable information and insight can and has been 
derived from the linear analysis, there are important questions for which it cannot 
provide answers. A nonlinear analysis will be presented here to consider two such 
questions: 

1. What is the extent of the allowable region of initial conditions from which 
subsequent system behavior will be acceptable? 

2. Assuming that disturbances acting on the system are in some sense bounded, 
how large a bound can be tolerated while maintaining acceptable system be- 
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havior? 


The nonlinear, Lyapunov based analysis, to be presented here 1 , shows that the an- 
swers to these questions depends to a great extent upon the overall shape of the 
compressor’s pressure-flow characteristic. This contrasts with the case of the linear 
stability analysis for which only the slope of the compressor characteristic is impor- 
tant. The fact that the overall shape of the compressor characteristic is important, 
suggest that a close-coupled resistance in series with the compressor could be used to 
modify the overall effective compressor characteristic to achieve more desirable sys- 
tem behavior. If a servo controlled valve is used for this purpose, rather than a fixed 
flow restriction, considerable freedom for shaping the characteristic can be obtained. 
The purpose of this appendix is to present some basic theory which can be employed 
to exploit this design freedom to best advantage. 


E.2 Dynamic System Model 


E.2.1 Basic Model 


The basic physical configuration to be considered consists of a compressor feeding 
a plenum volume which then discharges through a load resistance. This is shown 
schematically in Figure E-l which also shows typical non-dimensional pressure-flow 
characteristics for the compressor and load along with the assumed sign conventions. 


As indicated in Figure E-l , the typical non-dimensional pressure-flow characteristic, 

iThe material included in this appendix was originally presented by the author in a similar form 
at the 1991 American Control Conference [45], (For a closely related analysis also see Wyatt, [51J 

and Simon [44]) 
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^ c , of the compressor is non-monotonic. With the assumed sign conventions, it 
decreases monotonically in the second quadrant, has a single local minimum at a flow 
rate, and a single local maximum at a higher flow rate, <j> c2 , in the first quadrant 
and then decreases monotonically into the fourth quadrant. The flow through the 
load, given by the load characteristic, monotonically increases with increasing 

plenum pressure, V> p , and P asses onl y throu S h the first and third fl uadrants - In the 
remainder of this appendix, sufficiently smooth, continuous characteristics having 
these general shapes will be assumed. The following assumptions are also made: 

l The flow in the compressor duct is one dimensional, incompressible and inertia 
dominated (inviscid). 

2. An algebraic relationship exists between the compressor pressure rise and flow 
rate through the compressor at all times (sometimes called the quasi-steady 
compressor assumption) . 

3. The pressure in the plenum is uniform and the compression process in the 
plenum is isentropic. 

Application of the principles of conservation of momentum in the compressor duct and 
conservation of mass in the plenum then allows the system dynamics to be represented 
by a coupled set of nonlinear ordinary differential equations Following Greitzer [17], 
and choosing as state variables the non-dimensional flow rate through the compressor 
duct, <f> c , and the non-dimensional pressure in the plenum, V> p , leads to the governing 


221 


set of state equations: 


<t> c = B('Pc-V’p) 

p = ~ $<) 

Where <f> c , is the non-dimensional mass flow rate or flow coefficient, t/> p , is the non- 
dimensional plenum pressure or pressure coefficient and the parameter, B, is a non- 
dimensional ratio of plenum capacitance to compressor duct inertance, as described 

in [17]. 

E.2.2 Model with Control 

As will be shown, the stability of the basic compression system can be substantially 
modified by the introduction of a servo actuated, close-coupled control valve at the 
exit of the compressor. The term close-coupled here implies that there is no significant 
mass storage (fluid capacitance) between the valve and the compressor. The setup 
is shown schematically in Figure E-2. As indicated in this figure, the pressure drop 
through the valve is adjusted (by changing the valve area) as a function of the flow 
rate through the compressor to produce the desired functional relationship „ - With 
the introduction of the control action, ¥*, the state equations become: 

<j> c = B{V e -rJ> p-*«) ( E - 2) 

4> P = 
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In this form, the state equations explicitly exhibit the effect of the control action, 
which is to induce an acceleration of the fluid in the compressor duct. It is, however, 
convenient in the following analysis, to combine the compressor pressure rise, fl' c , and 
the control valve pressure drop \E r v , into a single function ^ e (</>), defined as. 


$„(</>) = ( E - 3 ) 

* 

The function, # e (<£), is interpreted as an equivalent compressor characteristic for the 
series combination of the compressor and control valve. This is shown graphically in 
Figure E-3. The state equations with control, Equation E.2, can then be rewritten 
using the equivalent compressor characteristic as. 

4 = BOM*)-*,) ( E - 4 ) 

V>P = - Qt) 

A pair of values, (<f > a 0 ,^ P0 ), which makes the time rate of change of the state vector, 
(the left hand side of Equation E.4) identically zero will be called an equilibrium 
operating point. In general, there may be more than one such point in the state 

space. 

Finally, it will be useful to transform the state equations into a new local coordi- 
nate system in which the origin is located at an equilibrium operating point, ( <f>co , 
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whose stability is being investigated. This is done by defining the transformations: 


4> 

— 4>c $c0 

(E.5) 


= Vv — VVo 

(E.6) 

Hi) 

= «e(^ + ^co)-*e(^o) 

(E.7) 

Hi) 

= + rp PQ ) - $(^ P0 ) 

(E.8) 


The resulting transformed state equations are: 

= - VO ( E - 9 ) 


d > 
dt 

dj>_ 

dt 


The state equations in this local coordinate system will be used in the remainder 
of this appendix. The state variables in this local coordinate system will be called 
the unsteady flow and unsteady pressure. The values of these unsteady state vari- 
ables measure deviations of the state from equilibrium, which are not assumed to be 
infinitesimally small. 
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E.3 Lyapunov Analysis 


E.3.1 Unsteady Energy and the Lyapunov Formalism 

A real, scalar, non-negative value, which will be called the unsteady energy , is assigned 
to each point in the state space by the Lyapunov function. 

= + (E.10) 


It is noted that this function is identically zero only at the origin (in the local coor- 
dinate system) and increases with increasing distance from the origin. Its graph is a 
bowl shaped surface whose level sets are ellipses enclosing the equilibrium point. The 
shape of these ellipses of constant varies with the value of the B parameter, 

however, the axes of the ellipses are always aligned with the pressure-flow axes in the 

state plane. 

The rate of change of V(i, VO along trajectories (the material derivative of V{j>, VO 
) is found by application of the chain rule to be. 

* 

V(V,^) = ^e(V)-^(W ( E - U ) 

which shows the remarkable fact that the time rate of change of unsteady energy along 
trajectories depends only upon the shape of the characteristics of the two resistive 
elements, namely, the equivalent compressor and the load. The inertial and capac- 
itive elements (through the B parameter), therefore, only determine the amount of 
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unsteady energy, V(<j>, 4>) , at any particular operating point but not the rate at which 
it changes. The two terms on the right hand side of Equation E.ll are each products 
of unsteady pressure and unsteady flow and are thus power-like. These products are 
interpreted as unsteady power flows, the first term representing the unsteady power 
production of the equivalent compressor (compressor and control valve in series) and 
the second term represents the unsteady power dissipation of the load. The unsteady 
energy of the system, V(4>, $) is increasing, that is, is positive, at any state 

where more unsteady power is produced by the equivalent compressor than can be 
dissipated by the load. Similarly, the unsteady energy is decreasing, that is, 0) 
is negative, whenever more unsteady power is dissipated in the load than is produced 
by the compressor. 

Because the compressor pressure rise does not, in general, decrease monotonically 
with increasing flow, the algebraic sign of the unsteady power produced by the com- 
pressor can vary with flow rate as shown in Figure E-4. Thus, for a typical equilibrium 
operating point, the compressor will produce unsteady power over some ranges of flow, 
thereby tending to make the sign of positive while, over other ranges, it dis- 

sipates unsteady power, thereby tending to make the sign of V{<f>,ip) negative. This 
contrasts with the case of the load resistance whose flow rate is assumed to increase 
monotonically with increasing pressure drop. The load, therefore, always dissipates 
unsteady power and so it always tends to make the sign of V{4>^) negative. This is 
also shown in Figure E-4. 
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E.3.2 Basins of Attraction and Ultimate Boundedness 

Ideally, the system state would return rapidly to equilibrium, without excessive ex- 
cursions along the way, from any initial point in the state space. As will be seen, this 
desirable property, which is called Global Asymptotic Stability , (for a precise definition 
see, for example, [41]) can be attained using the close-coupled control but only by 
incurring steady state losses across the control valve. If the required amount of loss 
is considered excessive, then some weaker form of stability must be accepted. In this 
case, all initial states will not return to the equilibrium point. The set of points that 
do return are called the basin of attraction. It will now be shown that the ellipsoidal 
level sets of unsteady energy, V(0,0), can be used to bound a subset of the overall 
basin of attraction as well as to bound the maximum state excursion which the state 
will attain as it approaches equilibrium. 

Suppose there is some open set, R u of points including the equilibrium point for 
which V{4>, tf) < 0 with equality holding only at the equilibrium point. Then, for 
a sufficiently small value of the unsteady energy, call it Vi, the ellipse of constant 
given by V(4>J) = V x will be entirely contained in the region R x . Since 
V(<i>, Tp) is negative for all points inside of this ellipse (Except for the origin where it is 
zero), the value of V{j>, $) can only decrease along the trajectories of any point which 

originate in the interior of the ellipse V{M) = * • Thus ’ P oint ori S inatin S 
inside of this ellipsoidal region can never leave it, and must eventually reach the 
equilibrium point as the value of V(M) decreases inexorably to zero. This ellipse 
of constant V$A) > therefore, provides both a conservative estimate of the basin of 
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attraction and as a bound on the maximum state excursion which will be experienced 
along any trajectory originating in its interior. This estimate is termed conservative, 
because there may be some initial states outside of this boundary which nevertheless 
return to equilibrium. It is noted that if there are no points in the state space for 
which V((j>, t/>) is positive then the bounding ellipse may be made arbitrarily large, 
that is, the basin of attraction will be the entire state space, which implies that the 

equilibrium point is globally asymptotically stable. 

It is also useful to be able to obtain some bounds on the transient excursions as 
well as the ultimate asymptotic behavior along trajectories which originate outside 
of the previously estimated basin of attraction. It will now be shown that level sets 
of the unsteady energy, F(^,V>), can also be used for this purpose. 

Suppose that a bounded set, R 2 , can be found which contains all the points in 
the state space for which V(M) is positive. Then, for a sufficiently large value, V 2 , 
of the unsteady energy, the entire region, R 2 can be contained in the interior of an 
ellipse defined by V(j>, $) = V 2 . For brevity, this bounding ellipse will be called 
V 2 . By the above construction, will be strictly negative outside of V 2 . The 

initial unsteady energy, call it V 3 , along any trajectory originating outside of this 
ellipse must therefore be decreasing and so the trajectory can never again leave the 

interior of the region bounded by the ellipse = V 3 determined by the initial 

> 

unsteady energy. Furthermore, the unsteady energy will continue to decrease until 
the state enters the interior of the ellipse V 2 . Once inside of the boundary, V 2 , the 
trajectory can never again leave. To see this, note that for the system’s state to leave 
the region bounded by V 2 its unsteady energy would have to exceed the value V 2 but, 
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since outside of the region bounded by the ellipse V 2 the unsteady energy is strictly 
decreasing, this would be impossible. A similar argument shows the impossibility of 
the trajectory ever leaving the region bounded by the ellipse V = K 3 set by the initial 
unsteady energy. Thus it can be concluded that: 

1. Any trajectory originating outside of the ellipse F 2 will eventually be inside it. 
This is ultimate boundedness (See, for example, [41]). 

2. The maximum excursion along any trajectory originating outside of the ellipse 
V 2 is bounded by the ellipse V$J) = V 3 determined by the initial unsteady 

energy, V 3 . 

3. Once inside of the region bounded by the ellipse V 2 the state can never again 
leave this region. 

4. Trajectories originating inside of the region bounded by the ellipse V 2 can also 
never leave this region 

The bounds developed in this section are illustrated schematically in Figure E-5. In 
this figure, there is a single isolated region in state space where V(jJ) is positive. In 
this region, the compressor is producing more unsteady power than can be dissipated 

by the load. 

E.3.3 Effect of B parameter 

As previously noted, the value of the B parameter does not affect the rate at which 
unsteady power is produced at a given location in state space. It does however, change 
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the shape of the level sets of and, thereby, changes the previously described 

ellipsoidal bounds on the transient excursion and estimates of the basin of attraction. 

The effect of the value of the B parameter on the shape of the ellipsoidal level sets 
of constant V{i, i>) can be determined by inspection of Equation E.10. Three distinct 
regimes can be identified corresponding to “large”, “intermediate and small values 
of the B parameter as follows: 

Large B For B > 1 the major axis of the ellipsoidal level sets lies along the unsteady 
flow axis. 

Intermediate B For B = 1 the level sets are circular. 

Small B For B < 1 the major axis of the ellipsoidal level sets lies along the unsteady 
pressure axis. 

Thus, as shown schematically in Figure E-6, for the same equilibrium operating point 
and resistive element characteristics, the bounds on the possible system behavior may 
vary markedly with changes in the value of B , even though the regions of unsteady 
power production remain the same. In particular, it is noted that as the B parameter 
increases, the estimated bound on the basin of attraction shrinks, indicating that the 
acceptable region of initial conditions will become increasing difficult to maintain as 
B becomes excessively large. 
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E.4 Regions of Positive V (</>, i/>) 

In order to find the bounds developed in the previous section, it is necessary to 
first determine the regions in state space where there is a net production of unsteady 
power, that is, where V(tJ) is positive. Further, in order to make this theory useful, 
it must be determined how the control action of the close coupled valve can be used to 
modify the location and extent of these regions of unsteady power production. This 

will be the subject of this section. 

E.4.1 General Features 

In general, for the compressor and load characteristics which are assumed here, regions 
of unsteady energy production will be be bounded, open sets, containing portions of 
the unsteady flow axis. In fact, V(W) will be maximal along the unsteady flow axis. 

These facts will now be justified. First of all, it was assumed at the outset that 
the compressor characteristic was continuous and, for sufficiently large positive or 
negative values of the flow through the compressor, <£, the characteristic falls in the 
second and fourth quadrant respectively. That is, for sufficiently large magnitudes of 
l the compressor dissipates unsteady power. Thus, from the sum expressed by Equa- 
tion E.ll, it can be seen that regions of positive V(jj) can only exist over a finite 
interval of unsteady flow values. Further, from the continuity of the compressor char- 
acteristic, it can be seen that that the unsteady power produced by the compressor 
is bounded. On the other hand, because of the assumed monotonicity and quadrant 
constraints on the load characteristic, the unsteady power dissipated by the load m- 
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creases without bound as the magnitude of the unsteady plenum pressure increases 
(in either a positive or negative sense). Therefore, for sufficiently large magnitudes of 
unsteady pressure, the bounded unsteady power production of the compressor will be 
dominated by the dissipation of the load resulting in a net negative value of ip) . 
Thus, it can be concluded that the possible regions of positive V{<j>A) are bounded. 
Since any nonzero value of the unsteady plenum pressure, 0, will cause dissipation 
by the load, the maximum value of V{j>A) must occur along the unsteady flow axis 
(that is the set of point in the state space where ip = 0 ). It can, therefore, be 
concluded that any regions of positive V{}, j>) must include portions of the unsteady 
flow axis and it will be maximal there. 


E.5 Effect of Bounded Disturbances 

In this section, the effect of persistent disturbances on the compression system will 
be considered. Up to this point, the compression system has been treated as if there 
were no persistent disturbances driving it; only the effect of initial conditions has been 
considered. An actual system will certainly be disturbed by local (sub-model scale) 
flow instabilities in the compressor and load flow as well as by external disturbances 
at the inlet and exit of the system. For acceptable operation, the compression system 
must be able to remain reasonably close to equilibrium in the face of these persistent 
disturbances. In this section, it will be assumed that such disturbance are bounded 
in magnitude, which is thought to be a reasonable physical assumption. No further 
assumption regarding the statistical nature or time structure of these disturbances will 
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be made. They could range from constant biases, reflecting some uncertainty in the 
compressor characteristic, to rapidly fluctuating forces due to turbulence in the flow 
through the compressor. Such a disturbance model is called an unknown but bounded 
model (see [43]). (Actually, to be completely rigorous, some additional restrictions 
probably must be made on the allowable class of disturbance functions in order to 
assure that the solutions to the differential equation representing the modeled system 
are well defined. It will be assumed here that the class of nonpathological functions, 
(ones for which the formulation used here is valid) is large enough to motivate the 
effort of studying them. The exact technical restrictions which are required can be 
investigated at some future time as deemed necessary.) It will now be shown that the 
previously defined unsteady energy, V(W), can be used to obtain useful bounds on 
the excursions from equilibrium which result from these bounded disturbances. 

E.5.1 Analytical Formulation 

Consider the effect of a bounded pressure disturbance, V d {t), tending to accelerate 
the flow in the compressor duct and a bounded flow disturbance, modeling 

unsteady plenum outflow. The previously derived state equations in the local (un- 
steady) coordinate system, Equation E.9, become with the addition of these distur- 

bance terms: 


d t 

ii = I(^-$(^)-‘M o) 

dt t> 
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The set of equations E.12 will be referred to as the disturbed state equations. Using 
the unsteady energy, as previously defined by Equation E.10, the chain rule 

is used compute the material derivative as in Section E.3.1. However, this expression 
is evaluated using the disturbed state equations E.12 to obtain the formula. 

V(i, tj>) = $*(]>) - + ^d(t) - ( E - 13 ) 

Comparison of Equation E.13 with the previously derived expression given for V(j>, ^), 
Equation E.ll, shows that the first two terms on the right hand side of both equa- 
tions are identical. The additional two terms on the right hand side of Equation E.13 
are unsteady power flows from the disturbances. They can be seen to have the effect 
of causing fluctuations in the overall unsteady energy balance. Thus, the regions in 
state space where V{d>A) is positive, now vary with time. Because the disturbances 
are bounded, the possible regions in the state space where V{hl>) can be positive 
also remain bounded. The boundaries of these regions however fluctuate with time. 
In particular, points in a neighborhood of the origin, where in the undisturbed case 
V(d>, was non-positive but close to zero, may now occasionally become positive. 
Because the algebraic sign of the disturbances may be either positive or negative, the 
disturbances may be either dissipating or producing unsteady power. As illustrated in 
Figure E-7, the envelope of these fluctuations gives a best case and worst case bound 
on the regions in state space where V(j>, rf>) may be positive, where the best case is the 
smallest possible region of positive and the worst case is the largest possible 

region. For a conservative design the larger (worst case) region must be used. 
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Although the magnitude of the disturbances is bounded, the unsteady power prod- 
ucts involving them are not. It, therefore, might not be clear at this point that these 
regions of positive V{<j>, $) remain bounded. To clarify this point, the unsteady power 
balance, Equation E.13 is rearranged to obtain. 

V (<j>, j>) = <}>[$'(]>) + ^ d(t )] - VWVO + MO] ( E - 14 ) 

It was assumed at the outset that |* e (#)I increases without bound as \}\ - oo . 
However, the disturbance term, V d (t) is bounded. Thus, for sufficiently large abso- 
lute values of $ e (#) the sum ’ + MO ’ wiU haVe th<? ^ algebraic sign &S 

* e (#). A similar argument shows that for sufficiently large absolute values of the 

sum $(i j>) + $ d (t) , will have the same algebraic sign as $(# This implies that, for 
sufficiently large distances from equilibrium, the overall sum, Equation E.14 repre- 
senting V{M) for the disturbed system, will have the same algebraic sign as that 
of the undisturbed system given by Equation E.ll. In particular, if beyond some 
distance from equilibrium the value of V&i) for the undisturbed system is nega- 
tive, then the same will hold for the disturbed system. Thus, the boundedness of the 
regions of positive for the undisturbed case, which was established earlier in 

Section E.4.1, implies the boundedness of these regions for the disturbed case. 

From Equation E.14, the disturbances can be seen as having the effect of putting 
a bounded range of uncertainty on the compressor and load characteristics. This 
interpretation is illustrated in Figure E-8 which shows the nominal characteristics 
along with the envelope of the bounded disturbances. These upper and lower bounds 
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will be denoted by ^P c (^) and ^ c ( (^) respectively where these terms are defined by. 

¥ c (^) = sup(W) + $ d (t)) ( E - 15 ) 

<>0 

i(i) = jnf («(#) + ^ rf (0) 

As indicated in this figure, this envelope on the characteristics, may cross into a 
quadrant of unsteady power production, while the nominal characteristic is dissipating 
unsteady power. This shows graphically that the possible regions of unsteady power 
productions will become larger as a result of the disturbances. 
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Figure E-l: Basic Compression System Configuration 



Figure E-2: System with Close-Coupled Control Valve 
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Compressor: 
active with 
with respect to 
operating point 



$ 





Figure E-4: Activity and Passivity with respect to Operating Point 
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Figure E-5: Constant Unsteady Energy - Ellipsoidal Bounds 
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Figure E-8: Envelope of Bounded Disturbance on Compressor Characteristic 
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